CERN-PH-TH/2004-140 



Theoretical tools for CMB physics 



Massimo Giovannini 



Department of Physics, Theory Division, CERN, 1211 Geneva 23, Switzerland 

Abstract 

This review presents, in a self-consistent manner, those analytical tools that are 
relevant for the analysis of the physics of CMB anisotropics generated in different the- 
oretical models of the early Universe. After introducing the physical foundations of the 
Sachs- Wolfe effect, the origin and evolution of the scalar, tensor and vector modes of the 
geometry is treated in both gauge-invariant and gauge-dependent descriptions. Some of 
the recent progresses in the theory of cosmological perturbations are scrutinized with 
particular attention to their implications for the adiabatic and isocurvature paradigms, 
whose description is reviewed both within conventional fluid approaches and within 
the the Einstein-Boltzmann treatment. Open problems and theoretical challenges for 
a unified theory of the early Universe are outlined in light of their implications for the 
generation of large-scale anisotropics in the CMB sky and in light of the generation of 
stochastic backgrounds of relic gravitons between few Hz and the GHz. 
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1 An inhomogeneous Universe 



1.1 Formulation of the problem 

What is the origin of the inhomogeneities we observe in the sky? What is their typical 
wavelength? How come that we are able to observe these inhomogeneities? The term in- 
homogeneity is rather generic and it indicates, for the purpose of the present introduction, 
fluctuations both in the background geometry and in the energetic content of the Universe. 
For simplicity, these inhomogeneities can be represented, for the purposes of the present in- 
troduction, by plane waves characterized by a comoving wave-number k. Since the evolution 
of the Universe is characterized by a scale factor a{t), a physical wave-number u = k/a is 
also customarily defined. The comoving wave-number does not feel the expansion while the 
physical momentum is different at different epochs in the life of the Universe. Conversely 
the value of a given physical frequency is fully specified only by stating the time at which 
the physical frequency is "measured" . 

If a given fluctuation has a momentum comparable with the present value of the Hubble 
parameter ^, i.e. cjq — -f^o — 5.6 x 10~^^ Mp, we will have that — 2.3 x 10~^^ Hz. 
Fluctuations with momentum smaller that Hq have a wave-length which is larger than the 
present value of the Hubble radius and, therefore, they seem to be impossible to detect 
directly since the distance between two maxima (or two minima) of the wave will be, in this 
case, larger than our observable Universe. 

At the decoupling epoch the Hubble rate is -ffdec — 6.7 x W~^^Mp. The decoupling fre- 
quency red-shifted today will be i^dcc — 3 x 10~^^Hz. Fluctuations of this typical frequency 
can be "detected" through the study of the of the temperature fluctuations of the Cosmic 
Microwave Background (CMB). The rationale for this statement is that the fluctuations in 
the background geometry and in the density contrasts of the various species composing the 
primeval plasma may induce tiny spatial variations in the CMB temperature. The standard 
lore for the formation of the structures in the Universe, is that the fluctuations detected 
by means of CMB temperature inhomogeneities will eventually collapse by gravitational 
instability to form galaxies and clusters of galaxies. 

At the time of decoupling and, anyway prior to the recombination of free electrons 
with protons, the temperature of the Universe was of the order of a fraction of the cV. 
The Universe, all along its history, reached much higher temperature (or, equivalently, 
much higher values of the expansion rate). At the epoch of the formation of light nuclei the 
temperature of the Universe was of the order of a fraction of the MeV and the corresponding 
value of the Hubble expansion rate was of the order of 8 x 10~^^ Mp. When (and if) the 
electroweak (EW) phase transition took place, the temperature of the Universe was of the 
order of 100 GeV leading, approximately, to Hew ^ 10~^^Mp. Comoving wave-numbers of 
the order of the Hubble rate at the EW or at the BBN epoch correspond, today, to physical 
frequencies of the order of 6<Jcw ~ 2x 10~^ Hz and — 9x 10~^^ Hz. Finally, in the context 
of the conventional inflationary paradigm the value of the Hubble rate during inflation was 
i?inf ~ 10~^Mp. Under the assumption the Universe was dominated by radiation right after 

^The definition of the Planck mass adopted here is Mp = (SttG)"^/^ = 2.4 x 10^® GeV 
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inflation, the present value of the correspoonding physical frequency was of the order of 10 
Hz. 

To "observe" fluctuations corresponding to present frequencies 10~^ Hz < w < MHz, 
CMB anisotropy experiments cannot be used. CMB anisotropy experiments are sensitive to 
present frequency scale only slightly higher than 10^^^ Hz. For frequencies in the interval 
10~^ Hz < u; < MHz, the only hope of direct observations are the experiments on direct 
detection of stochastic backgrounds of gravitational radiation. In some specific configura- 
tions these instruments could detect not only tensor fluctuations of the geometry but also 
scalar and vector modes. 

The cosmological inhomogeneities are usually characterized by their correlation func- 
tions. The simplest correlation function encoding informations on the nature of the fluctu- 
ations is the two-point function whose Fourier transform is usually called power spectrum. 
The two point function is computed by taking averages of the fluctuation amplitude at 
two spatially separated points but at the same time. From this quantity, with appropriate 
mathematical manipulations, it is also possible to deduce other relevant physical quantities 
like, for instance, the energy density carried by these fluctuations. 

From the physical scales discussed above, one is led to conclude that the power spectrum 
of cosmological inhomogeneities is defined over a huge interval of frequencies. For instance, 
the tensor fluctuations of the geometry which only couple to the curvature and not to the 
matter sources of the Universe have a spectrum ranging, in some specific models, from 10~^^ 
Hz up to the GHz for, roughly 27 powers of ten. 

From the figures obtained so far, legitimate questions may be raised: 

• given that the precise thermodynamical history of the Universe above 10 MeV is 
unknown, how it is possible to have a reliable framework describing the evolution of 
the cosmological fluctuations? 

• does it make sense to apply GR also up to curvature scales which are almost Planckian? 

• how is it possible to give sensible initial conditions to the various fluctuations of the 
geometry without having some handle on the theory of the initial singularity, i.e. the 
theory of the big-bang? 

We do not have any (even indirect) test of the thermodynamical state of the Universe for 
temperatures larger than the MeV. The four light isotopes D, ^He, ^He and ^Li are mainly 
produced at the big-bang nucleosynthesis below a typical temperature of 0.8 MeV when 
neutrinos decouple from the plasma and the neutron abundance evolves via free neutron 
decay. The abundances calculated in the simplest (homogeneous and isotropic) big-bang 
nucleosythesis model agree fairly well with the astronomical observations. This is the last 
indirect test of the thermodynamical state of the Universe. For temperature larger than 10 
MeV, collider physics and the success of the electroweak standard model offer a precious 
handle up to the epoch of the electroweak phase transition occurring for temperatures of 
the order of 100 GeV. There is then the hope that our knowledge of particle physics could 
help us to fill the gap between the nucleosynthesis epoch and the physics beyond above the 
electroweak curvature scale. 
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The evolution of inhomogeneities over cosmological distances may well be unaffected by 
our ignorance of the specific history of the expansion rate. The key concept is the one of 
conservation law. In the framework of Einsteinian theories of gravity (and in particular 
in GR) its is possible to show, on a rather general ground, that there exist conservation 
laws for the evolution of the fluctuations. These conservation laws are exact in the limit 
of vanishing comoving frequency, i.e. A; — > 0. This conclusion holds under the assuming 
that general relativity is a valid description throughout the whole evolution of the Universe. 
Similar properties also hold, with some modifications, in scalar-tensor theories of gravity (of 
Brans-Dicke type) like the ones suggested by the low-energy limit of superstring theory. The 
conservation laws of cosmological perturbations apply strictly in the case of the tensor modes 
of the geometry. In the case of the scalar fluctuations, the precise form of the conservation 
law depends upon the matter content of the early Universe since scalar fluctuations couple 
directly to the (scalar) matter sources. There is the hope that experiments aimed at the 
direct detection of stochastic backgrounds of relic gravitons will provide informations on 
the high-frequency behaviour of cosmological fluctuations. These informations could be, in 
principle, converted into valuable clues on the early evolution of the Hubble rate. 

Finally, there is no way of addressing the singularity problem in general terms. For 
instance, the conventional inflationary paradigm postulates that the evolution of the early 
Universe was driven by the potential energy of a single (scalar) degree of freedom named 
the inflaton. During inflation, the scale factor expands in an accelerated way (i.e. a > 
and a > 0). Inflationary cosmology does "predict" a singularity in the far past of the 
Universe, however, the effects of the singularity may well be invisible since the accelerated 
expansion of a quasi-de Sitter space-time is able to dilute all the spatial gradients which 
may eventually arise at the time of the cosmological singularity, prior to the beginning of 
inflation. In conventional inflationary models the singularity is not addressed but it is just 
removed beyond our observational capabilities. 

As a general remark, one can say that, in modern cosmology, when the description of 
the early Universe is approached, many authors like to use a simplicity principle. This 
means, for instance, that instead of introducing many degrees of freedom in the description 
of the early Universe, we are often biased towards the model containing the fewest number 
of degrees of freedom (like in the case of single-field inflationary models). This is, however, 
not the typical situation arising in the context of the low-energy string effective action where 
various massless modes are simultaneously present. It is here useful to recall the opinion 
of R. Feynman questioning the notion that simplicity should be a guiding principle in the 
search for truth about Nature pQ: "...the simplest solution, by far, would be nothing, that 
there should be nothing at all in the Universe. Nature is more inventive than that, so I 
refuse to go along thinking that it always has to be simple." 

1.2 A short introduction to a long history 

The history of the studies on the CMB anisotropies and on the cosmological fluctuations 
is closely linked to the the history of the standard cosmological model (see, for instance, 
the the textbooks by Weinberg t2^ and by Kolb and Turner j2j). In this introduction only 
few interesting points will be swiftly mentioned since the main scope of this review is not 
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historical. Furthermore, the second chapter of the excehent book j4, of B. Patridge describes 
in detail the excitement of experimental cosmology in the early sixties at the time when 
CMB was firstly discovered. In the he issue number 81 of the "Uspekhi Fizicheskikh Nauk" 
, on the occasion of the seventy-fifth anniversary of the birth of A. A. Friedmann, a number 
of rather interesting papers were published. Among them there is a review article of the 
development of Friedmannian cosmology by Ya. B. Zeldovich 5 and the inspiring paper of 
Lifshitz and Khalatnikov [S] on the relativistic treatment of cosmological perturbations. 

Reference [3] describes mainly Priedmann's contributions [7j. Due attention should also 
be paid to the work of G. Lemaitre 13121 that was also partially motivated by the debate with 
A. Eddington jJUj. According to the idea of Eddington the world evolved from an Einstein 
static Universe and so developed "infinitely slowly from a primitive uniform distribution 
in unstable equilibrium" jlUj . The point of view of Lemaitre was, in a sense, more radical 
since he suggested, in 1931, that the expansion really did start with the beginning of the 
entire Universe. Unlike the Universe of some modern big-bang cosmologies, the description 
of Lemaitre did not evolve from a true singularity but from a material pre-Universe, what 
Lemaitre liked to call "primeval atom" 19^. The primeval atom was a unique atom whose 
atomic weight was the total mass of the Universe. This highly unstable atom would have 
experienced some type of fission and would have divided into smaller and smaller atoms by 
some kind of super-radioactive processes. The perspective of Lemaitre was that the early 
expansion of the Universe could be a well defined object of study for natural sciences even 
in the absence of a proper understanding of the initial singularity. 

During the development of the standard cosmological model, a recurrent theme has 
been the explanation of structures we observe in the Universe. This problem goes under 
the name of the structure formation problem. The reviews of Kodama and Sasaki and 
of Mukhanov, Feldman and Brandenberger contain some historical surveys that are a 
valuable introduction to this subject. 

Among the various ideas put forward through the years we wish to mention two contri- 
butions that are suitable in order to introduce some of the problems treated in the present 
article: the classic paper of Sachs and Wolfe on the possible generation of anisotropics in 
the CMB ^S] and a paper of Sakharov ^[4, where the concept of Sakharov (i.e. acoustic) 
oscillations has been introduced. 

In jl3j the basic general relativistic effects leading to the anisotropy of the CMB were 
estimated assuming that the CMB itself had a cosmological origin. This effect is now one of 
the cornerstones of the physics of CMB anisotropics and, in light of its importance, it will 
be one of the first derivations of the present review. In the authors argued that galaxy 
formation would imply AT/T ~ 10^^. However, subsequent calculations ITT H [T7 | UK] 
brought the estimate down to 10~^, closer to what experimentally observed. 

The paper of Sakharov jl4j is rather intriguing if one thinks that it was written in 1965, 
i.e. just at the time when the possible cosmological origin of CMB was recognized. Sakharov 
was speculating on a model proposed earlier by Zeldovich. The origin and details of the 
model are, today, rather unimportant. Reference [TJ contained two ideas whose legacy, 
somehow unexpectedly, survived until the present day. The first idea is summarized in the 
abstract of Ref. |14j : " A hypothesis of the creation of astronomical bodies as a result of 
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gravitational instability of the expanding Universe is investigated. It is assumed that the 
initial inhomogeneities arise as a result of quantum fluctuations..." . The model of Zeldovich 
indeed assumed that the initial state of the baryon-lepton fluid was cold. So the premises 
were wrong, but, still, it is amusing to notice that today quantum fluctuations are what it is 
assumed as initial conditions in a variety of models present at a totally different curvature 
(or energy) scale as initial condition for the fluctuations of the geometry. 

In a theoretical perspective, Refs. and [2] underline three important problems: 

• the problem of initial conditions of the fluctuations of the gravitational inhomo- 
geneities (sections 4 and 6); 

• the problem of their evolution (sections 5 and 7); 

• the problem of their imprint on the CMB temperature fluctuations (sections 3, 8). 

The sections quoted in connection with each of the previous items indicate where each of 
the mentioned topic is treated in the present article. Due to the historical development of 
the subject some semantic ambiguities may also arise. For instance, the "initial conditions 
of the CMB anisotropies" could be described either in terms of the amplitude of metric 
fluctuations around the time of matter-radiation equality (i.e. in the late Universe), or in 
terms of the normalization of the cosmological perturbations in the early Universe. 

1.3 CMB experiments 

Even if the present review is theoretically oriented, it is appropriate to give a swift account of 
our recent experimental knowledge of CMB anisotropies. A new twist in the study of CMB 
anisotropies came from the observations of the COBE satellite and, more precisely, from 
the DMR (Differential Microwave Radiometer) instrument. The DMR was able to probe 
the angular power spectrum^ Q spectrum (see section 8) up to £ ~ 25 (see, for instance, 
|19| I2nj and references therein). As the name says, DMR was a differential instrument 
measuring temperature differences in the microwave sky. The angular separation of the 
two horns of the antenna is related to the maximal multipole probed in the sky according 
to the approximate relation t9 ~ vr/^. This will give the angular separation in radians that 
can be easily related to the angular separation in degrees. 

The angular separation explored by the COBE experiment were i}^?^. After the COBE 
mission, various experiments attempted the exploration of smaller angular separation, i. e. 
larger multipoles. A definite convincing evidence of the existence and location of the first 
peak in the Ci spectrum came from the Boomerang [21\ |221 , Maxima [211 Dasi |23| 
experiments. Both Boomerang and Maxima were balloon borne (bolometric) experiments. 
Dasi was a ground based interferometer. The data points of these last three experiments 

^While the precise definition of angular power spectrum wili be given in section 8, fiere it suffices to recaii 
tliat £{£+ l)Cf/(27r) measures tiie degree of inliomogeneity in tfie temperature distribution per iogaritlimic 
interval of £. Consequentiy, a given muitipoie £ can be related to a given spatial structure in the microwave 
sky: small £ will correspond to low wavenumbers, high £ will correspond to larger wave-numbers but always 
around a present physical frequency of the order of 10~^^ Hz. 
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Figure 1: Some of the most recent CMB anisotropy data are reported (figure adapted 
from ISH]): WMAP data (filled circles); VSA data (shaded circles) [2H]; CBI data (squares) 
[231301; ACBAR data (triangles) 

explored multipoles up to 1000, determining the first acoustic oscillation (in the jargon 
the first Doppler peak) for i ~ 220. Another important balloon borne experiments was 
Archeops [2S1 providing interesting data for the region characterizing the first rise of the Ci 
spectrum. Some other useful references on earlier CMB experiments can be found in |26j . 

The Ci spectrum, as measured by different recent experiments is reported in Fig. ^ 
(adapted from Ref. [2H1)- The WMAP (Wilkinson Microwave Anisotropy Probe) data 
(filled circles in Fig. ^ provided, among other important evidences: the precise determi- 
nation of the position of the first peak (i.e. i = 220.1 it 0.8 27 ); the clear evidence of the 
second peak; a perspective (once all the four-years data will be collected) for good resolu- 
tion up to ^ ~ 1000, i.e. around the third peak. The WMAP experiment also measured 
temperature-polarization correlations (see section 8) providing a distinctive signature (the 
so-called anticorrelation peak in the temperature-polarization power spectrum for I ~ 150) 
of primordial adiabatic fluctuations (see sections 4, 8). In Fig. ^on top of the WMAP data 
we illustrated the results of the Arcminute Cosmology Bolometer Array Receiver (ACBAR) 
as well as the data points of the Very Small Array (VSA) and of the Cosmic Background 
Imager (CBI). 

To have a more detailed picture of the evolution and relevance of CMB experiments we 
refer the reader to Ref. j37j (for review of the pre-1994 status of the art) and Ref. [HSI for 
a review of the pre-2002 situation) . The rather broad set of lectures by Bond jSS] may also 
be usefully consulted. 

In recent years, thanks to combined observations of Cosmic Microwave background 
(CMB) anisotropics Large scale structure [231 1311) supernovae of type la, big-bang 
nucleosyntheis some kind of paradigm for the evolution of the late time (or even 

present) Universe emerged. It is normally called by practitioners ACDM model or even, 
sometimes, "concordance model". The terminology of ACDM refers to the fact that, in this 
model, the dominant (present) component of the energy density of the Universe is given 
by a cosmological constant A and a fluid of cold dark matter particles interacting only 
gravitationally with the other (known) particle species such as baryons, leptons, photons. 
According to this paradigm, our understanding of the Universe can be summarized in two 
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sets of cosmological parameters: the first set of parameters refers to the homogeneous 
background, the second set of parameters to the inhomogeneities. As far as the homogeneous 
background is concerned, on top of the indetcrmination on the (present) Hubble expansion 
rate^, i.e. h, there are various other parameters such as: the dark and CDM energy densities 
in critical units, i.e. respectively i^cdm — 0.224 and Oa — 0.73; the radiation and baryon 
energy density in critical units, i.e. respectively Jlr — 8 x 10~^ and Jib — 0.046; the number 
of neutrino species and their masses; the possible contribution of the spatial curvature 
(ftii = in the concordance model) the optical depth of the plasma at recombination (of 
the order of 0.16 in the concordance model)... 

The ellipses stand for various other parameters which are normally either not included 
in the analysis or just set to a fiducial value such as the equation of state for the dark 
energy component, the possible effect of the dark energy sound speed; the possible variation 
(with the rcdshift) of the fine structure constant; some possible extra (warm) dark matter 
component and so on. The second set of parameters refers to the inhomogeneities. In the 
ACDM model the initial conditions for the inhomogeneities belong to a specific class of 
scalar fluctuations of the geometry that arc called adiabatic and are characterized by the 
Fourier transform of their two-point function, i.e. the power spectrum which is usually 
parametrized in terms of an amplitude and of a spectral slope. To these two parameters 
one usually adds the ratio between the amplitudes of the scalar and tensor modes of the 
geometry. Depending on how we count the minimal list includes from 10 to 12 parameters. 
Furthermore, few optional choices raise easily the number of parameters to 16 or even 18. 

Having said this it is important to stress that the present article will not deal with the 
problem of data analysis (or parameter extraction from the CMB data) . The purpose of the 
present review, as underlined before in this introduction, will be to describe the theoretical 
tools allowing the calculation of large-scale temperature fluctuations in a broad range of 
cosmological models whose predictions may not always fall in the concordance model. 

Finally a word of care concerning references. In reviewing this subject attention has 
been given to ideas emerged in the last ten years. For this reasons, various important 
works will only be indirectly quoted and the reader may refer to earlier review articles on 
the subject that have been already quoted and will be quoted again when needed. Among 
the theoretical ideas reported here, for reasons of space two important topics have been 
omitted, i.e. the problem of cosmological inhomogeneities induced by cosmic defects and 
the problem of the (late time) fluctuations induced by inhomogeneous quintessence models. 

* Recall that the present value of the Hubble expansion rate is customarily parametrized as Ho = h x 
100km/[sec x Mpc]. In the concordance model h ~ 0.72. 
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2 FRW Universes and their inhomogeneities 



2.1 Priedmann- Robertson- Walker Universes 

A Priedmann-Robertson- Walker (FRW) metric can be written in terms of the conformal 
time parametrization, leading to the line element 



g dx^dx'^ = a'^{T)\ dr 



(2.1) 



where k can take values 1, and —1 corresponding, respectively, to spherical, eucledian or 
hyperbolic spatial sections: r is called conformal time coordinate. Throughout this paper 
the signature of the metric will be, consistently, mostly minus. Therefore, the flat Minkowski 
metric will be ry^i, = diag(+, — , — , — ). 

In the standard cosmological model the evolution of the geometry is driven by the 
evolution of the matter sources according to Einstein equations 

B!'^-\51R = SitGT;, (2.2) 

where R^y and R are the Ricci tensor and Ricci scalar; T^y is the energy- momentum tensor 
of the sources. The form of the energy-momentum tensor may vary depending on the 
evolutionary stage of the Universe. 

After electroweak interactions have fallen out of thermal equilibrium (i.e. temperatures 
T ^ 1 MeV), the energy- momentum tensor of the plasma contains photons, baryons, elec- 
trons, neutrinos, cold dark matter particles. If the temperature of the plasma is around the 
eV, then the electrons and baryons will be in thermal equlibrium at a common tempera- 
ture and photons will be tightly coupled with the baryons because of Thompson scattering. 
Neutrinos will be collisionless. 

To these species, a dark energy component is also usually added. The dark energy is 
customarily parametrized either in terms of a cosmological constant A j4Ul I41j or in terms 
of some scalar degrees of freedom generically named quintessence (see for instance |42j for 
a review). As the CDM particles, the dark energy fluid is supposed to have negligible 
non- gravitational interactions with baryons, leptons and photons. 

For temperatures eV < T < MeV the Friedmann equations are 

W^ + K = ^aVE^A + ^a^ (2.3) 
-n' + K = 4TTGa^pJ2^x{l + wx), (2.4) 

A 

p'x + 3W(1 + wx)px = 0, (2.5) 

where the summation index A refers to each component of the plasma and wx = px/p\ is 
the corresponding barotropic index. Eqs. 1)2. 3(1 and 1)2.4(1 depend on Qxi i-e- the energy 
density, in critical units, for each component of the fluid. In Eqs. ()2.3|) - (|2.5|) TC = (In a)' 
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and the prime denotes a derivation with respect to the conformal time coordinate r. The 
cosmic time parametrization t is related to the conformal one by dt = a{T)dT. In the 
i-parametrization, H = a/a \s the Hubble expansion rate; the dot denotes a derivation 
with respect to t. The relation between the Hubble factors in the two parametrizations is 
Ha = 7i. In the present article we will switch, when needed, from one parametrization to 
the other. Using the relations between the two parametrizations, Eqs. p.3|l ~ (|2.5j) can be 
written in terms of the cosmic time coordinate: 

A 

ij = -47rGp^J^A(l + w^A) + 4' (2-7) 

A " 

px + ^H{l + wx)px = G, (2.8) 

where, now, p, is the total energy density. Equations (|2.H|) - (|2.5|) assume that, for each 
species, the energy-momentum tensor can be written as 

Tiu) = (Px + PxH^^4^^ - PxQi... (2.9) 

with the normalization condition = 1 valid for each species A. Since the ther- 

modynamical history of the Universe is progressively less understood as we go back in time, 
also the energy-momentum tensor of the plasma may take different forms. 

2.2 Fluctuations of the geometry in FRW Universes 

Given a conformally flat (k = 0) background metric of FRW type 

g^,{T) = a\T)rj^,, (2.10) 

its first-order fluctuations can be written as 

Sg^lu{T,x) = 5sg^,u{T,x) + 5^g^y{T,x) + (5t5^^(r, x), (2.11) 

where the subscripts define, respectively, the scalar, vector and tensor perturbations clas- 
sified according to rotations in the three-dimensional Euclidean sub-manifold. Being a 
symmetric rank-two tensor in four-dimensions, the perturbed metric dgfj^y has, overall, 10 
independent components whose explicit form will be parametrized as 

(5(700 = 2a2</., (2.12) 
5g,j = 2a\'ilj6ij - diOjE) - a^ij + a\d^Wj + djWi), (2.13) 
6goi = -a^d^B - a^Qi, (2.14) 

together with the conditions 

diQ' = diW' = 0, hi = dih] = 0. (2.15) 

The decomposition expressed by Eqs. H2.12() - (|2.14|) and ()2.15|) is the one normally employed 
in the Bardeen formalism (see also il2j). 
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According to Eqs. ()2.12|) - H2.14|) . the scalar fluctuations of the geometry are parametrized 
by the 4 scalar functions, i. e. (/>, ip, B and E. The vector fluctuations are described by the 
two (divergenceless) vectors in three (spatial) dimensions Wi and Qj, i.e. by 4 independent 
degrees of freedom. Finally the tensor modes are described by hij, leading, overall, to 2 
independent components because of the last two conditions of Eq. 1)2. 15(1 . 

Slightly different parametrizations exist in the literature. For instance, instead of de- 
noting 6sgij{k,T) = 20^ {T)[ip{k,T)6ij + kikjE{k,T)] the Fourier component of the spatial 
(scalar) fluctuation, the authors of Ref. [H] prefers to separate explicitely the longitudinal 
Hi^kjT) and transverse, HT{k,T) perturbations: 

1 

3^ 



SsQijik^r) = 2o?{t 



HL{k, T)6ij + Hrik, r) ( kikj - -6ij 



(2.16) 



where ki = ki/k. 



In the applications related to the synchronous gauge description the fluctuations are 
customarily denoted, in Fourier space, by |44| I45| I46j 



kikjh{k,T) + 6^{k,T)[ kikj 5i 



(2.17) 



^sgij{k,T) = a^ir) 

Under infinitesimal coordinate transformations, 

r ^ f = r + eo, x' ^ x' = x' + e' (2.18) 
the fluctuations of the geometry defined in Eqs. H2.12() - ()2.14() transform as 

Sgtiu SgtJ.u = Sgfiu - ^t^eu - ^u<^ti, (2-19) 

where is the covariant derivative with respect to the background goemetry and = 
a^(r)(eo, — Ci). The functions eo and are often called gauge parameters since the infinites- 
imal coordinate transformations of the type ()2.18|1 form a group which is in fact the gauge 
group of gravitation. The gauge-fixing procedure, amounts, in four dimensions, to fix the 
four independent functions eo and €{. The gauge parameters can be separated into their 
divergenceless and divergencefull parts, i.e. 

ei = die + Ci, (2.20) 

where diC = 0. The gauge transformations involving eo and e preserve the scalar nature 
of the fluctuations while the gauge transformations parametrized by d preserve the vector 
nature of the fluctuation. 

Since the covariant derivatives appearing in Eq. ()2.19|) (i.e. V = d^ey — T^^^eo-) 
are computed in terms of the unperturbed connections (see Eqs. ()A.7|) of the appendix), 
from Eqs. (|2.12j) ~ (|2.14j) . the fluctuations in the tilded coordinate system, defined by the 
transformation of Eq. ()2.18|) . can be written as 

^ <^ = - T^eo - eo, (2.21) 

^ = ^ + neo, (2.22) 

B^B = B + eo-e', (2.23) 

E ^ E = E - e, (2.24) 
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in the case of the scalar modes of the geometry. Under a coordinate transformation pre- 



serving the vector nature of the fluctuation, i.e. ^ = x* + (with = 0), the 
rotational modes of the geometry transform as 

Q,^Q, = Q^- (2.25) 

Wi^W, = W, + Ci. (2.26) 



The tensor fluctuations, in the parametrization of Eq. 1)2. are automatically invariant 
under infinitesimal diffeomorphisms. 

In the appendix useful technical results concerning the fluctuations of the geometry and 
of the sources are reported: to derive the evolution equations of the fluctuations in the 
different gauges employed in the present paper it is mandatory to have the perturbed form 
of the Christoffel connections and of the Ricci tensors to first order in the metric fluctuations 
of Eqs. (|2.12|1 - (|2.14|) . The fluctuations of the mentioned quantities have to be performed 
in general terms without a specific gauge choice. 

The perturbed components of the energy-momentum tensor can be written, for a single 
species, as: 

ST^ = 5px, 6Ti = -5p^l <5To^ = (pA + PA)5y^\ (2.27) 

where we defined 5np^ = diV^^\ Notice that the perturbed velocity field can be also written 
in a different way (see appendix) . The convention adopted here differs from the one of . 
Under the infinitesimal coordinate transformations of Eq. 1)2. 18|) the fluctuations given in 
Eq. (|2.27j) transform as 

5px 6px = Spx - p'xeo, (2.28) 
6px Spx = 5px - wxp'xeo, (2.29) 

yW ^ yW = yW + g'. (2.30) 

Using the covariant conservation equation ()2.5jl for the background fluid density, the gauge 
transformation for the density contrast, i.e. 6(^x) = Sp{\)/ P(\)-, follows easily from Eq. ()2.28|) : 

5(A) ='J(A) -3W(l + w;(A))eo. (2.31) 
There are now, schematically, three possible strategies 

• a specific gauge can be selected by fixing (completely or partially) the coordinate 
system; 

• gauge-invariant fluctuations of the sources and of the geometry can be separately 
defined; 

• gauge-invariant fluctuations mixing the perturbations of the sources and of the geom- 
etry can be employed. 
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2.3 Gauge choices and gauge-invariant variables 



The choice of the coordinate system and of the relevant gauge-invariant quantities is deter- 
mined by the peculiar features of the problem to be solved. The tensor modes of the geome- 
try, as introduced in the present section, are invariant under infinitesimal diffeomorphisms. 
The vector and the scalar modes do change for infinitesimal coordinate transformations. 
The gauge fixing of the vector fluctuations is closely analog (but technically simpler) to the 
one of the scalar modes and that is why vectors will be discussed before scalars. 

2.3.1 Gauge fixing and gauge-invariant fluctuations for vector modes 

Two different gauge choices appear naturally from the analysis of the vector modes. The 
first choice could be to set Qi = 0. In this case, from Eq. (|2.25|) . the gauge function Qi is 
determined to be 



Since, in this case, the gauge function is determined up to an arbitrary (space-dependent) 
constant, the coordinate system is not completely fixed. This occurrence is reminiscent 
of what happens in the synchronous coordinate system of scalar fiuctuations. The gauge 
choice = has been used, for instance, by the authors of Refs. [51115]. 

Another equally useful choice is the one for which Wi = 0. Starting from an arbitrary 
coordinate system, according to Eq. (|2.26|) . the gauge function is determined as Ci = 
and the gauge freedom, in this case, is completely fixed. The gauge Wi = has been 
exploited in gHl EOl EU • 

Instead of fixing a gauge, one may choose to work directly with gauge-invariant quan- 
tities. By inspecting Eqs. H2.25() and ()2.26|) it is easy to argue that the quantity + Qi 
is invariant under infinitesimal coordinate transformations preserving the vector nature of 
the fluctuation. This is the flrst example of a gauge-invariant potential and it is the vector 
counterpart of the Bardeen potentials deflned originally by Bardeen in 02] (see also [37j). 

2.3.2 Gauge fixing and gauge-invariant fiuctuations for scalar modes 

If the gauge freedom is used in order to eliminate all the off-diagonal entries of the perturbed 
metric (i.e. E and B according to the parametrization of Eqs. 1)2. 12() - (|2.14|) '). we are 
led to the conformally Newtonian (or simply longitudinal) gauge (see |^ E])- In the 
longitudinal gauge £'l = i?L = 0. To reach the conformally Newtonian gauge from an 
arbitrary coordinate system characterized by arbitrary values of E and B one has to require, 
according to Eqs. (|2.23() and ()2.24|) , that the gauge functions are determined as e = E 
and Eq = E' — B. In this gauge the only non- vanishing entries of the perturbed metric are, 
according to the parametrization of Eqs. (|2.12|) ~ (|2.14j) (p and ^/J. 

In a related perspective, a set of gauge-invariant generalization of the longitudinal de- 
grees of freedom can be constructed. This set of gauge-invariant fiuctuations can be written 




(2.32) 
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as 031112] 



$ = + (5 - E'Y + n{B - E') 

^ = ^ -n{B - E'), 



(2.33) 
(2.34) 



for the metric perturbations and 



6pg = 5p + p'{B 
Spg = Sp + p'{B 
= v' + d'E', 



E') 
E') 



(2.35) 
(2.36) 
(2.37) 



for the fluid inhomogeneities, written, for simphcity in the case of a single species. Using the 
transformations derived in Eqs. (|2.21|) - H2.24|) and H2.28() - ()2.3U() . it can be verified directly 
that the quantities defined in Eqs. (|2.33j) ~ (|2.37)) are gauge-invariant under infinitesimal 
coordinate transformations preserving the scalar nature of the fluctuation. Equation 1)2. 37() 
can also be written asVg = v + E' since, for scalar fluctuations, = d'^Vg and = d'^v. It 
is sometimes parctical to discuss not the peculiar velocity but its divergence, i.e. 9 = diV^ 
whose associates gauge- invariant quantity is = diV^. 

The off- diagonal gauge 1^31 (also correctly called uniform curvature gauge ESI 
inSl EZI) is characterized by the conditions E'od = and Vod = 0. Starting from and 
arbitrary coordinate system with non vanishing E and ip, the off-diagonal coordinate system 
is selected by fixing^ eo = —ijj/Ti and e = E. The two relevant metric fluctuations are then, 
in the off-diagonal gauge, B and (p. Since -00(1 = 0, the fluctuation of the spatial curvature 
is vanishing. 

The synchronous coordinate system stipulates that (ps = Bs = 0. Unlike the longitu- 
dinal and the off-diagonal gauges, the synchronous description does not fix completely the 
coordinate system. This aspect can be appreciated from Eqs. (|2.21|) and (|2.23|) . implying 
that, in the synchronous description, the gauge functions are fixed by the following two 
differential relations 



leaving, after integration over r, two undetermined space-dependent functions. This remain- 
ing gauge freedom is closely linked with the presence, in the synchronous description, of 
unphysical gauge modes whose features have been, however, completely understood i44i 115] . 
The parametrization usually employed in the literature is not the one based on i{j and E 
but rather the one already mentioned in Eq. 1)2. 17(1 . 

An important gauge choice is the so-called comoving orthogonal gauge. In the comoving 
gauge the quantity vq + Bq is set to zero. In this gauge the expression of the curvature 
fluctuations coincides with a relevant gauge-invariant expression whose evolution obeys a 
conservation law which will be discussed extensively in section 5. 

Finally, in the uniform density gauge the total matter density is unperturbed. The 
uniform density gauge is then defined from the relation (5/9d = 0. The latter implies that, 
starting from an arbitrary coordinate system where the total density fluctuations are non 

^To demonstrate this point, Eqs. 12.2211 and 12.2411 have to be used in an arbitrary coordinate system. 





(2.38) 
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vanishing, the gauge function has to be fixed as eo = {6p)/p' (see Eq. ()2.28|) '). The curvature 
fiuctuations on constant density hypersurfaces also obey a simple conservation law which will 
be discussed in section 5. If the total energy- momentum tensor of the sources is represented 
by a single scalar field one can also define the uniform field gauge, i.e. the gauge in which 
the scalar field fiuctuation vanishes (see section 5). 

3 Sachs- Wolfe effect 

One of the first concepts relevant for the theory of the CMB anisotropics is the so-called 
Sachs- Wolfe (SW) effect [13]^ • To introduce the SW effect, we recah that in the late- 
time Universe there are three rather close, but physically separate, time scales: the time of 
matter-radiation equality (denoted by Tgq); the time at which photons decouple from matter 
(denoted by Tdec); the time at which electrons and protons recombine to form Hydrogen 
(denoted by Tree)- Recalling the definition of redshift, i.e. z = oo/a — 1, the redshift of 
equality is Zeq ~ 2.4 x lO^r^m^^ while the redshift of decoupling is Zdoc ~ 1100. This 
occurrence implies that, at least, r^ec ~ 4.6req so that the decoupling time takes place 
during the matter dominated epoch. 

The geodesies of the photons that last-scattered at the decoupling time are sensitive to 
the degree of homogeneity of the geometry after equality. If, for sufficiently large length- 
scales after equality (but before decoupling) the geometry of the Universe was slightly 
inhomogeneous the energy of a photon will experience a shift. 

The shift in the energy of a photon as defined in the reference frame of a fiuid moving at 
some velocity relative to the background space-time should then be computed: this is the 
physical rationale for what is usually called the SW effect. The tensor modes, being only 
coupled to the curvature, will behave differently from the scalar modes which are sensitive 
to the evolution of the (scalar) matter sources. These differences will be reflected in the 
SW effect that has to be separately treated for scalar, vector and tensor fiuctuations. 

3.1 Scalar Sachs- Wolfe effect 

It is practical to recall that the geodesies of a photon are invariant under conformal rescal- 
ings. This observation implies that by defining the appropriately rescaled affine parameter 
and the appropriately rescaled metric tensor, some parts of the derivation of the SW ef- 
fect can be conducted as if the photons would propagate in a (slightly inhomogeneous) flat 
space-time. 

Consider then the line element of a spatially flat FRW metric and notice that it can be 
written as 

gaf3dx°'dx^ = a'^{T)gaf3dx°'dx^ , (3.1) 

where r is the conformal time coordinate; ga/s is the flat space-time (Minkowski) metric 
supplemented by the appropriate inhomogeneous part, i.e. in formulae, gap = rjap + ^(jap- 

®For a pedagogical derivation of the SW effect in the case of the standard adiabatic mode (see section 
4), Ref. 1581 can be usefully consulted. Here, however, the relativistic derivation will be presented. 
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Introducing the afiine parameter s, the geodesies of a photon, x^{s) in a spatially flat 
FRW metric (k = in Eq. (|2.1j) ') are the same as those in the perturbed Minkowski metric 
gap with null geodesies x"(r) and with affine parameters related by ds = a'^{T)dT. The 
latter property can be checked from the action of the geodesic of a massless particle 

d^a.Ms)]'^^. (3.2) 

whose specific form implies that under the rescalings 

9a(3 ^ a^{T)gaf3, ds ^ a^{T)dT, (3.3) 

S ^ S = S, i.e. the action is invariant. As a simple application of this observation consider 
the geodesic equation in the inhomogeneous Minkowski metric. As discussed in section 2, it 
is always possible to perform the derivation of a given physical quantity in a specific gauge. 
Being the final results physical, they will also be gauge- invariant. 

The gauge we ought to select is the longitudinal gauge which stipulates that, from Eqs. 
STT^ and f^TM , E = and B = 0, i.e. e = ^ and eo = {E' - B). Thus the perturbed 
Minkowski metric will be characterized by two independent scalar functions, i.e. 

4ffoo = 2(/', SsQij = 2ip6ij. (3.4) 

The perturbation of the photon geodesic can be written as — > x^{t) + SgX^, where 
x^{t) = n^T is the unperturbed photon geodesic; = (l,ra*) is the four-vector giving the 
direction of the photon in a coordinate system where the observer, located at the end of a 
photon world-line, is at x = 0. The perturbed photon geodesic, obtained by perturbing to 
first order in the amplitude of the scalar fluctuations Eq. (|3.2p . can be written as 

dHx^^ dx-dx^ 

where (^sT^^ are the perturbed connections of the inhomogeneous Minkowski metric ga(j. 
Since the connections vanish when evaluated on the background geometry (Minkowski 
space), in the second term at the left hand side of Eq. (|3.5() possible terms linear in 
5x^^ are absent. The explicit form of ^sT^^^ can be directly inferred from Eqs. ()A.in|) of 
the appendix by setting TC = (Minkowski space) and E = B = (longitudinal gauge 
condition). For the ends of the present derivation, the relevant component of the perturbed 
geodesic equation will then read 

'dSr^l _ _ _ 

-JsTgo™"™" - '^;r°.n*n^' - 25sT%n'n'^ = ^' _ ^' _ 2di(l)n'n^, (3.6) 



d_ 



dT 



where the second equality follows from ^gf qq = ^s^^j = —i^'^ij, ^s^oi = di(j). Notice, fur- 
thermore that the relation = n^r has been used. Equation (|3.6() can be easily integrated 
once^ with respect to r between the time Tj (coinciding with the decoupling time) and the 
time Tf (coinciding with the present time) 

^= r\^' + 4,')d^-2c^. (3.7) 



^Notice that integration by parts is necessary in order to integrate the term 2di(pn'. Recall, in fact that 
d(l>/dT = (j)' + di(j>n\ 
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The quantity to be computed, as previously anticipated, is the photon energy as mea- 
sured in the frame of reference of the fluid. Defining as the four- velocity of the fluid and 
P'^ as the photon four-momentum the photon energy to be computed is noting but 

£ = gf^.uf'P''. (3.8) 



The four-momentum of a photon can then be written by generalizing the known special 
relativistic relation 

pfi _ pO ' 



.dx^" P^dx^" El ,. d 



dX a? dr a? 



+ —5x^' 
dr 



(3.9) 



where E is a, parameter (not to be confused with one of the off-diagonal entries of the 
perturbed metric) defining the red-shifting photon energy. The four-velocity of the fluid 
satisfies, within the conventions adopted in this review, the condition 

g^,^u^'u^ = 1. (3.10) 

The logic will now be to determine n^, g^y and to first order in the fluctuations of the 
geometry. This will allow to compute the right hand side of Eq. (|3.8|) as a function of the 
inhomogeneities of the metric. The first-order variation of Eq. 1)3. lUj) leads to 

'5s5oo'"° = -2Js^t°fl'oo, (3.11) 

so that in the longitudinal coordinate system Eq. (|3.11j) gives, to first-order in the metric 
fluctuations, 

?xO = 1(1-0). (3.12) 
a 

Following the conventions of section 2 and of the appendix, the (divergencefull) peculiar 
velocity field is given by 

5y = - = -d'v. (3.13) 
a a 

The relevant peculiar velocity field will be, in this derivation, the baryonic peculiar velocity 
since this is the component emitting and observing (i.e. absorbing) the radiation. In the 
following this identification will be undesrtood and, hence, = v\^. 

The energy of the photon in the frame of reference of the fluid becomes, then 

£ = g^.u^P" = goou^'P'' + giju'P^, (3.14) 

Inserting now Eqs. 1)3. 9(1 and (|3.12|) into Eq. 1)3.14(1 and recalling the explicit forms of goo 
and gij to flrst order in the metric fluctuations we have 



a 



1 + - nj^t, + 



dr 



(3.15) 



Assuming, as previously stated, that the observer, located at the end of a photon geodesic, 
is at X = 0, Eq. ()3.15|) can be expressed as 

f = ^ 1 1 - - riivi + j^^ {il)' + (/-Odrj (3.16) 
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The temperature fluctuation can be expressed by taking the difference between the final 
and initial energies, i.e. 

5T _ af£{rif) - ajEj 
T QiEi 

In Eq. 1)3.17(1 Ei is the initial photon energy which can be expressed as 



(3.17) 



Ei — En 



1 + 



5T 



En 



1 + 



'^r(??i 



(3.18) 



where 5^ = bp-dpr and oc is the radiation density contrast. The final expression for 
the SW effect induced by scalar fluctuations can be written as 



(3.19) 



Sometimes, for simplified esitmates, the temperature fluctuation can then be written, in 
explicit terms, as 



in 



+ 



{iP' + (l)')dT. 



(3.20) 



Equation ((3. 20(1 has three contribution 



the ordinary SW effect given by the first two terms at the righ hand side of Eq. I|3.20() 
i.e. (5i./4 and cj); 

the Doppler term (third term in Eq, (|3.2U|) '): 

the integrated SW effect (last term in Eq, ((3.2U|) '). 



The ordinary SW effect is due both to the intrinsic temperature inhomogeneities on the 
last scattering surface and to the inhomogeneities of the metric. On large angular scales 
the ordinary SW contribution dominates. The Doppler term arises thanks to the relative 
velocity of the emitter and of the receiver. At large angular scales its contribution is 
subleading but it becomes important at smaller scales, i.e. in multipole space, for i ~ 200 
corresponding to the first peak in Fig. ^ The SW integral contributes to the temperature 
anisotropy if ip and (p depend on time. It will be shown that this may be the case for some 
classes of initial conditions (see section 4) but it is not the case for the standard adiabatic 
mode. In quintessence models off dark energy the longitudinal fluctuations may acquire a 
time-dependence in spite of the adiabaticity inducing a late time integrated SW effect. In 
section 8 Boltzmann equation will be used in order to derive, in a more refined way, the 
same formula given in Eq. ()3.20(1 . 



3.2 Vector Sachs- Wolfe effect 

In the case of the vector Sachs- Wolfe effect the shift in the photon geodesies is induced by 
the vector modes of the geometry defined in Eqs. (|2.13() and (|2.14() and named Qi and Wi 
(see also Eqs. ()A.17|) - ()A.19|) of the appendix). Since Wi and Qi are both divergenceless the 
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photon energy in the frame of reference of the baryonic fluid will be, to first-order in the 
amplitude of the vector fluctuations: 



£ 



E 
a 



-riiVi + 



d6^ 
dT 



(3.21) 



where is the rotational component of the baryonic peculiar velocity and, as anticipated 
in section 2, 5jV^ = 0; the notation 5^, as defined in the appendix, denotes the fluctuation 
induced by the vector modes of the geometry. 



From the geodesic equation we have 

d6^x^ 1 



(3.22) 



dr 2 

where, as before, the perturbed connections are obtained by setting Ti = in Eqs. HA.18|) : 
the quantity Vi is defined as 

Vi = Qi + Wl (3.23) 

Recalling the transformation rules for Qi and Wi as derived in Eqs. 1)2. 25(1 and (|2.26|) . it can 
be appreciated that Vi is invariant under infinitesimal diffeomeorphisms. Thus, the induced 
temperature fluctuations induced by the vector modes of the geometry can be written as 



(3.24) 



In the vector case there is no ordinary SW effect but the only terms are the ones connected 
with the Doppler effect and with the integrated SW effect. 

3.3 Tensor Sachs- Wolfe effect 



The tensor Sachs- Wolfe contribution, as anticipated, corresponds only to the integrated 
Sachs- Wolfe effect. In full analogy with the derivation of the scalar Sachs- Wolfe effect, the 
tensor contribution can be also obtained. Recalling that the tensor fluctuations of the metric 
are parametrized by a rank-two tensor in three-dimensions, hij which is divergenceless and 
traceless, from Eqs. (|2.13p and (|A.21|) the perturbed metric and the perturbed connections 
can be obtained. As before, the perturbed Minkowski metric will be considered. The photon 
energy is then 

^ E 
£ = - 1 + 



d6tx^ 



The geodesic equation to be integrated is then 



d_ 



d5fX^^ 



dT 



dT 



(3.25) 



(3.26) 



The explicit value of StT^j can be obtained from Eqs. (|A.21|) by setting Tl = and, 
consequently, 

/ A 1 PTr 

(3.27) 





--f 




2 Jr, 



h'^^n^'n^dT. 

The only contribution to the tensor Sachs- Wolfe effect is given by the Sachs- Wolfe integral. 
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4 Initial conditions for CMB anisotropics 



In the derivation of the SW the possible existence of metric inhomogeneities after equahty 
has been simply assumed. However, the metric inhomogeneities have to be solutions of the 
Einstein equations and this restricts their possible form. 

The effect of metric fluctuations on the temperature anisotropics can be discussed within 
two complementary approaches. The first approach could be called, in broad terms, model 
independent. In this case, the spatial and temporal variations of the metric inhomogeneities 

are determined as functions of a number of unknown initial conditions parametrized in terms 
of a suitable number of arbitrary constants. The second approach could be called model 
dependent. In this case the initial conditions of the fluctuations are given at an early stage 
in the life of the Universe (for instance during conventional inflation or in the context of 
other cosmological models inspired by superstring theory). The fluctuations will then be 
evolved in time using the techniques to be described in section 5 and the result of this 
procedure will then allow to compute the SW effect induced by the different models. 

The two approaches are complementary since the metric inhomogeneities computed in 
a specific model of evolution will always fall in one of the general solutions described in 
the model-independent approach. In other words, once a specific model of early evolution 
is adopted, the arbitrary constants appearing in the model-independent approach will be 
fixed to a specific value. Sometimes, in the literature, the model-independent study of the 
metric inhomogeneities relevant for the SW effect, is swiftly phrased as the problem of 
initial conditions. This terminology will be also adopted here with the caveat that there 
could be ambiguities. In fact, there is another problem of initial conditions which will be 
treated later on (in section 6) : namely the problem of initial conditions of the fluctuations 
in a specific model of the physics of the early Universe. 

To set initial conditions for the CMB anisotropies means to find large-scale (i.e. kr < 
1) solutions for the evolution of the metric inhomogeneities and of the relevant plasma 
quantities. When a given mode satisfies the relation A:r < 1, it is said to be outside the 
horizon. If this is the case, the associated physical frequency is smaller than the Hubble 
rate H, i.e. u = k/a{T) < H. In the opposite case, i.e. kr > 1 the mode is said to be 
inside the horizon and the corresponding physical frequency will be larger than the Hubble 
rate. For A;r ~ 1 the mode k is said to be crossing the horizon. The initial conditions of 
CMB anisotropies can be set, for numerical purposes, deep within the radiation dominated 
epoch. However, if swift estimates of the SW effect have to be obtained, one must recall 
that, since decoupling occurs already in the matter-dominated phase, the solutions of the 
evolution equations have to be modified accordingly. 

A relevant distinction playing a key role in the theory of the CMB anisotropies is the 
one between adiabatic and isocurvaturc ^ initial conditions. Consider, for simplicity, the 
idealized case of a plasma where the only relevant fluid variables are the ones associated 
with CDM particles and radiation. Incidentally, this is a rather good approximation for 

*In order to avoid misunderstandings it would be more appropriate to use the terminology non-adiabatic 
since the term isocurvature may be interpreted as denoting a fluctuation giving rise to a uniform curvature. 
In the following the common terminology will be however used. 
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S^-S,c^ -5c ^ -<5b- (4.2) 



r < Teq. The entropy per dark matter particle will then be given by ? = T^/uc where ric 
is the number density of CDM particles and pc = TncUc is the associated energy density. 
Recalling that 6^ = Sp^/pr and 6c = Spc/pc are, respectively the density contrast in radiation 
and in CDM, the fluctuations of the specific entropy will then be 

5 = ^ = 3^ - 5e = - 5e, (4.1) 

where the second equality follows recalling that pr oc T'^. If the fluctuations in the specific 
entropy vanish, at large-scales, then a chacteristic relation between the density contrasts 
of the various plasma quantities appears, i.e. for a baryon-photon-lepton fluid with CMD 
particles, 

-Sc — -( 
3 3 

Eq. (|4.ip can be generalized to the case of a mixture of different fluids with arbitrary 
equation of state as specifically discussed by Kodama and Sasaki For instance, in the 
case of two fiuids a and b with arbitrary barotropic indices r«a and W]j the fluctuations in 
the specific entropy are 

1 + -Wa l + Wb 

where 5a and (5b are the density contrasts of the two species. It is relevant to stress that, 
according to Eq. 1)2. 31() . giving the gauge variation of the density contrast of a given species, 
5a b is gauge-invariant. 

If the solution of the coupled system of the pertured Einstein and fiuid equations satisfies 
Eq. 1)4. 2|1 . then the initial conditions are said to be adiabatic. This possibility is the one 
favoured by conventional infiationary models. If the condition (|4.2j) is not satisfied by the 
initial conditions are said to be isocurvature. On top of the adiabatic mode, there are, 
in principle, four well characterized non-adiabatic modes: a baryon isocurvature mode, a 
CDM isocurvature mode and two neutrino isocurvature modes. While some of these modes 
are divergent in the longitudinal gauge description, they are all finite and well behaved in 
the synchronous gauge which, therefore, will be used to introduce them. In the following 
section after discussing the problem of the initial conditions in the longitudinal gauge, we 
will move to the synchronous gauge and complete the analysis. 



4.1 Longitudinal gauge description 

From the results derived in Eqs. (|A.14|1 - ()A.16|) . the perturbed components of the Einstein 
tensor in the longitudinal gauge are obtained by setting E = and B = 0. The perturbed 
Einstein equations can then be written as: 

ds^S = SttG^sTO, (4.4) 
5,gi = 87rG6,T^ , (4.5) 
SsG'o = 87TG6,n. (4.6) 
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The fluctuations of total the energy-momentum tensor are written as the sum of the fluc- 
tuations over the various species composing the plasma, i.e. according to Eq. (|2.27|1 . 



6,T^ = 6p = Y,pxdx, (4.7) 

A 

6sT^ = -6l6p + n| = -6{ wxPxSx + n|, (4.8) 

A 

6,T^ = {p + py = 5^(1 + wx)pxvl (4.9) 

A 

In Eqs. (|4.7() - (|4.9() the sum over A run over the different species of the plasma. We will 
be interested in the situation where the universe is deep within the radiation-dominated 
epoch for T < 1 MeV when neutrinos have already decoupled and form a quasi-perfect 
(collisionless) fluid. In this case the species are baryons, photons, neutrinos and CDM 
particles. Notice that if the energy-momentum tensor may include a dark-energy component 
parametrized, for instance, in terms of a cosmological term. In this case the fluctuations 
of the dark-energy are negligible. This situation may change, after decoupling, if the dark 
energy is parametrized in terms of scalar field(s). In this case contributions to the integrated 
SW effect are naturally expected. 

In Eq. (|4.8)) an anisotropic stress 11^ has been included to account for the possible 
presence of collisionless particles, like, in the present case, neutrinos. Equations (|4.4() - ()4.6() 
lead then to the following system of equations: 

- 3W(W0 + = A7rGa'^6p, (4.10) 
V2(W(/> + ip') = -AnGa'^ip + p)e, (4.11) 



< + W(2V' + 0') + (2W' + W^)^ + iv2(0 _ ^) 



-^did^i(P - V) = ATTGa\5p5i - U{). (4.12) 
where the divergence of the total velocity field has been defined as: 

(P + P)0 = Y.^Px + Pa)^a (4.13) 

A 

with 9 = div'^ and 9x = div\. Equations (|4.1U|) and 1)4.11(1 are, respectively, the Hamiltonian 
and the momentum constraint. The enforcement of these two constraints is crucial for the 
regularity of the initial conditions. Taking the trace of Eq. (|4.12|1 and recalling that the 
anisotropic stress is, by definition, traceless (i.e. 11* = 0) it is simple to obtain 

+ 2Hip' + H>p' + {2n' + 7^2)0 + -V2(</. -i))= AirGa^dp. (4.14) 

3 

The difference of Eqs. (|4.12|1 and (|4.14|) leads to 

lv\^ - iP)6{ - Idid^icl) - V) = -A7Ta^GU{. (4.15) 

^ 
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Going then to Fourier space and defining (see appendix) 



^ - i<5i)n| = {p + p)a = J2{px + px)ax, (4.16) 

A 

Then, Eq. (|4.15|) becomes : 

V2((/>-^) = 127rGa2(p + p)cT = 127rGa2^(pA+PA)<TA. (4.17) 

Since, as elaborated before, the only collisionless species for T < 1 MeV are neutrinos, the 
relevant contribution to Eq. (|4.17l) will mainly come from the neutrinos. 

Equations H4.1U() - ()4.12() may be supplemented with the perturbation of the covariant 
conservation of the energy-momentum tensor, i.e. Eqs. ()A.55|) and ()A.56|) : by selecting the 
longitudinal coordinate system (as explained in the appendix), Eqs. ()A.55|1 - HA.56|) become, 
for a single species, 

Spx 



6'^ = {i + wx){3i^' -ex) + m 



WX e 

opx 



6x, (4.18) 



e'x = {3wx - i)nex - ^^0x - -^^^^'h + vVa - v% (4.i9) 

Wx + l Wx + 1 opx 

In the following, the evolution equations of the fluctuations in the longitudinal gauge will 
be exploited for the derivation of the large-scale adiabatic mode. 

4.1.1 The adiabatic mode 

Defining as 5^, S^, 5^ and 6c the neutrino, photon, baryon and CDM density contrasts, the 
Hamiltonian constraint of Eq. (|4.1U|) can be written as 

- 3W(W(/> + V') - A;V = ^n^iiRuSu + (1 - Ru)S^) + ^^b'^b + ^cSc], (4.20) 



where, for N^, species of massless neutrinos, 

11 j ' ^"'l + R' 



^=o^A — ] ' ^'^ = -^^^ R^ = l-Ry, (4.21) 



so that Ry and i?^ represent the fractional contributions of photons and neutrinos to the 
total density at early times deep within the radiation-dominated epoch. Eq. 1)2. 3|) has 
been used into Eq. (|4.2n|) in order to eliminate the explicit dependence upon the total 
energy density of the background. Notice that, in the following the index k appearing in 
the Fourier transformed quantities will be omitted since it could be confused with the other 
indices defining the various species present in the plasma. 

From the momentum constraint of Eq. (|4.11|) . and from Eq. (|4.14() the following pair 
of equations can be derived: 



4 

— {RuOu + R^O^y) + ^b^b + ^C^^C 

o 



(4.22) 

ijj" + (2^' + (t)')'H + {2H' + n^)^ - —((A - ^) = —{Ry6y + 6^R^), (4.23) 
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where following Eq. ()4.13|) the divergence of the (total) peculiar velocity field has been 
separated for the different species. Furthermore, in Eqs (|4.22|) - (|4.23|) . Eq. (|2.4j) has been 
used in order to eliminate the explicit dependence upon the (total) energy and pressure 
densities. Finally, according to Eq. (|4.17|) the neutrino anisotropic stress fixes the difference 
between the two longitudinal fluctuations of the geometry, i.e. 

fc2(</,-^) = -6?t:V^, (4.24) 

completing the set of perturbed Einstein equations written in explicit form. 

Equations (|4.18|) and ()4.19|) should now be written in the case of the various species. In 
the case of CDM particles (i.e. Wc = 0), Eqs. (|4.18j) and (|4.19|) imply 

e'^ + nec = k^cj), (4.25) 

6'^ = 3ip' - 9c. (4.26) 

Photons an baryons are tightly coupled through Thompson scattering and, therefore, are 
treated as collisional species. The fluid equations for photons 

5'^ = -^9., + Ai;' , (4.27) 



2 



9' — __ 

^~ 4 



aneCrT{9h — 9^). (4.28) 



and baryons 



5; = 3V''-eb, (4.29) 



i = -n9b + 4bA:2jb + feV + --^aneXeaT{9^ - 6'b), (4.30) 



are again obtained from Eqs. (|4.18)) and (|4.19|) by recalling that 



i, w^ = 0, C3% = |^. (4.31) 



In Eqs. Km and 

CTT = — = 0.665246 barn, ro = (4.32) 
3 



is the Thompson scattering cross section |59| . Xe is the ionization fraction of the plasma 
and 

Ue = 1.13 X 10"^ X (n^h^) (1 + zf cm-3 (4.33) 

is the electron density. From Eqs. 1)4.32(1 and H4.33() the mean free path of CMB photons 
can be computed prior to recombination as a function of the redshift z: 

£y ^ — - ~ 5 X 10^(rJb/i^)"^(l + z)~2 Mpc. (4.34) 

Xe(TTne{z + 1) ^ 

For comoving scales shorter than the photon mean free path photons and baryons are 
tightly coupled and the characteristic time for the synchronization of the photon and baryon 
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velocities is small with respect to the expansion time and to the oscillation period. In the 
first approximation, ut ^ oo and O-f ~ ^b- This result holds to zeroth order in the 
tight-coupling expansion. Higher orders are also interesting in the context of the CMB 
polarization and will be treated after having introduced the full Boltzmann hierarchy for 
the brightness functions (see section 8). 

The last species to be introduced is given by neutrinos. Unlike baryons and photons 
neutrinos are not collisional and therefore the fluid approximation, for this species is not 
strictly justified. The usual approach is to characterize neutrinos with a set of fluid equations 
where, however, the anisotropic stress contributes directly to the evolution equation of the 
peculiar velocity. Since, in this approach, the neutrino anisotropic stress is dynamical, to 
be fully consistent the evolution of has to be included. The full set of relevant equations 
will then be 

5'y = -Uy + Ai^', (4.35) 
91 = —by - k^Gy + k^(\>, (4.36) 

< = ^Oy - ^kj^y^. (4.37) 

Equations ()4.35|) and 1)4. 36(1 are directly obtained from Eqs. (|4.18|) and 1)4. 19(1 in the case 
Wy = 1/3 and 0",^ 7^ 0. Eq. (|4.37j) is not obtainable in the fluid approximation and the full 
Boltzmann hierarchy has to be introduced. The quantity J^y^ introduced in Eq. 1)4. 37() . is 
the octupole term of the neutrino phase space distribution (see section 8). 

The system of the above equations will now be solved for the case of the adiabatic mode 
deep within the radiation dominated epoch. Inspection of the Hamiltonian and momentum 
constraints (|4.2fl|) - (|4.22j) shows that if the adiabaticity relation of Eq. (|4.2)) is satisfied, 
then the fluctuations, in the conformally Newtonian gauge, should be constant in the limit 
kr < 1. Inspection of the momentum constraint also implies that the peculiar velocity fields 
should be of order k'^r and that they should all be equal to lowest order. Thus, the solution 
of the system can be found as a perturbative expansion in kr and the starting trial solution 
can be written as: 



6y = -2(j)o + Ak^r"^ = S^ = = pc 


(4.38) 




(4.39) 


V' = V'o + C^k'^r'^, 


(4.40) 


0^ = + D^k\\ 


(4.41) 


0^ = ^k'r + Dyk'r^ 


(4.42) 


9^ = ^k\ + D^k\^ 


(4.43) 


e, = hk^r + D,k^r\ 


(4.44) 




(4.45) 
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Notice that the equahty between 6-^ and 5u up to ©(/c^r^) has not been postulated but it 
is a direct consequence of Eqs. (j4.35|) and (|4.27jl . Furthermore, the leading term of the 
velocity fields has been determined by consistency of the trial solution with Eqs. ()4.26|) . 
(ICTl (P3n|l and (jOHl) . 

Consider now the evolution equations for the neutrinos as written in Eqs. (|4.,'-{5|l " (|4.l-i7j) . 
and assume, for the moment, that the octupole moment of the neutrino phase-space density 
vanishes, i.e. J^^z = to the wanted order in kr. In this case, Eqs. (|4.42j) and (|4.45j) imply, 
to lowest order, 

-0 = % (4.46) 

Since, from Eq. (|4.24|) . the value of the anisotropic stress is fixed in terms of the difference 
between the two longitudinal fluctuations of the metric to 

Eq. (|4.46j) implies that, to lowest order, the relation between -00 and 4>q is: 

V^o = (l + ^^.)'Ao. (4.48) 

The remaining arbitrary constants appearing in Eqs. H4.38() - (|4.45|) are fixed by requiring 
the compatibility with the Einstein- fluid system of equations to the relevant order kr. From 
the Hamiltonian and momentum constraints we get, respectively,: 

+ 2C^ = - ^— + ^ 
C^ + 2C^ = 2R^D^ + 2{l - R^)D^. (4.49) 
Equation ()4.23|) implies, after the use of Eq. ()4.48|1 . 

QC^ + C^ + ^Rucpo - 4 = 0- (4-50) 
From Eqs. (|4.3(ij) and (|4.28|) the following algebraic relations can respectively be obtained: 

3D^-^-C^ + ao = 0, 

3D^-^-C^ = 0. (4.51) 
From Eqs. ()4.35|) and 1)4. 27(1 the relation between the constants is 

2A = 8C^ - ^(/.Q. (4.52) 



Finally, from Eqs. (|07)) and we do get 



Ci/i ~ C(f> — '^RiyDiy, 

a, = ^. (4.53) 
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The equations for the baryon and CDM velocity fields, i.e. Eqs. (|4.29|) and 1)4. 25(1 . simply 
imply that 4:D^ c = C'</< both for baryons and CDM particles. By solving the algebraic 
system of equations provided by all the conditions obtained so far in Eqs. (|4.49I) - H4.53() the 
various constants are determined to be ^ 



A 



(525 + 188 + 16 Rl) 



C,/, = 



45 (25 + 2i?^) 

(75 + 14 i?^ - 8 Rl) 
' ~ 90 (25 + 2i?j.) 
{75 + 79R^ + 8Rl) 



90 (25 + 2i?^ 
(25 + 8i?^) 
20 (25 + 2i?^)'^ 
(65 - 
36 {25 + 2R^)' 



^ (65 + 16i?^) , ,^ 

= 15cJi = - ' , ^ '. (Po, (4.54) 



which leads to the solution 



<5b = 


^ 3^ (525 + 188i?^ + 16^2) ^ 
2^" 60 (25 + 2i?^) 


(4.55) 


6y = 


^ (525 + 188i?, + 16i?2) 
45 (25 + 2i?.) 


(4.56) 




^" 90 (25 + 2i?^) 


(4.57) 




^ , 2„ ^^ (75 + 79i?, + 8i?2) 
l^+5^^j^° 90 (25 + 2i?.) 


(4.58) 




<Ao^^2_ (65 + 16i?,) ^4_3 
2 36 (25 + 2 i?^)^" 


(4.59) 


0h = 


'Ao.2_ (75 + 14 i?. - 8 ,_3 
2 360 (25 + 2 i?^) 


(4.60) 


Oc = 


0o,2_ (75 + 14 i?. - 8 ii^) ^ ^^^_3 
2 360 (25 + 2 i?^) 


(4.61) 




<^o^.2_ (25 + 8i?.) 
2 20 (25 + 2i?j.)^^ 


(4.62) 




Cj)0 2^2 (65 + 16i?,) 

15^^ 540 (25 + 2i?,)'^°^^ ' 


(4.63) 



This is the expression of the adiabatic mode. Clearly, the adiabatic mode is determined, 
by fixing the numerical value and spectral dependence of 4>o{k). Having fixed the value of 
SgoQ, the value of ip follows from the value of the anisotropic stress which is proportional to 
the fractional contribution of the neutrinos, i.e. Eq. (|4.48|) . It should also be noticed that, 
to lowest order, all the velocity field are equal. However, as kr approaches 1 the values of 
9^, 6^, 9h and 9c will start becoming different. 

^Notice that not all the conditions obtained from Eq. 14.491 to Eq. 14.531 are independent. 
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4.2 Synchronous gauge description 



There are examples in the hterature of isocurvature modes that are divergent, at early times, 
in the longitudinal gauge, but which are perfectly physical and regular in the synchronous 
coordinate system. Moreover, it is relevant to consider the formulation of the problem 
of initial conditions in the synchronous gauge, since various numerical codes solving the 
Boltzmann hierarchy use, indeed, the synchronous description. 

The synchronous gauge condition stipulates that = and B = 0. The perturbed form 
of the FRW line element is then parametrized in terms of ip and E. For the problem at 
hand, the parametrization given in Eq. (|2.17j) is customarily chosen. The relation between 
the two parametrizations is simply, in Fourier space 

2k'^E = /i + 6^, ij = -C- (4.64) 

Inserting Eq. (|4.64() into the perturbed Einstein tensors of Eqs. (|A.14|) - ()A.16|) and setting 
(j) = B = the fluctuation of the left hand side of Einstein equations can be obtained in the 
synchronous gauge parametrized as in Eq. (|2.17|) . 

By performing an infinitesimal gauge transformation the 6goo and 6goi parts of the met- 
ric (which vanish in the synchronous gauge) also transform. The gauge modes can be exactly 
identified by requiring that Sgoo = and 5gij = in the transformed coordinate system, 
i.e. for gauge transformations that preserve the synchronous nature of the coordinate sys- 
tem. The first gauge mode corresponds to a spatial reparametrization of the constant-time 
hypersurfaces. As a consequence, in this mode the metric perturbation is constant and the 
matter density unperturbed. The second gauge mode corresponds to a spatially dependent 
shift in the time direction. 

With the notation of Eqs. (|4.64|) . the perturbed Einstein equations become 

k^C - ^h' = hi\R^5^ + (1 - Ry)5^ + 17b5b + l^c^c], (4.65) 



k\' = -\n 



URu9u + (1 - Ru)9^) + ^^b^b + ^c0c 



(4.66) 

h" + 2nh' - 2k^i = [Ry5y + {1- Ru)5^] , (4.67) 
h" + 6^" + 2nh' + 12W^' - 2k'^^ = Un'^R^a^, (4.68) 

where the background equations (|2.3|) - (|2.5|) have already been used to eliminate the energy 
and pressure densities^ ^. 

By combining appropriately Eqs. 1)4. 65(1 and 1)4. 67(1 . it is possible to obtain a further 
useful equation 

h" + nh' = [2R^6u + 2(1 - i?^)(5^ + Q^6i, + Qc^c] ■ (4.69) 



^"Notice that sometimes in the literature the variable called ^ is sometimes called -q. This notation will be 
avoided since rj will be reserved for one of the slow-roll parameters in the context of conventional inflationary 
modes. 

^^In this section the density contrasts and the peculiar velocity field are named in the same way as in the 
longitudinal gauge. It is understood that the density contrasts and the peculiar velocity fields are not equal 
in the two gauges and are related by the transformations listed below (see Eqs. 14.801 and (14.821 '). 
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The evolution equations of the pecuhar velocities and density contrasts of the various species 
of the plasma can be obtained by perturbing the covariant conservation. The result of this 
calculation is reported in Eqs. HA.57|) and HA.58|) without any gauge fixing. Hence, setting 
cp = B = 0m Eqs. (TT^TIl and (Tr^S]) and using Eq. I|4.H4|) the following set of equations 
can be readily obtained: 



5l = -^9, + h', 


4 2 
6', = --9j + -h', 


(4.70) 


S'b = -^b + \h', 




(4.71) 


el = -k^a^ + —5y, 


^2 


(4.72) 


d'h = — W6'b, 




(4.73) 




15 5^ 


(4.74) 



As anticipated the synchronous gauge modes can be made harmless by eliminating the 
constant solution for h and by fixing, for instance, the CDM velocity field to zero. These 
two requirements specify the coordinate system completely. 

The adiabatic solution can be also obtained in the synchronous gauge, and in praticular. 



5 + _^ 



k^T\ (4.75) 



.6(15 + 4i?^) 

h = -Ck^T^. (4.76) 

The constant introduced in Eqs. (|4.75|) and (|4.76|) has been defined in order to match the 
standard notation usually employed in the literature (see for instance jlHl) to characterize 
the adiabatic (inflationary) mode. The solution for the density contrasts will be 



2 



5y = 5, = -^Ck\^, (4.77) 



6i,=Sc- ^kW (4.78) 



Setting, as explained, 9c = 0, the peculiar velocity satisfy instead: 

ft --^ ft __'^1±^ 

18' 15 + 4i?/ 



(4.79) 



Any solution obtained in the synchronous coordinate system can be studied in the confor- 
mally Newtonian gauge by using the appropriate transformation, i. e. 

ci^i. = -^m+h)" +-H{Qi+h)'], 

9i^^=9i^^-^{h' + 6a, (4.80) 



31 



where the subscripts refer to the quantities evaluated either in the longitudinal or in the 
synchronous gauge. Using Eq. (|4.8Ujl . from Eos. (|4.75jl - (|4.76jl and (|4.77jl - (|4.79|1 . we recover 
the standard adiabatic mode with the constant C determined as 

C = ^*„. (4.81, 

where 4>o is the value of the adiabatic mode in the longitudinal gauge discussed previously. 
Conversely, the solutions of the longitudinal gauge can also be directly transformed into the 
synchronous gauge using the appropriate transformation, which can be easily derived: 

n 



^ = -V'L / acl)LdT, 

a J 

= 6Vl + 6- / a(j)Ldr, -2k' I / a{T')cl)LiT')dT' , 



4.2.1 Isocurvature modes 

It is useful to recall that, in looking for perturbative solutions (both in the longitudinal and 
in the synchronous gauge), there are various small parameters. One is certainly fcr, which 
is small outside the horizon. However, also fib and f^c are small parameters deep within 
the radiation-dominated epoch. Let us make this statement more precise by considering 
the scale factor 



a(r) 



(4.83) 



interpolating between the radiation-dominated phase for r <C ri and the matter dominated 
epoch for r ^ ti. In this case we can also write 

^b,c = Hb.c ^ , (4.84) 

a(r) + 1 

where the subscripts in Eq. ()4.84|1 refer either to baryons or to CDM and where we took, 
for simplicity, ri = 1. 

Eqs. (|4.83j) and (|4.84jl can be inserted into the evolution equations for the perturbations 
in the synchronous gauge, regarding l^b^c as small paramters deep within the radiation- 
dominated epoch (r 0). The result of this procedure is summarized by the following 
solutions: 

/i~ (-4nbr + 6nbr2), ^ ^Ijl^r -Tl^t'^ , (4.85) 



for the metric perturbations and 

8 
3 



5^ ~ ( --17br + 417b ), (5b =^ (1 -217bT + 317br^), (4.86) 
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5u - (^- 3^br + 40br' j , '^c ^ ^U^t + 3n^T^, (4.87) 

6^ ~ -^nbfcV2, ~ -^HbfcV^, (4.88) 

0e = O, a. ~ -- fcV^. (4.89) 

for the fluid quantities. Clearly, for this mode, the adiabaticity condition is not satisfied. 
Furthermore, by transforming the solution to the Newtonian gauge, it is easy to check that 
the longitudinal fluctuations of the metric vanish for r ^ 0. This is the baryon isocurvature 
mode already discussed in [221 ■ ^ similar solution can be obtained by changing Ob ~^ ^ci 
i.e. the CDM isocurvature mode. The baryon isocurvature mode has been studied in detail 
inRefs. IHl EH ESI iH] • 

On top of the baryon and CDM isocurvature modes, there are other two modes which are 
physical and which may parametrize, together with the modes discussed so far, the most 
general initial conditions for the CMB anisotropies. They are the neutrino isocurvature 
density mode and the neutrino isocurvature velocity mode. These modes were analyzed in 
j64j and partially introduced in |65| . Following the same procedure described in the case 
of the baryon and CDM density modes it can be shown that, for the neutrino isocurvature 
density mode, the metric inhomogeneities are given by 

^"-lOi?;^"' ^" 6(4i^. + 15) ^"' 
while, setting 9c = 0, the fluid quantities are determined to be 

Ru ,2^3 X. r^^^" U2 2 



e,^9,^-\^^k\ + hl,^k'r^ (4.93) 

(Tv - —, r, ^ • (4.94) 

2(15 + 4i?^)' 4 ^ ' 

From the solution of Eq. (|4.94|) we can see that R^6i, + R^5^ ~ 0. So, the initial conditions 
are such that the total energy density is essentially unperturbed. However, as the mode 
enters the horizon, the neutrinos free-stream while photons are tightly coupled with baryons. 

In the case of the neutrino isocurvature velocity mode, the solution for the metric 
inhomogeneities can be written as 

h ~ -Ifih^kr', (4.85) 

^ ' + _ .„^°^: .„ \kr'. (4.96) 



3{5 + AR,y) V Rj (5 + 4i?^)(15 + 4i?j. 
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The fluid variables will instead be given by: 



-^n^^kr^ 5b ^ ^kr - m^?^^^^^^kr\ (4.97) 

4 inx^ 



5, ^ --kr - n^^kr' 5, ^ ^^kr - n^^^^i^J^^kr' , (4.98) 



d^^e^^ -^k + ^nh^kT + 3|f 17b(3i?^ - 3n^)kT^ + ^k'T\ (4.99) 

J9 + 4R^ (4.100) 
6(5 + 4i?^) ' ^ ' 

" ^(i^^T^'" + (5 + 4i?!)(^5 + 4i?.)'^'' ^^-^^^^ 

^--7(5T4i?::)^'^'- (^-^^^^ 

As in the previous examples of non-adiabatic modes, the CDM peculiar velocities are set 
to zero (i.e. 6c = 0) to fix completely the coordinate system. 

In the longitudinal coordinate system the neutrino isocurvature velocity mode leads 
to divergent perturbations in the limit kr — > 0. This point can be easily demonstrated 
either by solving directly the equations in the longitudinal gauge or by transforming the 
metric fluctuations of the neutrino isocurvature velocity mode from the synchronous to the 
longitudinal gauge. Inserting then the two expression of Eq. 1)4. 96() at the right hand side 
of the first two transformation rules given in Eq. (|4.80j) it can be easily obtained that 

which diverge in the limit kr 0. The fact that metric fluctuations diverge in a specific 
description should not be regarded as a physical problem but rather as a mathematical 
difficulty. We shall get back to this point later on in section 6. 

The possibility of having a mixture of modes (adiabatic and isocurvature) as initial con- 
dition for the CMB anisotropics led to a number of mode-independent analysis of various 
sets of data. In [HEl it was argued, on the basis of Maxima data, that mixture of isocur- 
vature and adiabatic modes could not be excluded and in jBT] Enqvist and Kurki-Suonio 
pointed out that the Planck could detect isocurvature modes thanks to its foreseen sensi- 
tivity to polarization. CMB polarization is basically induced by scattering processes (either 
last Thompson scattering or late reionization scattering). This occurrence allows to elimi- 
nate possible contamination with line-of-sight (integrated) effects and makes polarization a 
valuable tool in order to constrain isocurvature modes. In [HB] CDM isocurvature modes in 
open and closed FRW backgrounds have also been considered. 

In |()9|l7()j the importance of polarization was also stressed in view of the Planck planned 
sensitivities. A general (model independent) analysis of the recent WMAP data (combined 
with other sets of data including 2dF galaxy redhift survey) was recently performed in 
[711 172j . The results can be roughly summarized by saying that the isocurvature fraction 
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allowed by the present data ranges from about 10 % (when only one isocurvature mode is 
present on top of the adiabatic mode) up to 40 % (when two isocurvature modes are allowed 
together with their cross-correlations). Finally the fraction of isocurvature modes rises to 
about 60 % when three isocurvature modes are allowed. Notice, as a remark, that in j22] the 
authors indeed allowed for general correlations between adiabatic and isocurvature modes 
while the power spectrum was parametrized by means of a power-law. From the analyses 
reported so far it seems that the constraints on isocurvature modes are rather dependent 
upon the possible correlations with the adiabatic modes. In |73j . following the analysis of 
Ref. inni) the possibility of (correlated) CDM and adiabatic modes has been discussed in 
light of WMAP and large-scale structure data. Different values of the spectral indices for 
the power of each mode (and of their correlation) have been scrutinized. As a consequence, 
the bounds reported in [BS] became more stringent by, roughly, a factor of 1/6. 

One can also address the question on how some correlated mixture of adiabatic and 
isocurvature modes may arise in the early Universe. This subject will be also partially 
addressed in section 6 where some general ideas on how isocurvature modes may be excited 
in the early Universe will be introduced. As an example, consider Ref. where a model 
dependent analysis has been performed. In the context of minimal curvaton models (see 
section 5) the authors found that correlated isocurvature modes seems to be disfavoured 
with respect to pure adiabatic modes. However, the presence of a (totally correlated) baryon 
isocurvature mode does not seem to be ruled out (see also |75| 176]). 

Various sorts of different initial conditions may also be phenomenologically constrained. 
In the following we are going to list some possible interesting departure from the generalized 
(model-independent) analysis of the initial conditions. In the present discussion the dark 
energy has been always parametrized in terms of a cosmological constant. However, if the 
quintessence fluid is described in terms of one (or more) scalar degrees of freedom, then 
the picture may change. Recalling the form of the integrated SW effect, the presence of 
a quintessence field may induce a time evolution in the longitudinal degrees of freedom of 
the metric inducing an important integrated contribution (see for instance |77|). In [JH] 
(see also |79p it was argued that if the primordial fluctuation of the quintessence has a 
correlation with the adiabatic density fluctuations, the large-scale temperature fluctuations 
may be suppressed. This may have implications for the so-called WMAP anomaly for the 
low multipoles of the CMB angular power spectrum. 

Other possible modifications on the initial conditions include the presence of decaying 
modes j80] as well as the presence of other physical modifications such as a primordial 
magnetic field [HB]- Notice that in the latter case, not only vector and tensor modes are 
excited [811 182j but also scalar modes. In this case all the solutions described up to now 
(i.e. the adiabatic mode and the four isocurvature modes) are modified in a computable 
way (53] . 

4.3 Simplified views on the problem of initial conditions 

For numerical purposes it is safer to set initial conditions in the radiation dominated phase 
and integrate across the radiation-matter transition by using, for instance, an interpolating 
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solution of the background Einstein equations (|2.3() - H2.5() (like the one discussed in Eq. 
(|4.83jl . However, since the decoupling time follows the equality time and last scattering 
occurs when the Universe is dominated by matter, it is also possible to study directly the 
problem of initial conditions during the matter epoch. After introducing some general ana- 
lytic considerations on the problem of initial conditions during the matter-dominated epoch 
the adiabatic and isocurvature contributions to the ordinary SW effect will be computed. 

The fluctuation in the total pressure density can be connected to the total fluctuation 
of the energy density as 

5p = cl6p + 5pr,^A, (4.104) 

where 



2 



(4.105) 



is the speed of sound computed from the variation of the total pressure and energy density 
at constant specific entropy, i.e. 6<^ = 0. The second term appearing in Eq. (|4.1U4|) is the 
pressure density variation produced by the fluctuation in the specific entropy at constant 
energy density, i.e. 

' 5p 



5pn.d = [^j^j <5?, (4.106) 

accounting for the non-adiabatic contribution to the total pressure perturbation. 

If only one species is present with equation of state p = wp, then it follows from the def- 
inition that = w and the non-adiabatic contribution vanishes. As previously anticipated 
around Eq. (|4.1() . a sufficient condition in order to have 5pnad 7^ is that the fiuctuation in 
the specific entropy 5<; is not vanishing. Consider, for simplicity, the case of a plasma made 
of radiation and CDM particles. In this case the speed of sound and the non-adiabatic 
contribution can be easily computed and they are: 

4 = l = i±4^i(^J^^), (4.107) 



P' P'.y + Pc 3\3pc+4:p, 



"7 

f5p\ 4 5p^ _4/ pc/>7 \_ 2 ^/i ^ns^ 

Pj Pc 

To obtain the final expression appearing at the right-hand-side of Eq. I|4.1()7|) the conserva- 
tion equations for the two species (i.e. p'^ = —ATLp^ and p'^ = —STipc) have been used. In 
Eq. (|4.1()8|) . the first equality follows from the fiuctuation of the specific entropy computed 
in Eq. 1)4. 1|) ; the second equality appearing in Eq. H4.108|) follows from the observation that 
the increment of the pressure should be computed for constant (total) energy density, i.e. 
6p = 6p^ + 6pc = 0, implying 6pc = —Sp-y. The third equality appearing in Eq. H4.108j) is a 
mere consequence of the explicit expression of Cg obtained in Eq. (|4.107j) . 

As in the case of Eq. (|4.3|) . the analysis presented up to now can be easily generalized 
to a mixture of fiuids "a" and "b" with barotropic indices tUa and t(;b |llj . The generalized 
speed of sound is then given by 

^2 ^ Wa(Wa + l)Pa + Wb(Wb + l)ph ^^g. 
{W^ + l)/3a + {Wh + l)ph 
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From Eq. ()4.1U7() and from the definition of the (total) barotropic index it fohows that, 
for a CDM-radiation fluid 



Using Eqs. H4.1U6|) - H4.108() the non-adiabatic contribution to the total pressure fluctuation 
becomes 

4 S 

Wd = oPcTT— 7' (4-112) 
6 6a + 4 

where the definition given in Eq. (|4.1|1 . i.e. S = ((5?)/?, has been used. 

Using the splitting of the total pressure density fluctuation into a adiabatic and a non- 
adiabatic parts, Eq. ()4.1U|1 can be multiplied by a factor Cg and subtracted from Eq. 
(|4.14j) . The result of this operation leads to a formally simple expression for the evolution 
of curvature fluctuations in the longitudinal gauge, namely: 

i;" +n[<j)' + i2+3c^,)^'] + [n\l + 2c'^,) + 2n'](l)-c^,V^^+\v\^-i^) = 47TGaHp,,^d, (4.113) 

which is indipendent of the speciflc form of 5pnad- The left hand side of Eq. (|4.113() can be 
written as the (conformal) time derivative of a single scalar function whose speciflc form is, 

Taking now the flrst (conformal) time derivative of TZ as expressed by Eq. (|4. 114(1 and using 
the deflnition of Cg we arrive at the following expression 

^' = ' ATrGa^ip + p) ^"^" + ^^^^ + + '^'^ + t^^' + ^^"^^ ^ l)^']'^}- (4-115) 

Comparing now Eqs. (|4.115|1 and (|4.113() . it is clear that Eq. ()4. 115(1 reproduces Eq. ((4.113() 
but only up to the spatial gradients. Hence, using Eq. (|4.115() into Eq. (|4. 113(1 the following 
flnal expression can be obtained: 



/<- = : OPnad + 



p + p ATTGa?{p + p) 



(4.116) 



There are slightly different ways to express the terms at the right hand side of Eq. ()4. 116(1 . 
For instance the term ((/> — ij:) could be replaced by the anisotropic stress a. However, 
the physical content of the equation remains unchanged. The speciflc geometrical meaning 
of IZ will be discussed in detail in section 5 and Eq. (|4. 116(1 should be regarded, for the 
moment, as a useful mathematical simpliflcation allowing, in some cases, a swifter reduction 
to quadratures of the evolution of 'ip for kr -^1. 
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Disregarding the complication of an anisotropic stress (i.e., from Eq. 
from Eqs. (|4.18|) - (|4.19|) the covariant conservation equations become 



Si 



Si 



= ^) 

(4.117) 
(4.118) 

(4.119) 



(4.120) 



4.3.1 Adiabatic contribution to the SW effect 

Since in the adiabatic case 5pnad = 0, Eq. 1)4.1131) becomes 

iP" + m^p' + —^p = 0, (4.121) 
o 

where it has been used that Cg = 1/3. 

Since during radiation a(r) ~ r, Eq. H4.121() can be solved as a combination of Bessel 
functions of order 3/2 which can be expressed, in turn, as a combination of trigonometric 
functions weighted by inverse powers of their argument 

. .V cos y — sin y ^ ,-, sin y + cos y , . 

^P{k,T) = Ai{k)^ ^— ^+Bi{k)^ ^— ^, 4.122 

y6 yi 

where y = kr / -v/3. If the solution parametrized by the arbitrary constant Ai[k), ip ^ ipj- for 
fcr <C 1 ( where V'r is a constant). This is the case of purely adiabatic initial conditions. If, 
on the contrary, Ai{k) is set to zero, then ip will not go to a constant. This second solution 
is important in the case of the non-adiabatic modes and will be discussed later. 

In the case of adiabatic fluctuations the constant mode ipr matches to a constant mode 
during the subsequent matter dominated epoch. In fact, during the matter dominated 
epoch and under the same assumptions of absence of anisotropic stresses Eq. H4.113() . in 
the case Cg = 0, leads to the following simple equations 

tp" + snip' = 0, (4.123) 

Since, after equality, a^r) ~ r^, the solution of Eq. H4.123() is then 

iP{k,T)='tP^ + Di{k)(^^y, (4.124) 

where ipm is a constant. Typilcally the term multiplying Di(k) is of the order of 10~^ at 
T = Tdec- The values of tpr and V'm are different but can be easily connected. In fact we 
are interested in wave-numbers kr < 1 after equality. In this limit Eq. H4.116|) implies that 
TZ' = 0. More specifically, if the background is dominated by a perfect fluid with barotropic 
index w, Eq. H4.114|) leads to 
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where Eqs. (|2.3|) - H2.5|) have been used. 

From Eq. ()4. 125(1 we obtain that, during radiation (i.e. w = 1/3), TZ^ = —(3/2)^/^1.) while 
during matter (i.e. w = 0) we have T^m = — (5/3)V'm- Since, by virtue of Eq. ()4.116|1 . TZ is 
conserved in the hmit A;r <C 1, then TZ^i = 7?.r, implying 

V'm = ^A- (4.126) 

The same result can be also obtained by integrating, from Eq. ()4. 114(1 . the evolution of 
tp directly in terms of the scale factor. From Eq. (j4.114() . changing variable from r to a the 
equation to be solved will be 

dip 5a + 6 3(3a + 4) 

-T + TT^TT^ = V'r .' ^ 4.127 
da 2a[a + 1) 4a(a + 1) 

whose integration, after imposing adiabatic initial conditions for a — > ( deep within the 
radiation-dominated phase), leads to 

^(«) = T^{16(\/^TT - 1) + a[a{9a + 2) - 8]}; (4.128) 
the limit for a — > oo (matter-dominated phase) of the right hand side of Eq. ((4.128(1 leads 

to (9/10)V'r- 

Combining Eqs. ((4.119() and ((4.12()() in the presence of the constant adiabatic mode tpj^i 
and during the matter-dominated phase 

6" + k'^c% = -4c2A;Vm (4.129) 
where Cs = l/\/3. The full solution of Eqs. (I4.n8l) - (|4.'I2()|1 will then be 

5c = -2i:^ - ^^T^ (4.130) 



0c = ^V'm, (4.131) 

6r = ^i^m[cosikc,T)-3], (4.132) 

= ^ sin (kcsT). (4.133) 
v3 

As in the case of the adiabatic mode during the radiation-dominated epoch, also in the 
matter epoch, for kr <^ 1, 9c — 9^- 

The ordinary SW effect can now be roughly estimated. Consider Eq. ((3.19(1 in the case 
of the pure adiabatic mode. Since the longitudinal degrees of freedom of the metric are 
roughly constant and equal, inserting the solution of Eqs. ((4.13fl(l - ((4.133(l into Eq. ((3.19(1 
the following result can be obtained 

^ ^(t"^^) = ^COs(A;Csrdec) = ^V'rCOs(A:CsTdec), (4.134) 

where the third equality follows from the relation between the constant modes during radi- 
ation and matter, i.e. Eq. ((4.126(1 . Concerning Eq. ((4.134() few comments are in order: 
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• for superhorizon modes the baryon peculiar velocity does not contribute to the leading 
result of the SW effect; 

• for k Cs Tdec ^ 1 the temperature fluctuations induced by the adiabatic mode are 
simply ; 

• even if more accurate results on the temperature fluctuations on small angular scales 
can be obtained from a systematic expansion in the inverse of the differential optical 
depth (tight coupling expansion, to be discussed in section 8), Eq. (|4.134|) suggests 
that the first true peak in the temperature fluctuations is located at A;Cs?7dcc — 



In this discussion, the role of the baryons has been completely neglected. In section 8 a 
more refined picture of the acoustic oscillations will be developed and it will be shown that 
the inclusion of baryons induces a shift of the first Doppler peak. 



4.3.2 Non-adiabatic contribution to the SW effect 

In the simplistic case of CDM-radiation plasma a rather instructive derivation of the gross 
features of the non-adiabatic mode can also be obtained. If Sp^ad / 0, the evolution of TZ 
is, in the limit kr <^ 1 becomes: 

Eq. H4.135|) can be easily obtained inserting Eq. ()4.112|) into Eq. ()4. 116(1 and recalling 
that, in the physical system under consideration, {p + p) = Pc + (4/3)/Or- In the case of 
the CDM-radiation isocurvature mode, the non-adiabatic contribution is non-vanishing and 
proportional to S. Furthermore, it can be easily shown that the fluctuations of the entropy 
density, S are roughly constant (up to logarithmic corrections) for kr <^ 1, i.e. for the 
modes which are relevant for the SW effect after equality. This conclusion can be easily 
derived by subtracting Eq. (|4.117|) from 3/4 of Eq. (|4.119|) . Recalling the definition of S 
the result is 

S' = -{e,-e^). (4.136) 

Since 0r and 9c vanish in the limit kr <^ 1, S is indeed constant. The explicit form of 9^- 
and 9c can be obtained (in the limit kr ^ 1) from the asymptotic solution (obtained in 
the synchronous gauge, Eqs. ()4.88|) - (|4.89|) ') applying the gauge transformation (|4.8U() to 
re-express the fluctuations in the longitudinal description. 

Eq. (j4.135|) can then be integrated in explicit terms, across the radiation-matter tran- 
sition 

7^ = _45 r — ^ -S^-, (4.137) 
io (36 + 4)2 3a + 4' ^ ^ 

implying that 7^ — > for a ^ (radiation-dominated phase) and that TZ — > —5/3 for 
a — > cxD (matter-dominated phase). 
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Recalling again the explicit form of TZ in terms of il^, i.e. Eq. ()4. 114(1 . Eq. (|4.137|) leads 
to a simple equation giving the evolution of ip for modes /cr <C 1, i.e. 

dip 5a + 6 , S / , N 

-r- + — 7 rV' = r- 4.138) 

da 2a{a + l)^ 2(a + l) ^ ^ 

The solution of Eq. (|4. 138(1 can be simply obtained imposing the isocurvature boundary 
condition, i.e. tpiO) — > 0: 



^(a) = ;^{16(1 - \/^n) + a[8 + a{a - 2)]}. (4.139) 

Eq. H4.139|) is similar to Eq. (|4. 128(1 but with few crucial differences. According to Eq. 
((4.139() (and unlike Eq. ((4.1281) '). ip{a) vanishes, for o — > 0, as Sa/8. In the limit a — > oo 
V'(a) — >■ 5/5. This is the growth of the adiabatic mode triggered, during the transition from 
radiation to matter, by the presence of the non-adiabatic pressure density fluctuation. 

Having obtained the evolution of the evolution of the total density contrasts and of 
the total peculiar velocity field can be immediately obtained by solving the Hamiltonian 
and momentum constraints of Eqs. ((4.1U() and ((4.11() with respect to 6p and 6 

' +5e^ = -2fV' + -^y (4.140) 



p a + 1 a + 1 \ din a 

2/c2(a + l) / dip \ , , 



(3a + 4) V dlna 
which also implies, for fcr <C 1, 

Equation ((4.142() is indeed consistent with the result that the total velocity field is negligible 
for modes outside the horizon. Inserting Eq. ((4.139() into Eqs. ((2.28() - ((4.19() it can be easily 
argued that the total density contrast goes to zero for a ^ 0, while, for a ^ oo we have the 
following relations 

2 1 

6c ^ — cS ~ -2ip ~ — 6r (4.143) 
5 2 

The first two equalities in Eq. ((4.143() follow from the asymptotics of Eq. ((4.14U() . the 
last equality follows from the conservation law (valid for isocurvature modes) which can be 
derived from Eq. ((4.119() . i.e. 

6r ~ 4ip. (4.144) 

Thanks to the above results, the contribution to the scalar Sachs- Wolfe effect can be 
obtained in the case of the CDM-radiation non-adiabatic mode. From Eq. ((3.19() we have 



(^Y^=(^-L + i) ^2iP^.,^^ls, (4.145) 



nad 

where the second equality follows from Eq. ((4.144() and the third equality follows from Eq. 
((4.139() in the limit a — > oo (i.e. a S> acq)- 

The following comments are in order: 
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• as in the case of the adiabatic mode, also in the case of non-adiabatic mode in the 
CDM-radiation system, the pecuUar velocity does not contribute to the SW effect; 

• for k Cs Tdcc ^ 1 the temperature fluctuations induced by the adiabatic mode are 
simply 2V'nad (unlike the adiabatic case) ; 

• Equation (|4.145jl suggests that the first true peak in the temperature fluctuations is 
located at fecsr/dcc — vr/2. 

The last conclusion comes from an analysis similar to the one conducted in the case of the 
adiabatic mode but with the crucial difference that, in the case of the isocurvature mode, 
V' vanishes as r at early times. This occurrence implies the presence of sinusoidal (rather 
than cosinusoidal) oscillations. This point will be further discussed in section 8. 

5 Evolution of metric fluctuations 

In the following the evolution of the tensor, vector and scalar modes of the geometry will be 
summarized. The tensor modes of the geometry are, as repeatedly stressed, invariant under 
gauge transformation. The vector modes are not gauge-invariant, however, their description 
is rather simple if the Universe expands and in the absence of vector sources. The scalar 
modes are the most difficult ones: they are not gauge- invariant and they directly couple 
to the (scalar) sources of the background geometry. The analysis will be conducted in the 
case of a spatially flat background geometry and in the case of different models of early 
Universe. 

5.1 Evolution of the tensor modes 

For the tensor modes the perturbed Einstein equations imply that ^tRj = 0. Hence, 
according to Eq. HA.23|) . the following equation should be satisfied, in Fourier space, by the 
two tensor polarizations: 



where and obey the same evolution equation (|5.1|) and will be denoted, in the 
remaining part of this section, by h. Equation (|5.1|) can be written in two slightly different 
(but mathematically equivalent) forms: 



Recall that St, as discussed in the appendix, denotes the first-order fluctuation with respect to the tensor 
modes of the geometry. Similarly 5s and 5v denote, respectively, the first-order scalar and vector fiuctuations. 





(5.3) 



k /iyfc = 



(5.4) 
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where fik = ah^. Equation H5.3|l can be interpreted as a conservation law for the amphtude 
of the tensor mode. 

In the hmit kr <^ 1 Eq. (|5.3|) imphes that the amphtude of the gravitational wave 
evolves as 

hu{r) ^M + B,j ^^=A, + B,J (5.5) 

where Aj. and Bj. are integration constants to be determined from the analysis of the (early) 
initial conditions which will be discussed in section 6. 

By parametrizing the scale factor of a spatially flat Universe as a{t) ~ t", the Universe 
expands and accelerates (i.e. d > and o > 0) for a > 1 and decelerates (i.e. d > and 
d < 0) for < a < 1. Consequently, Eq. 1)5. 5|) implies 

/ife(t) + (5.6) 

According to Eq. (|5.6)) . if the Universe inflates (i.e. a > 1), then hk will have a constant 
mode (proportional to A^) and a decaying mode (proportional to B^). If the Universe 
decelerates the situation can be different since, on top of the constant mode, there could 
be a growing mode (for a < 1/3). The case a = 1/3 corresponds to the case of a fluid with 
stiff equation of state, namely a fluid where the speed of sound coincides with the speed of 
light In this case the solution is not given by Eq. (|5.(i|) but rather by — Ak + Int 
j84j . This case is particularly relevant for quintessential inflation. 

Eq. ()5.4() has the same physical content of Eq. 1)5. 3|) . For k"^ ^ \a" /a\ the solutions of 
Eq. (|5.3|) are oscillatory, and the amplitude of the tensor modes, /i, decreases as 1/a, if the 
background expands. For fc^ <C \a" /a\ the solution is non-oscillatory and is given by 

r dr' 

fikir) ^ Aka{T) + Bka{T) ^— . (5.7) 

For instance, in the case of de Sitter expansion, the previous equation implies that h is 
constant in this regime. If k"^ <C \a" /a\, the mode is said to be super-adiabatically amplified 
[851 1861 187j . The time-dependent function a" /a is customarily named pump field. This 
terminology is borrowed from quantum optics where a time dependent electric field may 
produce coherent (or squeezed) states of the electromagnetic field |HH1 lEHl IHOl 1^ • The 
evolution of the tensor modes of the geometry can be deduced from the time dependence of 
the pump field. An example is illustrated in Fig. |21where the evolution of \a" /a\ is reported 
for a model where an infiationary phase ("i" stage) turns into radiation at a conformal 
time Tr after an intermediate phase ("s" stage). At Tm the Universe becomes dominated by 
matter. Both during a de Sitter stage and a matter-dominated stage a" /a = 2/t'^. During 
radiation a" /a = since a(r) ~ r. When a given mode with comoving wave-number k 
gets under the potential barrier the corresponding amplitude is amplified during a time 
inversely proportional to the magnitude of k. The amplitude corresponding to /c ~ 
will be minimally amplified. The Fourier amplitudes corresponding io k < will be 
maximally amplified. For k S> Tj^^ the Fourier amplitudes will not be amplified. 

^^By solving the Friedmann equations with barotropic fluid sources (p — wp), a = 2/(3(w + 1)) leading, 
for w = 1, to Q = 1/3. 
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Figure 2: The time evolution of the pump field appearing in Eq. (|5.4|) is illustrated for 
a background evolution where the inflationary stage is followed by a radiation dominated 
stage after an intermediate regime. The second barrier on the right corresponds to the 
radiation-matter transition. 

When a given mode is in the super-adiabatic regime, the corresponding physical fre- 
quency is also, grossly speaking, smaller than the Hubble rate: indeed for an expanding 
FRW background is of the order of H and a" /a is of the order of a^H^. In Fig. |31 
the evolution of the Hubble rate is illustrated as a function of the logarithm (in ten basis) 
of the scale factor. The knee corresponds to the inflation-radiation transition. The small 
ankle at the right corresponds to the radiation-matter transition. Superimposed to the 
plot of the Hubble rate the evolution of two physical frequencies is also reported. The first 
square marks the moment when the given mode becomes smaller than the Hubble rate. 
This moment is called first crossing as opposed to the second crossing (second square at the 
right) when the mode becomes, again larger than the Hubble rate. There exist a number 
of equivalent terminology. For instance, when the given mode is smaller than the Hubble 
rate, its corresponding physical frequency is larger than the Hubble radius (i.e. H^^) or, 
in the jargon, superhorizon sized. 

5.1.1 Hamiltonians for the tensor problem 

The variable //^ = ahj. has a special status in the theory of initial conditions since it is the 
canonical variable diagonalizing the action perturbed to second-order in the amplitude of 
the tensor fiuctuations. The quadratic action for the tensor modes of the geometry can be 
obtained by perturbing the Einstein-Hilbert action to second-order in the amplitude of the 
tensor fiuctuations of the metric. To this end, it is practical to start from a form of the 
gravitational action which excludes, from the beginning, the known total derivative. In the 
case of tensor modes this calculation is rather straightforward and the result can be written 
as 





where r]aj3 is the Minkowski metric. 
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Figure 3: The evolution of the Hubble rate in Planck units is illustrated in a schematic 
model of post-inflationary evolution. The evolution of two different physical frequencies is 
also illustrated in the case when the duration of inflation is non- minimal (about 80 efolds). 
With the dotted line some fundamental physical scale (smaller than the Planck mass) is 
also illustrated. 



After redefinition of the tensor amplitude through the Planck mass 
the action 1)5. 8() becomes, for a single tensor polarization, 



(5.9) 



(5.10) 



whose canonical momentum is simply given by 11 = a^h'. The classical Hamiltonian asso- 
ciated with Eq. (|5.1fl|) is then 



1 



(5.11) 



The Hamiltonian 1)5. 11() is time-dependent: it is always possible to perform time-dependent 
canonical transformations, leading to a different Hamiltonian that will be classically equiv- 
alent to ()5.11() . Defining the rescaled field, // = ah, the action ()5.1U|) becomes 



./2 



while the associated Hamiltonian can be written as 



2n|l^l' + TTn^ - {dill) 



tt' + 2ninT + {^^^ly 

* Recall that in the present conventions £p = Mp ^ — \/8tvG. 



d?x 



(5.12) 



(5.13) 
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in terms of /x and of the canonically conjugate momentum tt = fi' — Tifi. A further canonical 
transformation can be performed starting from (|5.13|) . Defining the generating functional 
of the canonical transformation in terms of the old fields fi and of the new momenta vf, 

J^2^3ifi, n,T)= J d^x (^fi^ - ^/x^^ , (5.14) 

the new Hamiltonian can be obtained by taking the partial (time) derivative of ()5.14|) . with 
the result 

(5.15) 



where, recalling the definition of vr, from (|5.14j) we have vf = /x'. 

The possibility of defining different Hamiltonians is connected with the possibility of 
dropping total (non covariant) derivatives from the action of tensor fluctuations. 

To illustrate this point consider indeed the action given in Eq. 1)5. 1U() for the single 
polarization which can be rewritten as: 

5gw = d^xdTa'^{T)[h'^ - {d^hf]. (5.16) 

By looking at Eq. (|5.16j) it is natural to define the canonical field as = ah. Recalling that 
from fi = ah we also have ah' = (fi' — Tifi), Eq. 1)5. 16() becomes 

5'gw = ^ y" d^xdT[fi'^ + - 2H/X;u' - (5i/x)2]. (5.17) 

Recalling now that the action can be always written in terms of an appropriate Lagrangian 
density as 

5g^ = f dTL{T), L{t) = f d^xC{T,x). (5.18) 



where L{t) is the Lagrangian and C{t,x) is the Lagrangian density. In the case of Eq. 
1)5. 17() . the canonically conjugate momentum is obtained by functional derivation of the 
associated Lagrangian density with respect to /i' and the result will be vr = /x' — Tifj,. Hence, 
the Hamiltonian will be 

^gw(r) = J d^x[TTfi' - C^{t, x)]. (5.19) 

Consequently, from Eqs. (|5.17j) and recalling the expression of tt the Hamiltonian will 
exactly be the one defined in Eq. (|5.13|) . 

Going back to Eq. ()5.17|) and noticing that the term TCfxfj,' can be rewritten as 



W/x2 



(5.20) 



the action 



Sg^ = lJ d'^xy^ + (w' + n')^l^ - {dii^f] (5.21) 

will be classically equivalent to the action of Eq. ()5.17|) since the same Euler-Lagrange 
equations can be obtained from both actions. Now, however, the conjugate momentum will 
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be vf = /i'. Hence, the Hamiltonian can be easily obtained, again, from Eq. H5.19() but with 
TT replaced by n. After performing this simple calculation the result will exactly be the one 
fiven in Eq. ()5.15() . 

The possibility of defining different Hamiltonians is inherent to the time-dependent 
nature of the problem. The same conclusion also applies for the canonical treatment of 
scalar fluctuations as it will be later illustrated in this section. In the Lagrangian approach, 
this aspect reflects in the possibility of defining different classical actions that differ by 
the addition (or the subtraction) of a total time derivative. While at a classical level the 
equivalence among the different actions is complete, at the quantum level this equivalence 
is somehow broken since the minimization of one Hamiltonian or the other may lead to 
computable differences if some recently debated normalization prescriptions are adopted 
(see section 6). 

Up to now the evolution of the tensor modes of the geometry has been discussed in the 
case of expanding FRW backgrounds. Several extensions of this analysis have been studied 
in the last ten years. For instance the present discussion can be extended to the context 
of non-Einsteinian theories of gravity like in ,!)2 93, .94, to the case of cosmological 
models containing compactified extra-dimensions j96| 19 7j and higher derivative extension of 
general relativity [HSI • All these techniques are relevant for the calculation of the stochastic 
background of relic gravitons in unconventional models of the early Universe such as the 
pre-big bang ^92, ^93^ and the ekpyrotic j99i models. 

5.2 Evolution of the vector modes 

The evolution of the vector modes of the geometry can be obtained by perturbing the 
Einstein equations and the covariant conservation equation with respect to the vector modes 
of the geometry, i.e. 

6.0"^ = 87rG,5,T;, (5.22) 
V^d^T^'" = 0. (5.23) 

Recalhng Eqs. (|A.18|) - (|A.19I) of the appendix, the exphcit form Eqs. (f^?^ and (|^^ can 

be obtained 

V^Qi = -16ttG{p + p)a^Vi, (5.24) 
Q', + 2HQi = 0, (5.25) 
V'i + (1 - 3'u;)WH = 0, (5.26) 

where Vj is the divergenceless part of the velocity field, and where as discussed in section 
2, the gauge has been completely fixed by requiring Wi = 0, implying, according to Eqs. 
(|2.25j) and p.26)) . that Q = —Wi. Notice, furthermore, that according to the conventions 
used to represent the velocity field (see Eq. (|A.30|) '). 5^Tq = {p + /o)V*. 

Eqs. 1)5. 24() . (|5.25|) and (|5.26l) follow, respectively, from the (Oi)-(zj) components of 
Einstein equations and from the spatial component of the covariant conservation equation. 
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Eqs. ()5.24|) - H5.25|) are not all indipendent. In particular it is easy to check that using Eq. 
(|5.24j) . Qi can be eliminated from Eq. I|5.25jl and Eq. (|5.2tj|) is easily obtained. 

By Fourier transforming Eqs. ()5.24|1 and (|5.25l) . the following solution can be found: 

0.(r) = §H, (5.27) 

where Ci{k) is an integration constant. 

Two distinct situations are now possible. If the Universe is expanding Qi is always 
decreasing. However, Vj may also increase. Consider, indeed, the case of single barotropic 
fluid p = wp. Then, since according to Eq. H2.3() - ()2.5() . p ~ a"'^'-'"^^-*, 

Vi oc a^""-^. (5.29) 

li w = 1, Vi increases as r (since a(r) ~ V^)! while Qi{T) ~ r^^, and it is decaying for 
large conformal times. This observation was diuscussed, in the past, in connection with the 
idea that the early Universe could be dominated by vorticity jlUUl llUlj . Recently a similar 
observation was also put forward in Ref. jlU2j in the framework of the so-called holographic 
proposal which implies the presence of a barotropic fluid with stiff equation of state during 
the post-inflationary evolution |1U31 11U4] . 

If the Universe contracts the evolution of the vector modes is reversed. In this case Qi 
may even increase j491 15U| . Finally, the situation may change even more radically if the 
theory is not of Einstein-Hilbert type [H^ or if the theory is higher-dimensional [SII- Notice 
that if the energy-momentum tensor of the fluid sources possesses a non-vanishing torque 
force rotational perturbations can be copiously produced |3H|. 



5.3 Evolution of the scalar modes 

In analogy with the case of the tensor modes of the geometry, large-scale scalar fluctuations 
follow a set of conservations law that are only approximate for fcr <C 1 and which are exact 
only in the limit kr — > 0. Here, various important gauge-invariant quantities will be studied 
and interpreted in light of specific gauge-dependent descriptions. 

In section 4 the relevant set of equations describing the evolution of the metric fluctu- 
ations has been derived in the longitudinal gauge. This set of equations is formed by Eqs. 
(|I?TU)) - (|irn|) and Eqs, (|^mi) ~ (|I?T7)) . Covariant conservation of the evolution of the fluid 
sources leads to Eqs. ()4.18() ~ H4.19(1 . Always in section 4, the CDM-radiation isocurvature 
mode has been treated by defining the non-adiabatic pressure density fluctuation, i.e. (5pnad- 
The introduction of this quantity allowed to derive Eq. i|4.11H|) . i.e. the evolution equation 
for TZ that can be regarded, in a broad sense, the scalar analog of Eq. 1)5. 3() that holds 
in the case of tensor modes of the geometry. Eq. (|5.l-ij) . does not have a source term. On 
the contrary, Eq. H4.116|) does have a source term, proportional to (5pnad- The following 
comments are then in order: 
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• the conservation of TZ is enforced in the A; — > hmit provided the non-adiabatic 
variation of the total pressure vanishes; 

• the absence of (5pnad is not guaranteed a priori, unless the content of the plasma is 
only given by a single species (i.e. a single scalar field, a single fluid); in all other 
situations 5pnad will not be vanishing and its effects should be taken into account; 

• various models share the feature of leading to a vanishing (5pnad likcj for instance, the 
case of single field inflationary models. 

The evolution equations for the fluctuations of the geometry and of the sources can be 
first written without fixing a specific coordinate system. Then, by using the definitions 
of the Bardeen potentials, i.e. Eqs. 1)2. 33() - (|2. 37(1 . the wanted gauge-invariant form of the 
equations can be obtained. 

Inserting the definitions (|2.33|) - (|2.37j) into Eqs. (|A.5fl|) and ()A.51|) . the gauge-invariant 
forms of the Hamiltonian and momentum constraints becomes 



where, as remarked after Eq. (|2.37|) . 6 = diV^ is the divergence of the gauge-invariant 
(scalar) peculiar velocity field. 

Inserting the expressions of the gauge-invariant fluctuations given in Eqs. (|2.33|) - (|2.37jl 
into Eqs. HA.53|) and HA.54|) the following set of equations can be also derived: 



Finally, the gauge-invariant form of the perturbed covariant conservation equations can be 
written as: 



as it can be checked by using Eqs. ^rm - ^TT7\i into Eqs. (|X371) - (|X38l) that are written 
without imposing any gauge condition. 

Equations 1)5. 3U() - ()5. 31(1 and ()5.32|) - H5.33|) have the same form of the evolution equations 
in the longitudinal gauge. All the manipulations leading to the evolution of TZ (see, in 
particular, the derivation preceding Eq. (|4.116|) ') can be repeated with the same results. 
The explicit expression of TZ in gauge-invariant terms can then be written as 



(5.30) 
(5.31) 




(5.32) 
(5.33) 



5p'g - 3(p + p)^' + {p + p)e + 3n{6pg + 5pg) = 0, 

{p + p)Q' + [{p' + p') + An{p + p)]Q + VHpg + ip + p)V^<^ = 0. 



(5.34) 
(5.35) 



7^ = -{^ - nVg) = - 



H(^' + w«>) 



(5.36) 
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5.3.1 Physical interpretation of curvature perturbations 



In the comoving gauge the three velocity of the fluid vanishes, i.e. vc = 0. Since the 
constant (conformal) time hypersurfaces should be orthogonal to the four-velocity, we will 
also have that Bq = 0. In this gauge the curvature perturbation can be computed directly 
from the expressions of the perturbed Christoff'el connections bearing in mind that, unlike 
the appendix, we want to compute here the fluctuations in the spatial curvature, namely 

= <5s7'^'^!f + 7'^'<5si?Sf = ^^'^c, (5.37) 

where the subscript C refers to the fact that the calculation has been conducted in the 
on comoving hypersurfaces. The curvature fluctuations of the comoving gauge can be 
connected to the fluctuations in a different gauge characterized by a different value of the 
time coordinate, i.e. rc — > r = rc + eo- Under this shift 

V-c ^ = ■00 + neo, (5.38) 
{vc + Bc)^B + v = {vc + Be) + eo, (5.39) 

Since in the comoving orthogonal gauge Bq + vq = 0, Eqs. H5.38() - (|5.39() imply, in the new 
coordinate system, 

^^ = ^l;-n{v + B) = {^ -nVg), (5.40) 

where the second equality follows from the deflnitions of gauge-invariant fluctuations given 
in Eqs (|2.33|) - H2.37|) . Prom Eq. 1)5. 40(1 it can be concluded that TZ corresponds to the 
curvature fluctuations of the spatial curvature on comoving orthogonal hypersurfaces. 

Other quantities, defined in specific gauges, turn out to have a gauge-invariant inter- 
pretation. Take, for instance, the curvature fluctuations on constant density hypersurfaces, 
V'D- Under inflnitesimal gauge transformations 

V'D ^ = V'D + Weo, 

6pB 5p = 5pi) - peo. (5.41) 

But on constant density hypersurfaces dpD = 0, by deflnition of uniform density gauge. 
Hence, from Eq. 1)5. 41(1 . 

^^ = ^p + n^- = ^ + n^, (5.42) 
p p 

where, again, the second equality follows from the deflnitions of gauge-invariant fluctuations 
given in Eqs. (|2.33|) - H2.37|) . Hence, the (gauge-invariant) curvature fluctuations on constant 
density hypersurfaces, can be deflned as 

C = -(^+W^)- (5-43) 

Equation ()5.43|) deflnes the gauge-invariant curvature fluctuations on constant density hy- 
persurfaces. The last sentence seems to contain a contradiction, but such an expression 
simply means that C coincides with the curvature fluctuations in the uniform density gauge. 
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In a different gauge (for instance the longitudinal gauge) (" has the same value but does not 
coincide with the curvature fluctuations. 



The values of C and TZ are equal up to terms proportional to the laplacian of This 
can be shown by using explicitely the definitions of TZ and ( whose difference gives 

C-n = -n(v^ + ^-^y (5.44) 

Recalling now that from the Hamiltonian and momentum constraints and from the conser- 
vation of the energy density of the background 

= - . ^ + ^'), 



^P?. 1 rv72 



p' A'irGa'^{p + p) 

we obtain 



[V2^-3W(W$ + ^')]> (5-45) 



The density contrast on comoving orthogonal hypersurfaces can be also defined as |431l47j 
= ^ ^ + ^'(^ + ^) = 5p^-m{p + p)V^ 

where the second equality follows from the first one by using the definitions of gauge- 
invariant fluctuations. Again, using the Hamiltonian and momentum constraints, 

^ V^^ (5.48) 



that means, according to Eq. (|5.46j) . that ((" — 7^) oc e 



5.3.2 Curvature perturbations induced by scalar fields 

Scalar fields are an almost universal ingredient for the discussion of potentially physical 
effects in the early Universe since, together with relativistic fluids, they have the virtue 
of leaving unbroken the spatial isotropy of the homogeneous background geometry. Conse- 
quently scalar flelds may appear in diverse (and sometimes unrelated) cosmological contexts 
ranging from conventional inflationary models to minimal |114|lllf^ and non minimal pre- 
big bang models jll6llll71lll8j . Ekpyrotic models jll9lll21Lll22j . together with their cychc 
extension |123j are formulated in terms of scalar degrees of freedom. Finally, quintessential 

^^Various conventional inflationary models have been developed so far. A partial list should probably 
include: old inflation 105 , new inflation 1107) . cahotic inflation |108j . natural inflation |109) . eter- 
nal inflation |110| . hybrid inflation |lll| . quintessential inflation |112| . By the word conventional we mean 
that the above models share the feature of being formulated in terms of a quasi-de Sitter phase of acceler- 
ated expansion driven by a phenomenological master field generically named inflaton. See also |118| for a 
comprehensive review of inflationary models. 
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models of dark energy j42j are connected with the presence of a scalar degree of freedom 
whose potential is carefully chosen in such a way that the scalar field energy density domi- 
nates just today, say for z between and 3. 

The physical meaning attributed to the scalar degree of freedom changes from model 
to model. For instance, in pre-big bang models the scalar fields will be the string theoretic 
dilaton or the Kalb-Ramond axion (which is a pseudo-scalar field in four space-time dimen- 
sions). In the ekpyrotic model the scalar field will parametrize, at least in some versions 
of the model, the effective distance between two colliding domain-walls defined in five di- 
mensions. It is important to remark that the scalar degree of freedom appearing in models 
inspired by string/M theory is the result of the compactification (achieved with diverse 
methods) of a more fundamental degrees of freedom which are typically defined either in 10 
or 11 dimensions. It is difficult, sometimes, to attribute a specific meaning to the infiaton 
field of conventional infiationary models: the infiaton is often a master field (light during 
the infiationary phase) that drives a phase of quasi-de Sitter expansion. 

Some of the considerations developed earlier in this section will now be specialized to 
the case of scalar field sources characterized, rather generically, by a potential W{ip). In a 
spatially flat FRW geometry the background Einstein equations read, in the single scalar 
fleld case, 

W2 = ^f^ + W^a'y (5.49) 



3 V 2 

^" + 2n^' + = 0. (5.51) 

dip 

The fluctuation of i-p changes, for inflnitesimal coordinate transformations, as 

X ^ X = X - '/''eo; (5.52) 

consequently the associated gauge-invariant scalar field fiuctuation is given by 

X = x^v\B-E'). (5.53) 

The two Bardeen potentials and <I>) and the gauge-invariant scalar field fiuctuation X 
define the coupled system of scalar fluctuations of the geometry: 

dW 



- 3n{n<S> + ^') = 4ttG 



-ip''^^ + ip'X' + ^a^X 



(5.54) 



7^$ + ^' = AttG^'X, (5.55) 



.^'2^ + ^'X' - ^a^X 

dip 



(5.56) 



where Eqs. (|5.54|) - (|5.56|) are, respectively, the perturbed (00), (Oi) and {ij) components 
of Einstein equations. Eqs. H5.54() - (|5.56() are sufficient to determine the evolution of the 
system, however, it is also appropriate to recall the gauge-invariant form of the perturbed 
Klein-Gordon equation: 

fl2TIA ATI/ 

X" + 2nX' - V^X + ^-^a^X + 2$— - ip'(^' + 3^') = 0. (5.57) 

dip 
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Eq. H5.57|) can be obtained from the perturbed Klein-Gordon equation without gauge fixing, 
reported in Eq. HA.49|) of the appendix. The gauge-invariant fluctuations of the metric of 
Eqs. (ESSll-dSSil) can be inserted into Eq. (|A.49|) together with Eq. Eq. ((^37)) 

will then appear as a result of the simplifications. 

If the perturbed energy-momentum tensor has vansihing anisotropic stress, the {i ^ j) 
component of the perturbed Einstein equations leads to <I> = ^' 

The gauge-invariant curvature fluctuations on comoving orthogonal hypersurfaces are, 
for a scalar field source, 

X mn^ + ^'') 

where the equality follows from the use of Eq. (|5.55jl and of Eq. (|5.5Uj) . 

The definition of (|5.58|) and a linear combination of Eqs. 1)5. 54p and (|5.56|) leads to the 
following simple equation 

n' = 2^^*' (5-59) 



implying the constancy of TZ for modes /cr <C 1 |43l 1471 11241 1125j . The power spectrum of 
the scalar modes amplified during the inflationary phase is customarily expressed in terms 
of 7^, which is conserved on super-horizon scales (this statement will be made more specific 
at the end of this section while discussing the Weinberg theorem). 

Curvature perturbations on comoving spatial hypersurfaces can also be simply related 
to curvature perturbations on the constant density hypersurfaces, denoted by C, 

C = -^-W^ = -^ + ^. (5.60) 

Taking the difference in the definitions (|5.36l) and (|5.6UI) . and using Eq. ()5.54|) : 

iip 6 if' 

TZ and C differ by Laplacians of the Bardeen potential. 

Taking the time derivative of Eq. ()5.59|) and using, repeatedly, Eq. ()5.58)) and Eqs. 
H5.54p - (|5.56() . the following second-order equation can be obtained: 

TZ" + 2-TZ' - V^TZ = 0, (5.62) 
where the function z{t) is given by 

z{t) = ^. (5.63) 

^®It would not be correct to assume that ip — (f> already at the leven of Eq. l|A.49fl giving the perturbed 
Klein-Gordon equation without gauge fixing: by doing so, the correct form of Eq. 15.5711 cannot be obtained 
since ip and (j) change in a different way under infinitesimal coordinate transformations and hence lead to 
different gauge-invariant fiuctuations. 
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z" 



Eq. H5.62() is formally similar to Eq. (jS-lf) . i.e. the evolution equation of a single tensor 
polarization. This similarity becomes somehow deeper in the specific cases when the scale 
factor has a power-law dependence. In this case, let us suppose that the scale factor has a 
power-law dependence such as a{t) ~ in cosmic time. 

The case of power-law inflation is characterized by a scale factor evolving as a(t) ~ 
where p > 1. Exact solutions of Eqs. 1)5. 49() ~ ()5. 51(1 can be found with exponential potentials 
[12fi| 11271 112S| [T29] . In this case the solution can be parametrized, in cosmic time, as 

a(t) r^tP, ip= z{t) = J-Mpa{t) (5.64) 

t y p 

where, as usual, the dot denotes a derivation with respect to the cosmic time coordinate t. 
Clearly, in this case z{t) and a{t) are proportional and, therefore, Eqs. ()5.62|) and 1)5. 1() are 
exactly equivalent. This equivalence surfaces again if we eliminate the first time derivative 
appearing in Eq. 1)5. 62() by introducing the gauge-invariant variable |13L)l I131j 

q = -zU = aX + z^i, (5.65) 

which leads, once inserted into Eq. (|5.62() . to 

r 

-q - V^g = 0, (5.66) 

z 

that is the scalar analog of Eq. (|5.4|) . In the case of power-law infiation given in Eq. (|5.64|) . 
Eqs. H5.4|) and Eq. 1)5. 66() share the same pump field. This occurrence implies that Fig. |31is 
often interpreted, in some vulgarization of the subject, as holding both for tensor and scalar 
modes. Figure 01 also describes approximately the evolution of curvature perturbations in 
the case of a single scalar field and in the absence of any fluid sources or of any other 
spectator flelds which do not drive infiation but whose presence generate a Sp^ad, i.e. a 
non-adiabatic fiuctuation of the pressure density. The related caveat is that, unlike the 
evolution of the tensor modes, the evolution of the scalar modes is constrained from both 
the (00) and (Oi) components of the perturbed Einstein equations. 

As discussed above, in the case of power-law inflation the simplicity of the (exponential) 
potential leads to a privileged relation between the scalar and tensor pump flelds that turn 
out to be exactly equal, i.e. a" /a = z"/z. For more general inflationary potentials (but 
always in the single fleld case) this relation is not verified anymore. The pump fields of the 
scalar and tensor modes for more general potentials leading to a quasi-de Sitter expansion of 
the geometry can be discussed in terms of the slow-roll parameters e and t] , i.e. jl,S2|ll,S;-S] 



^'^Notice that in |133| r; is defined with opposite sign with respect to our conventions. Moreover some 
authors |132| like to define Stgij = —2a'^hij while we define StQij = —a^hij. So the tensor fluctuation defined 
in |132| is twice the one defined here. This choice is arbitrary provided everything is done consistently (for 
instance, as far as the tensor fluctuation is concerned, the tensor Sachs- Wolfe effect will be sensitive to the 
conventions chosen in order to write the perturbed line element). In spite of the arbitrariness of the definition 
this may lead to confusions as far as the normalization of the power spectra is concerned (see section 6). In 
the following, all the conventions previously established will be carefully followed mentioning, when needed, 
the major differences with the other choices made in the literature. 
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In the case of power-law potentials the slow-roll parameters have the property of being 
constants and equal in modulus, i.e e = 1/p = — r/. 

To find the expressions of the pump fields in terms of the slow-roll parameters it is useful 
to write the background equations as a first-order set of non-linear differential equations 

una insi inni nszi 



2M2 2M4 ' 



dip 

a IT 

^ = -2M2 — , (5.69) 

where, as usual, the adopted definition for the Planck mass will be the one already made 
explicit in the introduction and after Eq. 1)5. 8|) . 

Bearing in mind Eq. ()5.67|) . the explicit relation determining the form of the tensor 
pump field follows from the chain of equalities 

— = T-C"^ n' = 2a^H^ + a^H = 2a'^H^ f 1 - (5-70) 



where in the last equality, Eq. ()5.67|) has been used. In order to find a suitable expression 
we have to determine aH as a function of r. To do so we must recall that the relation 
between cosmic and coformal time depends upon e. In fact, if e is constant the following 
equalities hold thanks to integration by parts: 

dt r da 1 f da 

The third equality in Eq. (|5.71|) follows from integration by parts and by taking e to be 
constant. Equation (|5.7H) allows to determine aH as a function of r, i.e. 



Equation ()5.72|) implies: 

a" /i2 - 1/4 



/^=7T7^ (5-73) 



e)2' ^ 2(1 -e) 



Recall now the definition of z as discussed in Eq. 1)5. 63(1 . The same calculation leading 
to Eq. (|5.73|) implies, for the scalar pump field 

z" 1^^ - 1/4 3 + e + 2r? 

The calculation is rather straightforward: z" /z has to be written in terms of the explicit 
derivatives ofTC, a and ip' as they appear in Eq. (|5.63|) : then the conformal time expression 
have to be translated in theur cosmic time analog; then the definitions of Eq. 1)5. 67|) must 
be used. Ntice that the indices fj, and u of Eqs. (|5.73|) and (|5.74jl have been introduced for 
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later convenience. Moreover, fj, = v m the case of power-law inflation when e = —77. If e is 
not constant Eq. I|5.72|) is generalized to 

' ' '^<' + '')+ (5.75) 



aHl-e aH 

as it can be concluded from integration by parts without assuming the constancy of e; the 
ellipses stand for higher orders in the slow-roll expansion. 

To conclude, it is useful to point out that different decoupled equations for the gauge- 
invariant fluctuations can be obtained, in the single scalar fleld case. The first equation 
of this sort is The other equation to be introduced is a decoupled equation for the 

Bardeen potential. Taking the difference of Eqs. ()5.56() and 1)5.54(1 we obtain rather easily 

^" + QH^' + {2n' + m^)^ - = -SvrG^a^X (5.76) 

dip 

To eliminate X, the momentum constraint of Eq. ()5.55() can be used in the form 

^' + 

Inserting Eq. (|5.77|) into Eq. (|5.76|) and using Eq. 1)5. 51|) to eliminate the derivative of W 
with respect to (p we obtain the following decoupled equation: 



^" + [m - ^' + 2\^n' - ^j"^ - V^^ = 0. (5.78) 

The first derivative of a second order equation with time-dependent coefficients can be 
always eliminated by a suitable field re-definition. It is the easy to show that, defining a 
new variable, be it \I' = (a/p')^, Eq. ()5.78|) can be recast in the form 

^" _ - = 0, (5.79) 

where z has been defined in Eq. ()5.63|) . Let us remark that, apparently, Eqs. ()5.78|) and 
()5.79j) are rather simple. However, the price for simplicity is, in this case, that the derivation 
of Eq. (|5.78|) assumes that ip' never goes to zero, since, by solving the momentum constraint 
(see Eq. ()5.77|) ) we divided by (p' . The effect of this manipulation is evident from the 
expressions of the pump fields of both Eqs. (|5.78|) and ()5.79|1 that contain (explicitly or 
implicitly) inverse powers of 99'. 

Up to now we mainly dealt with the case of a single scalar field. Extensions of some 
of the ideas discussed so far are certainly possible. Consider, for instance, the case of two 
scalar fields ip and a characterized by a potential term W{p,a). In this case the relevant 
Friedmann equations (|5.49() and (|5.5U() can be simply obtained by doubling the kinetic 
terms, i.e. by replacing ip''^ with {p''^ + o"'^). The spatial curvature perturbation can be 
written, in this case (see for instance |1381 113i?l ll4Uj ). as 
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where 



lip 



H 



(5.81) 
(5.82) 



are, respectively, the canonically normahzed fluctuations of if and a, i.e. the generalization 



of the variable q obeying eq. 
can be written as 



. As discussed in jl38j the relevant evolution equations 



+ 3Hq^ 
qa + 3Hqa 



1 



+ 



1 



1 d 



M^a? dt\H 
1 d (a^ 



H 



d fa^ 
H 



Mia? 



dt\ 



1 



d (a^ 



(5.83) 



Equation (|5.83jl . as well as the whole set of formulae derived in the two-field case, can 
also be generalized to the situation of an arbitrary number of scalar fields |141j . This is 
particularly useful from the point of view of the Hamiltonian description of the problem. 
The most notable aspect of the previous two equations is that, unlike the single field case, 
the evolution of q^p and is coupled. This already signals a fundamental physical difference 
with respect to the single field case. 

In fact, in the two-field case the evolution of Tt develops a source term which is the analog 
of the non-adiabatic fluctuation of the pressure density already introduced in section 4 when 
talking about the longitudinal description of isocurvature modes. Direct calculation jl38| 
shows that the evolution of TZ is, in this case. 
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where <5|?nad is now 



5Pnad = (Cs - 1)^(9^ + f^) + (1 - cl){!JpX^ + dX,) - (1 + C^^ 



dW 



where Cg is the generalized sound speed, i.e. 

2 

P 



c? = ^ = l + 



(dW dW 



3H{ip'^ + a2) V da 



(5.86) 



In the first equality of Eq. (|5.86|) p and p are, respectively, the total pressure and energy 
density of the system written in terms of the two background fields <p and a. 



5.3.3 Hamiltonians for the scalar problem 



In light of possible applications, it is desirable to treat the evolution of the scalar fluctuations 
of the geometry in terms of a suitable variational principle. On this basis, Hamiltonians 
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for the evolution of the fluctuations can be defined. By perturbing the action of the scalar 
fluctuations of the geometry, the final form of the action can be expressed in terms of the 
curvature fluctuations 



S. 



(1) 

seal 



/2 



{dillf 



(5.87) 



Defining now the canonical momentum tt-ji = z^VJ the Hamiltonian related to the action 
1)5. 87() becomes 
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seal 
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and the Hamilton equations 

tt' = z^V^n, n' = ^. (5.89) 

Z'^ 

Equations (|5.89|) can be combined in a single second order equation so that Eq. (|5.62|) is 
recovered. 

The canonically conjugate momentum, ■n-ji is related to the density contrast on comoving 
hypersurfaces, namely, in the case of a single scalar field source I124j . 



bp^ + m{p^ + p^)v _ a^bp^ + m^'x 



Pip 



a?p^ 



(5.90) 



where the second equality can be obtained using that p^p+Pip = '^''^ /o? and that the effective 
"velocity" field in the case of a scalar field is 1^ = X/ip' . Making now use of Eq. (|5.55|) into 
Eq. 1)5. 54() . Eq. ()5.9U|) can be expressed as 



aP'P^ 



3 ^2 



(5.91) 



where the last equality follows from Eq. (|2.3I) . From Eq. 1)5. 91() . it also follows that 

VTT^ = z^lZ' = -Qa^He^, (5.92) 

where Eq. (|5.91j) has been used together with the expression of IZ' coming from (|5.59|) . 
Hence, in this description, while the canonical field is the curvature fluctuations on co- 
moving spatial hypersurfaces, the canonical momentum is the density contrast on the same 
hypersurfaces. 

To bring the second-order action in the simple form of Eq. (|5.87jl . various (non-covariant) 
total derivatives have been dropped. Hence, there is always the freedom of redeining the 
canonical fields through time-dependent functions of the background geometry. In particular 
the action (|5.87|) can be rewritten in terms of the variable q defined in Eq. 1)5. 65(1 . Then 
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whose related Hamiltonian and canonical momentum are, respectively 
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(5.93) 



(5.94) 
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In Eq. H5.93() a further total derivative term can be dropped, leading to another action: 



-(3) 
seal 



- / d^X 



r - {^^qf + 



and another Hamiltonian 



(5.95) 



(5.96) 



where vr = q' . As in the case of the Hamiltonians for the tensor modes, Eqs. 
(|5.94j) and (|5.96|) are all related by canonical transformations. Furthermore, notice that, 
classically, the scalar end tensor Hamiltonians have exactly the same form in the case of 
power-law inflation. 



5.3.4 Weinberg's theorem 

Concerning the the evolution of scalar fluctuations, a plausible question to ask is how many 
adiabatic solutions and how many non-adiabatic solutions are compatible, for instance, 
with the conservation of TZ and This question, usually approached within the separate 
Universe picture |47j (see also |142) for a reintroduction of some of arguments given in 
[17]). was recently addressed by Weinberg on a more solid ground in a series of papers 
jl43lll4jlll45j . The separate Universe picture amounts to stipulate that any portion of the 
Universe that is larger than {H being the Hubble rate) but smaller than the physical 
wavelength on the perturbation a/k will look like a separate unperturbed Universe. In the 
following the separate Universe picture will not be invoked but the way of reasoning put 
forward in refs. |143l 11441 1145j will be briefly outlined. The hypotheses of the theorem are 
the following: 

• suppose that the theory of gravity is, for simplicity, of Einstein-Hilbert type; 

• suppose that the background geometry is given by a FRW metric which we can take 
for simplicity to be spatially flat; 

• suppose that the Universe is always expanding (possibly with different rates); 

• suppose that the sources of the geometry are relativistic fluids. 

The last assumption is quite general and it allows for the presence of anisotropic stresses. 
From the mentioned assumptions it also follows that Eq. (j4.116jl is fully valid. This is 
true for two separate reasons. The first reason is that Eq. ()4.116|) is deduced within an 
Einsteinian theory of gravity. The second reason is that, since the Universe is assumed to 
be always expanding, the time-dependent coefficients of the terms at the right hand side 
of Eq. ()4.116() are all non-singular. Moreover, since fluids with different barotropic indices 
can be simultaneously present, the non-adiabatic fluctuation of the pressure density, (5pnad 
can be present. 

The thesis of the theorem is that the evolution equations, in the longitudinal gauge 
description, have always a pair of physical solutions for which Jpnad and TZ approaches 
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a constant for A; ^ 0. In other words, foUowing the usual terminofogy, there will be always 
at least a pair of adiabatic solutions with (5pnad = and TZ constant in the limit k ^ 0: 
one solution with TZ ^ and the other with TZ = 0. In order to understand correctly the 
thesis of the theorem the limit A; — > is crucial. The second point to be borne in mind 
is that there could also be other solutions with TZ' ^ 0. These solutions, if allowed by the 
Hamiltonian constraint, will be non-adiabatic. 

The theorem can be demonstrated in general terms, as illustrated in Ref. jl4;-ij . For 
purposes of presentation it is sufficient to discuss it in the simple case of a scalar field source. 
In this case, the evolution equations in the longitudinal gauge can be simply obtained from 
Eqs. H5.54() - (|5.56() and 1)5. 57() by substituting the gauge-invariant fluctuations with the 
corresponding longitudinal fluctuations, i.e. <I> — > 0, ^ ^ ^ and X ^ x- 

Having done so, consider now the evolution equations in the limit A; ^ (or, equivalently, 
the limit of vanishing spatial gradients). The key observation in order to prove the theorem 
in general terms is that these equations have an extra accidental symmetry for a particular 
class of coordinate transformations, i.e. 

r^f = r + eo(T), ^ = - Xx\ (5.97) 

where the parameter A is a space-time constant. This coordinate transformation clearly in- 
duces a transformation in the fluctuations defined in the longitudinal gauge, in particular, 
using Eq. (|2.19l) and recalling the definitions of longitudinal fluctuations, the transforma- 
tions will be: 

(j) ^ ^ = (P - T-Ceo - e'o, tj; ^ = + neo - X, X^X = X-^'(^o- (5-98) 

While the gauge transformation for (p is exactly the one derived in Eq. (|2.21|) . the trans- 
formation for ip derived on the basis of 1)5. 97(1 differs from Eq. ()2.22() . There is no surprise 
for this occurrence: in fact we are here taking the limit A: ^ 0. Hence, the part of the 
perturbed metric proportional to kikj (which would be generically induced by a coordinate 
transformation) is not induced by the specific coordinate transformation given in Eq. (|5.97|) . 
Therefore, from Eqs. (|5.98|) it is clear that a solution of the longitudinal gauge equations 
can be obtained, in the limit A; — > 0, by setting 

<^ = -Co - Ti-eo, ip = Tieo - A, x = -^W ■ (5-99) 

This solution can be generalized to the case when the sources are more general, namely, 
for instance, several fluids, several interacting scalar fields and so on. The solution given in 
(|5.99j) would simply be a pure gauge. However, it can be extended to the case of non-zero 
k by looking at the conditions implied on eo by the components of the perturbed Einstein 
equations which are trivially satisfied in the limit A; ^ 0: these are the off-diagonal terms 
in the (ij) equation, implying ^ = {^u the absence of an anisotropic stress) and the 
momentum constraint of Eq. (|5.55j) implying ^' -|- 7^-0 = ^.'kGi^'x- The condition (j) = ip 
implies, according to Eq. 1)5. 99(1 that eo satisfies the following simple equation 

eo + 2Weo = A, (5.100) 
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whose solution can be written as eo = Cq^^ + ^q'^ 



^ f\\r')dr', 4'^ = ^, (5.101) 



" a2(T) 7 ^ ' ' " a2(r) 

are, respectively, a particular solution of the inhomoheneous equation and the general so- 

(2) 

lution of the homogeneous equation. The solution expressed by eg can be also understood 
from the freedom of shifting the lower limit of the integral appearing in e^^^ . 

With this observation, the system of the longitudinal fluctuations can be solved. In 
order to appreciate the significance of the results to be obtained, it is useful to write down 
the system as a third order differential system in the variables ij^k, Xk ^i^d = fk- Using 
Eqs. (|5.55j) and (|5.57j) the system to be solved is, in Fourier space 



V'fc = -WVfc + 47rG(^'Xfc, (5.102) 

A = fk. (5.103) 

d'^W dW 

fk = -2nfk - k\k - ^r^a\k - 2^IJk—a' + A^j'^^' . (5.104) 



In this case the coefficients of the various terms appearing in the system are continuous in 
a neighborhood of A; = 0. This condition guarantees that the would-be-gauge mode can be 
lifted to the status of physical solution in the limit k ^ Q. In the general case (when the 
sources are more complicated than a single scalar field) this condition may not be satisfied 
in general but it is certainly rather plausible given the typical forms of energy-momentum 
tensors customarily employed in cosmological model building. 

After repeated use of the background equations, (and in particular of Eq. (|5.51|) ). Eq. 
(|5.in4j) implies for A; = that 

^k = \(^a). (5.105) 

Now, since 



.(1) _ ..(2) _ (2) 



xr = -4^V', X'u' = -4^V', (5.106) 



we will also have that 



^ = -\ + XH a^{T')dT\ (5.107) 
V'f = C2^. (5.108) 

Knowing what are Xk \ Eqs. (|5.107|) and (|5.108|) allow to determine what is TZk in the 
limit — > from the, by now familiar, relation 

7^1^'^) = -^^^)-^X^^). (5.109) 



^*In the following the index denoting Fourier transformed quantities will be restored since no ambiguity 
can arise with different subscripts. 
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In particular, inserting Eqs. ()5.1U6() and H5.1U7|) - ()5.1U8|) into Eq. H5.109|) we will have 



K = ^' K = 0- (5-110) 

The result is: for A; — > there always two adiabatic modes, one with TZ constant and the 
other with TZ vanishing. Since the system of Eqs. H5.1U2|) - H5.104|) is a third-order system 
ther will be 3 independent solutions. However, the Hamiltonian constraint eliminates one, 
the non-adiabatic, i.e. the one leading to 7^' 7^ for A; — > 0. If instead of one scalar 
field there are Ng scalar degrees of freedom the theorem applies as well. The system of 
Eqs. (|5.1()2j) - (|5. 104)1 can be generalized to this case with the result that the number of 
independent equations will be (2A's + 1). Of these {2Ns + 1), 2 are adiabatic and a third 
one is eliminated by the Hamiltonian constraint. Thus, the remaining 2(A's — 1) are non- 
adiabatic modes. 

This perspective was further developed in jl45j (see also jl44p with the purpose of 
discussing the fate of non-adiabatic modes in multi-field inflationary models. As remarked 
in |146| . in a related context, the conservation laws for TZ and its analog for the tensor modes 
of the geometry bears some similarities with the Goldstone theorem |147| I148j normally 
employed in relativistic quantum field theories. The modes for which TZ or hij are constant 
outside the horizon take the place here of the Goldstone bosons that become free particles 
for long- wavelength. 

As a final side remark, it is interesting to mention that the derivation of the present 
results has been conducted in the longitudinal coordinate system. What would happen 
in a different coordinate system, like, for instance, the synchronous? The answer to this 
question is that different gauges may suggest that different quantities are constant outside 
the horizon jl43j . However, this does not imply that different gauges are inequivalent as 
fare as the evolution of physical quantities are concerned. The reason is that the limit /c — > 
may have different meanings in different gauges jl43j . 

5.4 Divergences of specific gauge descriptions 

In the study of the evolution of gravitational perturbations a problem often encountered is 
represented by possible "divergences" in a specific set of variables used to describe a given 
physical system. This often indicates that the description of the fluctuations should be 
performed in an alternative coordinate system where divergences do not necessarily appear. 

The first example is the neutrino isocurvature velocity mode. As it was analytically 
shown in section 2, the metric fluctuations of the neutrino isocurvature velocity mode are 
perfectly regular in the synchronous gauge description (see Eq. ()4.96|) ) while the longitudinal 
fluctuations of the metric are singular in the limit /cr — > (see Eq. ()4. 103(1 ') for a non- 
vanishing fractional contribution of neutrinos to the radiation background, i.e. Ri, 7^ 0. 
The neutrino isocurvature velocity mode is perfectly physical. However, the longitudinal 
coordinate system is not suitable to treat this mode. As discussed before, a natural set of 
gauge-invariant variables seems to be, the Bardeen potentials $ and ^. In the mentioned 
example they will diverge too. 
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Gauge-invariant quantities are useful, however, their physical interpretation has always 
to be discussed within a specific gauge. Consider, for instance, the simple case of a scalar 
field (with quadratic potential) oscillating in a FRW geometry with zero spatial curvature. 
The background equation for (/? follows from Eq. (|5.51|) and it is 

V?" + 2H(^' + m^aV = 0. (5.111) 

The solution of Einstein equations for an oscillating scalar in a quadratic potential are well 
understood (see, for instance, |149[ ll5L)j ). Prom the solution of the system it is easy 
to argue that a'^{pip + p^p) = tp''^ goes to zero with ip'. For conformal times TL < ma the 
solution of Eq. (|5.11H) is simply given by a sinus (or a cosinus) damped by the expansion of 
the Universe as a^/^ so that, averaged over many oscillations, (p^) = 0. This implies that 
the evolution of the scale factor will be with oscillating corrections which also show up 
in Ti. and TC'. If we would insist to study the fluctuations of this system with a decoupled 
equation for the Bardeen potential (see Eqs. (|5.78|) - (|5.79|1 ). we would soon realize that 
because of the divergences in the pump field when the solutions will also diverge. 

Exactly the same occurrence can be verified in the longitudinal gauge where the decoupled 
evolution equation for ip exactly coincides with Eq. (|5.78j) . This problem has been discussed 
in the literature in various frameworks. In jl51l 1152) this patology of Eq. ()5.78|) has been 
pointed out and analyzed. In |15,Sj and jl54| I155| the analysis of the fluctuations in the 
case of multiple oscillating (and interacting) fields has been performed. The knowledge 
of this pathology and of its resolution is now common knowledge in pre-heating studies 

dsni csa csHi cnni cnni nsi cni nnsi umi iibsj. to solve the puzzie, the "trick" is to 

integrate directly the evolution equation for q (see Eq. (|5.66|1 ) the practical observation 
is then that the pump field of q, i.e. z"/z, can be re-expressed, by repeated use of the 
background equations, in a form which does not contain inverse powers of (p' |15H I152j . 
This allows easily to compute, for instance, the curvature perturbation, i.e. TZ, in simple 
cases like such as the ones of quadratic or quartic potential 156, 157 . Again, the rationale 
for the resolution of this pathology can be traced back to a gauge choice. In the longitudinal 
gauge ^/J (and its gauge-invariant generalization, \I') measures the curvature fluctuations on 
zero-shear hyper surf aces. We are then forcing the shear to be zero in a case where this 
choice is questionable. 

Motivated by the divergence of the longitudinal gauge description, let us discuss the 
problem in a different gauge, such as the off-diagonal (or uniform curvature) already intro- 
duced in section 2. The evolution equations of the fluctuations in the off-diagonal gauge 
can be simply obtained from the results of the appendix by setting ip = E = in Eqs. 
HA.14|1 - HA.16|) and ()A.44|) - (|A.46|) . The resulting evolution equations couple then together 
(f>, B and x which are, respectively, the (00) component of the perturbed metric, the off- 
diagonal entry (Oi) of the perturbed metric and the scalar field fluctuation. The relevant 
set of evolution equations can then be written, in the case off a quadratic potential and in 
Fourier space, as: 

n[^' - V^B] = SvrGxV, (5.112) 

Similar problems may arise with quartic potentials and with more general polynomial potentials in the 
vicinity of a local minimum. For instance, in the case of a quartic potential Eq. IS.lllH admits oscillatory 
solutions which are not sinusoidal but rather given by elliptic functions. 
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nicj)' + V^B] + 2{H' + 271^)4) = SirGm^a^x, (5.113) 
Hcj) = 'inGif'x, (5.114) 
(j) = -B' -2nB, (5.115) 



x" + 2nx' - y\ + 



X = 0. (5.116) 



Equations (|5.112j) and (|5.113|) are given, respectively, by the sum and by the difference 
of the (00) and (u) components of the perturbed Einstein equations. Eqs. ()5. 114(1 and 
(j5.115|) are the momentum constraint and the {i ^ j) component of the perturbed Einstein 
equations. Finally Eq. (|5.116|) is the perturbed Klein-Gordon equation (see Eq. (|A.49|) 
where the other equations have been used in order to get rid of the dependence upon the 
metric fluctuations. 

Equation ()5.116|) . being decoupled, can be swiftly solved since its pump field does not 
contain any dependence upon inverse powers of ip'. Then the result can be inserted into 
Eq. (|5.114() that is, again, free of divergences and this leads to the explicit value of (p. 
Finally, direct integration of Eq. (|5.115|) leads to the evolution of the remaining perturbation 
variable, i.e. B. Notice that, in the off-diagonal gauge, q = ax- this can be explicitly verified 
since by repeated use if the background equations, Eq. (j5.116|) can be brought to the same 
form of Eq. (|5.66|) . Thus, the oscillating scalar field is a perfectly physical problem that 
requires, however, a suitable coordinate system to be treated. 

The third example to be discussed is the one of a contracting Universe whose dynamics 
is driven by the kinetic energy of a scalar field. The evolution of the background can be 
easily obtained from Eqs. (|5.49|) - (|5.51|) by setting W{if) = 0. In this case it is easy to show 
that Ti' = —2Ti^. The scale factor will describe then an accelerated contraction for r < 0, 
i.e. a(r) ~ (— r)^/^. Also now, as in the previous two cases, the longitudinal description 
breaks down, in particular, the evolution equation of -0 becomes, in this case 

V-" + 6WV'' - = 0. (5.117) 

The solution for ^p diverges as |r|~^ for r — > 0~. This is the typical situation of the case 
of minimal pre-big bang models in four space-time dimensions. As in the two previous 
situations, one can easily argue that the description of contracting backgrounds is not 
unphysical but should be treated with some care. A detailed analysis in this direction 
performed in the context of pre-big bang models shows indeed [SSI that in the off-diagonal 
gauge the evolution can be treated perturbatively. Setting = in Eq. (|5.116|) and 
recalling that Tl' = —21i?, we have that the evolution of the fluctuation of the scalar field 
fiuctuation is simply 

X" + 2Hx' - V^x = 0, (5.118) 

which also implies that cuvature fiuctuations are constant (up to logarithmic corrrections). 

The example of the accelerated contraction is particularly useful because, from a more 
accurate analysis, it mixes two different problem: 

• the first problem is how (i.e. in which coordinate system) one may describe the 
evolution of the curvature fluctuations in contracting geometries; 
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• the second problem, more difficult and also more intriguing, cocnerns the singularity 
inherent to contracting backgrounds. 

To motivate further the second statement we recall that in the limit r — > 0~, the 
accelerated contraction driven by a minimally coupled scalar field leads to a singularity in 
the curvature. This aspect can be simply appreciated from Eq. ()5.5U|) . implying, in the 
cosmic time description, that H = —AnGcp'^. If the background makes a transition from 
contraction to expansion, H has to change sign and this is impossible to achieve just with a 
minimally couples scalar field. This kind of problem appeared in the literature thanks to the 
development of pre-big bang models (see |1141 1115] ). If the specific evolution of the model 
is not known in all the details, it is rather plausible the the perturbations in the different 
regimes are matched in such a way that the super-Hubble conservation of TZ is enforced. 
Various reprises of these matching conditions can be found in the literature jl661 [TB7j for 
the case of conventional inflationary models. 

The same type of conservation laws obeyed in the ordinary inflationary models was also 
suggested for the pre-big bang case [SH]. In |1161 1117] a specific model for the resolution 
of the singularity has been studied in detail in order to test the arguments put forward 
in [SH]. The model is based on the low-energy string effective action supplemented by a 
non-local dilaton potential whose effect is to regularize the evolution of the gauge coupling 
and of the curvature at short distances. In spite of the fact that the evolution equations 
of the background and of the linearized fluctuations are different, they reduce to the usual 
general relativistic equations for large times away from the maximal curvature regime. By 
studying both analytically and numerically the evolution of the fluctuations, it is found that 
the evolution of the curvature perturbations obeys an equation that is similar to the one 
derived in a general relativistic context, i.e. Eq. (|5.59|) . with an important difference: in 
ordinary general relativity AirGip''^ = Ti^ — Ti' while in case of Ref. |11()| I117j 



in natural string units. The first and the second term are the analog of the relation stem- 
ming from Friedmann equations, while the third term arises from the presence of a non local 
dilaton potential V(^) where Tp is the shifted dilaton, which is invaraint under scale factor 
duality |115j . Thus, when ili? — Ti') = 0, i.e. at the would be singularity, the coefficient at 
the right-hand side of Eq. (|5.59j) does not blow up but it is regular and finite. Thus curva- 
ture perturbations are indeed nearly constant as specifically checked in various numerical 
examples [116) flTT] . The evolution of the tensor modes does also obey a generalization of the 
analog conservation law described earlier in this section. The conservation of super-Hubble 
curavature perturbations has also been checked, recently, in the context of a regularization 
obtained by means of higher derivatives. Similar results can be obtained jl7()| I17H I172j 
(see also, more recently, jl73j ) if the singularity problem is addressed by supplementing the 
low-energy string effective action with higher derivatives terms jl68( I169j . 

The ekpyrotic-cyclic models |119| 112(71 H^lj face situation which is, in some way, similar 
to the one of pre-big bang model and some debates have been recorded in the literature 
all along the last few years. The proponents of the ekpyrotic scenario argue that curvature 




(5.119) 
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perturbations are not continuous across the transition from contraction to expansion jl2Ul 
11741 1175] . If curvature perturbations are conserved on super -Hubble scales, then, as in 
the pre-big bang case the spectrum of fluctuations of the minimal ekpyrotic model will 
be violet thought with a slightly different spectral index jl76| I177| I178| I179j (see section 
6). The proponents of the ekpyrotic scenario, on the contrary, suggest that the curvature 
perturbation is not the conserved quantity that should be used in order to match across the 
transition from contraction to expansion. More specifically the idea is that when branes 
collide curvature perturbations are not conserved across the bounce and they match to the 
nearly scale-invariant fluctuation of the Bardeen potential. This type of arguments has been 
further amplified by some authors |18n[ I181j . 

As a consequence of the debate on what should be continuous across the transition vari- 
ous models of general relativistic bounces have been studied. In general relativity there are 
essentially two ways to get bounces depending on the presence (or absence) of spatial curva- 
ture. In the presence of (positive) spatial curvature Friedmann equations (supplemented by 
a single scalar field with suitable trigonometric potential) may lead to some sort of de Sitter 
bounce jl82[ llS^-i] where the scale factor, in cosmic time, behaves as a{t) cosh(t/to)- In 
the absence of spatial curvature, the only way to get a bounce is to violate the null energy 
condition. In practical terms {p + p) has to change sign. This requirement is impossible to 
achieve with a single scalar field since {p + p) = ^"^^ which is always positive. 

In the context of the ekpyrotic models some variation on the theme of de Sitter bounces 
has been considered for instance in jl84j and in jl85j . The path of the violation of the null 
energy condition has been taken by another set of investigations (see, for instance, jl86j 
and references therein). Just to give an example, in Ref. |186j the bounce is implemented 
through a system of two minimally coupled scalar fields, say fi and ip2- One of the two 
scalars (for instance 992) is a ghost, i.e. it has negative kinetic energy. Thus, according to 
Friedmann equations H will now be the difference if two positive definite quantities, i.e. 
{ip2 — ifi). In this situation H changes sign and a bounce is possible provided ghost fields 
are allowed. Notice that the model studied in jl86j is rather similar to the one discussed in 
jl87j . The conclusions on the evolution of the fluctuations seem, however, to be different. 

6 Normalization and amplification of metric fiuctuations 

Cosmological fluctuations could be either classical or quantum mechanical. Classical and 
quantum cosmological fluctuations have common features hiding physical differences. Dif- 
ferences and analogies can be summarized as follows: 

• in the linearized approximation, classical and quantum fluctuations obey the same 
evolution equations; 

• classical fluctuations are given once forever (on a given space- like hypersurface); 

• quantum fluctuations keep on reappearing all the time during the inflationary phase 
because of the zero-point fiuctuations of the various fields (metric inhomogeneities, 
perturbations of the inflaton, fluctuations of some spectator field); 
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• the evolution and normalization of quantum mechanical fluctuations can be consis- 
tently described both in the Heisenberg and in the Schrodinger representations. 

A natural question to ask is wether the fluctuations affecting the SW effect are classical or 
quantum mechanical. The answer to this question cannot be given in general terms but it 
depends on the model describing the early evolution of the Universe. 

Consider, for simplicity, the framework of conventional inflationary models driven by a 
single scalar degree of freedom. The minimal duration of inflation is between 60 and 65 
efolds. The rationale for such a statement is that, today, the total curvature of the Universe 
receives a leading contribution from the extrinsic curvature and a subleading contribution 
from the intrinsic (spatial) curvature. The ratio between the intrinsic and extrinsic curvature 
(see, for instance, Eq. (|2.6j) ) goes as l/a?. During an epoch of decelerated expansion (i.e. 
a < 0, a > 0) such as the ordinary radiation and matter-dominated phases, l/a? can become 
very large. The role of inflation is, in this context, to make very minute at the end 
of inflation, so that it can easily be of order 1 today. The minimal duration of inflation 
required in order to achieve this goal is about 60-efolds. 

If the duration of inflation is minimal (or close to minimal) classical fluctuations, which 
were super-horizon sized at the onset of inflation will be affected neither by the inflationary 
phase nor by the subsequent post-inflationary epoch and can have computable large scale 
effects like, for instance, the so-called Grishchuk-Zeldovich effect jl88| ITKflj . If the fluctua- 
tions are classical, there are, virtually no ambiguities in normalizing them: it is sufficient 
to assign the values of the various inhomogeneities over a typical scale and at a given time. 

When the duration of inflation is much longer than 60-efolds, the large scale fluctuations 
are probably all of quantum-mechanical origin, at least in the case of inflationary models 
driven by a single inflaton field. Quantum-mechanical fluctuations result from the zero- 
point energy of the metric inhomogeneities (and possibly of other quantum fields) present 
during the inflationary epoch. The distinction between classical and quantum cosmological 
fluctuations, can be appreciated by comparing Figs. [51 and ^ 

In both figures, a schematic evolution of the logarithm (in ten basis) of the Hubble rate 
(in Planck units) is reported both for the case of non-minimal duration of the inflationary 
phase (Fig. ^ and in the case of minimal duration (Fig. On the horizontal axis 

the logarithm of the scale factor is reported. The inflationary region corresponds, in both 
figures, to the quasi-flat plateau lasting 82 efolds (in Fig. ^ and 65 efolds (in Fig. |^. 
These numerical values are purely illustrative. 

If the duration of inflation is non-minimal a classical fluctuation present very close to the 
beginning of inflation will make the second crossing of the Hubble radius in the far future 
(see Fig. 01 thin full Une). Today, H/My> ~ 10"^^ while during inflation H/M^ ~ 10"^- 
10~^. Then it is clear that a hypothetical classical fluctuations crossing the horizon close 
to the beginning of an inflationary phase (lasting much more than 60 efolds) will not cross 
the horizon the second time prior to the present epoch (see Fig. ISJ. Conversely (see Fig. |1} 
if the duration of inflation is close to minimal the latter occurrence is not excluded. In the 
following, after a detailed discussion of the standard lore, some unconventional approaches 
to the normalization of quantum fluctuation will be reviewed. 
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Figure 4: The same quantities illustrated in Fig. [S] are reported in the situation when the 
flat plateau lasts for about 65 efolds. 

6.1 The standard lore for the normalization 

The evolution of the metric fluctuations can be described either within the Schrodinger 
description, or within the Heisenberg description. In the Schrodinger description, the evo- 
lution can be pictured as the spreading of a quantum mechanical wave-functional. The 
initial wave-functional will be constructed as the direct product of states minimizing the 
indetermination relations for the different harmonic oscillators forming the quantum field. 
The initial vacuum state will have zero momentum and each Fourier mode of the field will 
evolve into a (two mode) squeezed quantum state where the indetermination relations are 
still minimized but in such a way that one of the two canonically conjugate operators will 
have a variance much larger than the quantum limit, while the other canonically conju- 
gate operator will have a variance much smaller than the quantum uncertainty limit. The 
description of the evolution of cosmological fluctuations by means of the squeezed states 
formalism is rather instructive and has been exploited by various authors jl9()| I19H I192j . 
In this formalism, some inspiring connections with quantum optics emerge (see for instance 
[90]). Squeezed states are a powerful generalization of the concept of coherent state pT] . 
As the pumping action of a background electromagnetic field (a laser) is able to produce, 
under some circumstances, squeezed states of photons, the pumping action of the classical 
gravitational field (for instance the curvature) is able to produce squeezed states of gravi- 
tons (for the tensor modes) or of phonons (in the case of the scalar modes). In the context 
of the squeezed states formalism a natural definition of a coarse grained entropy of the 
quantum fluctuations emerges. This entropy, in the analysis of jl931ll9^ . is associated with 
the spreading of the wave-functional in an appropriate basis. The coarse graining procedure 
can be performed within different basis with consistent results jl951 11961 11971 H^Sl 11991 000] . 
The concept of entropy of cosmological fluctuations leads naturally to the problem of the 
transition from the quantum to the classical regime which can be in turn related to the 
decoherence of the density matrix of the cosmological fluctuations. The Schrodinger de- 
scription proves also very useful for the self-consistent evolution of the inflaton and of its 
fluctuations [201] . 
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In fully equivalent terms, the evolution of the fluctuations can be described within the 
Heisenberg representation. This is the description which will be adopted for the purposes 
of the presentation. The Hamiltonians of the scalar and tensor modes of the geometry will 
then be minimized for r — oo, which is a physical limit, not a mathematical one. In fact, 
inflation cannot last indefinitely in the past. Since de Sitter space-time is not geodesically 
complete in the past^'' , there is no reason to assume that r should indeed be going to — oo. 



6.1.1 Large-scale power spectra of the tensor fluctuations 



Consider the Hamiltonian given in Eq. 1)5. 15() . dropping, for simplicity, the tildes in the 
momenta: 



vr^ fi^ + idifiY' 



(6.1) 



where = a/i^Mp/y^ and identically for the other polarization. 

After imposing the commutation relations for the canonically conjugate quantum fields 
(units h = 1 are adopted) 

[fiix,T),n{y,T)]=i6^^\x-y), (6.2) 

the operator corresponding to the Hamiltonian 1)6. 1|) becomes, then, 

1 
2 



H{t) 
tum fi 

fi{x,T) 



d^x 



a 



(6.3) 



In Fourier space the quantum fields can be written as 

1 



7r(x, r) 



2(27r)3/2 
1 



^ „—ik-x I ±.] Ak-i 
k 



(6.4) 



2(27r)3/2 

Inserting now Eqs. (|6.4|) into Eq. ()6.2() and demanding the validity of the latter implies 
the following canonical commutation relations for the Fourier components of the quantum 
operators: 

[Afc(r),4(T)]=i5(3)(^_p-^, 

[(i\{r),iT^{r)]=i6^^\k-p), 
[Afc(r),vrp^(T)]=i<5(3)(^ + p-), 

[Al:W,4(T)]=i<5(3)(fc+p-). (6.5) 

Inserting now Eq. (|6.4|) into Eq. (|6.3|1 we get the Fourier space representation of the 
quantum Hamiltonian ^^: 

1 (/^fe/^? + 



H(t) 



(Pk 



a 



(6.6) 



^"The analysis of the past geodesic (in)completeness of conventional inflationary models has been discussed, 
in a series of papers, by Borde and Vilenkin PU^HTHIIWI and by Borde, Guth and Vilenkin ;208' . 

Notice that in order to derive the following equation, the relations /i^ ^ = /tg and tt^ ^ = vrg should be 
used . 
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In the Heisenberg representation the field operators obey: 



ITT 



(6.7) 
(6.8) 



Using now the mode expansion ()6.4|) and the Hamiltonian in the form ()6.6|) the evolution 
for the Fourier components of the operators is 



a 
a 



The general solution of the system is then 

T^k{r) = aj:{To)gk{T) + aK{To)gl{T), 



where the mode function fk obeys 



fk = 0, 



(6.9) 
(6.10) 

(6.11) 
(6.12) 

(6.13) 



and Qk = f'k- In the case when the scale factor has a power dependence, in cosmic time, as in 
Eq. H5.64() , the scale factor will be, in conformal time a(r) = {—t/ti)~^ with /3 = p/{p—l). 
The solution of Eq. 1)6. 13(1 is then 



fk{r) 
9k{T) 

where x = kr and 



2k 

f'k 




H 



(1) ( 



-x) + 



:i - 2^l) 



2{-x) 



The functions H. 



AT 



J„ 
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(M+l/2)7r 



fl = P + 



(6.14) 
(6.15) 

(6.16) 



-x)+ iY^{-x) is the Hankel function [20111011 of ^rst kind 

and the other linearly independent solution will be H^^ = H^^^*^. Notice that the phases 
appearing in Eqs. (|6.14|1 and (|6.15|1 are carefully selected in such a way that for r ^ — oo, 
fk^e-'^yV2k IMI- 

In the case of more general inflationary potentials the solution of Eq. 1)6. 13(1 is also 
reasonably simple and it has the same form given in Eq. (|6.14)) and (|6.15)) . However, this 
time, using the slow-roll formalism, the index /i will be given by Eq. 1)5. 73() . 

In Eqs. ()6.11|) and 1)6. 12(1 the operators a^{TQ) and a_g(To) annihilate the vacuum, i.e. 
the state minimizing the Hamiltonian (|6.6|1 for tq — > — oo which is chosen, as anticipated, 
close to the onset of inflation. The time tq is the same for different modes k. The second 
remark is that, the vacuum, i.e. the state such that a_^(To) ag(To)|Ooo) = 0, has zero total 
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momentum. This can be understood in terms of a process of pair production from the 
vacuum triggered by the pumping action of the gravitational field. Inserting Eqs. 
and (|6.12|) into Eq. 1)6. 4() the following Fourier expansions can be simply obtained 



, r) 

Tf{x, t) 



1 



(27r)3/2 

1 

(27r)3/2 

ik-x 



a.(ro)/fc(r)e-'=- + at (ro)4*(r)e*'=- 



«,Kro)<7fc(r)e-'=- + 4(ro)5|(r)e^^- 



(6.17) 



since, upon integration, a_re = are 



-ik-x -t ik-x ^ at gifc-x_ 

k 



-k k -k 

The two point function for the tensor fluctuations can then be obtained by computing 



(0_oo|/i(x,r/)/i(y,7?)|0_ 



24 



2 ^—ik-f 



(27r)3 a2(r) 

where r = {x — y). After angular integration, the previous expression becomes 

sin kr 



where 



(O_oo|/i(x,?7)/i(y,?7)|0_oo) = 
24 k^ 



dln/c Vt 



kr 



a2(r) 27r2 



(6.18) 



(6.19) 



(6.20) 



is the tensor power spectrum. Taking the small argument limit |204j of the Hankel functions 
appearing in Eq. l|6.2Ujl . the square root Eq. (|6.20j) becomes 



1/2 



_2f_ 

3/2 



VT' 



2 ' 



(6.21) 



where Eq. (|6.2U|) has been multiplied by 2 to account for the 2 poalrizations. 

Since, by definition, ki = \aiHi\ = /3/ri, and bearing in mind that /5 = p/{p — 1), Eq. 
(Unn leads to 



Mp VA;i 



7r3/2 ^V2(p-1 



p — 1 
p 



p/{p-l) 



(6.22) 
(6.23) 



In the case of generic slow-roll inflation we will have 



V. 



T -^7^r(/.)(i-e)'^- 

7r3/2 \MpJaH=k 



where /i is now determined by Eq. 1)5. 73(1 . 



(6.24) 



Notice that as a consequence of the different (conventional) definitions of the perturbed line element 
Ref. |132| arrives at a power spectrum which is 1/2 of our result (see also previous footnote). 
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In the study of slow-roll dynamics the tensor spectral index is usually defined as 



UT = = -2e - (3 + ci)e2 - (1 + ci)er/, (6.25) 

where the second equality follows from the expansion of the obtained result in terms of 
the slow-roll parameters which are small deep within the inflationary phase; the constant 
is biven by ci = 0.0814. It is an exercise to perform directly the derivative of InT^T with 
respect to In/c. To obtain the correct result it is useful to employ Eqs. (|5.68|) and (|5.69jl . 
For instance, recalling that k = Ha, it can be easily shown that 

dink _d\D.H ^d\na _lda M^H ^^^^^ 
dip dip dip a dip 2 if^^ 

where the comma followed by ip denotes a derivation with respect to ip and where the second 
equality follows from Eq. (|5.69|) . Introducing the explicit definition of slow-roll parameters 
it can be also shown that 

^ = ^ie-l). (6.27) 
dip V2e 



6.1.2 Large-scale power spectra of the scalar fluctuations 



The calculation of the scalar power spectra follows exactly the same algebra discussed in the 
case of the tensor modes. In full analogy with the calculation of the tensor power spectrum 
the Hamiltonian given in Eq. 1)5. 96() can be chosen. Promoting the normal mode q and its 
conjugate momentum to quantum operators, the relevant commutation relations are now: 



(6.28) 



Recalling the relation between the normal modes and the gauge-invariant curvature fiuctu- 
ations, obtained in Eq. ()5.65|1 . we can then write the scalar two-point function as 



(O_oo|7t(r,x)7t(r,y)|0_ 



1 



d^k 
(2^ 



l/sfc(r)|^e 



2-ik-r 



(6.29) 



where /^(t) are the mode functions pertaining to the scalar problem and obeying the 
equation 



fsk + 



s k 



0. 



(6.30) 



While in the power-law case the solutions of Eq. (|6.3U|) are exactly the same as Eq. (|6.14() , 
in the in the generic case /s t will be 



/s k{T) 

N 

where, as usual, x = kr. 



A4 



2k 



^gl(i.+l/2K 



3 + e + 2r? 
2(1-6) • 



(6.31) 
(6.32) 
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After having performed the angular integration in Eq. H6.29|l the scalar power spectrum 
is then 



1/2 
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/^2V(p-i)^ 
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(6.33) 
(6.34) 



in the case of power-law inflation and 
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P / aH=k 



(6.35) 

for the general slow-roll case whose related pump field has been computed in Eq. (|5.74)1 . 



Notice that, within our notations "Pr 



1/2 ,^1/2 



T 



2V2e. 



As in the case of the tensor power spectrum, the spectral index is conventionally written 
to lowest order as^^ 

^^"^^ - 2ry. (6.36) 



1 



dink 



-4e 



By comparing the scalar and tensor power spectra with the tensor spectral index in 
the slow-roll approximation we find that nx = —{1/4:)Vt/'Ptz- These types of relations are 
often called consistency relations and play a relevant role in the attempts to reconstruct the 
inflaton pootential in single field inflationary models. Notice that some authors argue that 
it is possible to redefine the power spectra I133j in such a way that the ratio of the newly 
defined spectra indices equals exactly e (without the numerical pre-factor obtained within 
the present conventions). 



6.2 Transplanckian problem(s)? 

Consider, for simplicity the case of power inflation parametrized in terms of the conformal 
time coordinate r as a(r) ~ (— t)~^ with (3 > 0. If the standard normalization prescription 
is interpreted in strict mathematical terms, then it will happen that if the initial normal- 
ization time To is sent to — oo, a given physical frequency at the time tq, oj = k/a{To) 
will become much larger than the Planck mass, or as often emphasized, transplanckian. 
A common theme of various investigations in this direction of research is the observation 
that in the "transplanckian" regime the precise description of the evolution of the metric 
fluctuations could be foggy. These statements leave the room to various proposals which 
can be summarized as follows: 

• the dispersion relations can be modified by transplanckian effect; 

• the indetermination relations are modified; 

Notice that slightly different definition of the slow-roll parameters are possible. Our results follow the 
conventions previously established in Eq. (15.6711 and are consistent with the ones of |132| . Some authors like 
to define, for instance, the second slow-roll parameter as 77 = MpV^ipip/V. In terms of our definitions we will 
then have rj — e — rj. This also implies that ris = 1 — 6e -f 2rj. 
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• transplanckian physics induces the presence of a new fundamental scale M smaller 
than the Planck mass; observable effects in the scalar and tensor power spectra can 
be expected. 

The modifications invoked in the dispersion relations, as often acknowledged by the authors, 
have always some ad-hoc feature. This possibility was extensively investigated, in a series 
of papers, by Brandenberger and Martin j2U9l I21UL I211j as well as by other authors [2121 
12131 12l4l 12151 [^TBl 12171 f218j . The observation may be summarized by saying that the wave- 
number of the fluctuation is modified in such a way that k"^ — > k'^g{k,T) |2U9l I21U1 1211j . 
Such modifications should be derived from a suitable variational principle |215j (but this is 
not completely obvious) and they can well be non-linear |213j . Some effort has been made in 
order to justify the existence of modified dispersion relations in various frameworks ranging 
from analog modifications arising in black-hole physics [215j to quantum- Poincare algebras 
applied to a cosmological setting [21 4j , to non-commutative geometries |21f)| I217j . 

The modifications of the indetermination relations are sometimes motivated by the 
scattering of superstrings at Planckian energies |219l I22UI 12211 I222j . However, in that 
context the modifications in the dispersion relations occurs, in critical dimensions ( 26 in 
the case of the bosonic string); the "position" operator is an impact parameter of two 
colliding strings. In the transplanckian context, on the contrary, the modifications of the 
indetermination relations are studied in a four-dimensional context. Tachyonic instabilities 
then occur (this is also a problem, in some cases, when dispersion relations are modified). 
This approach has been followed by various authors [2231 12241 12251 1226j . 

One possible approach to the problem would be to assume that there is a fundamental 
scale M < Mp and that fluctuations should be normalized as soon as they "exit" from the 
physical regime characterized by the scale M. This aspect is partially illustrated in Fig. 
El where with the dotted line the new fundamental scale is reported. According to these 
types of proposals [2271 12281 1246j the modes of the quantum field should be normalized (at 
a finite conformal time tq) as soon as the physical frequency (with dashed line in Fig. ^ 
crosses the dotted line. This crossing clearly occurs at different times for different comoving 
frequencies. In the case of power-law inflation, a given physical frequency u;(r) = k/a{T) 
will "cross" the scale M at a characteristic time To{k) determined as 

where a(r) ~ (— r/ri)^^. The second equality in the second equation follows from the rela- 
tion between the cosmic and conformal times in power-law backgrounds (see, for instance, 
before Eq. (|6.14|) ). In the expanding branch of de Sitter space-time /? = 1 and p — > oo. 
Consequently /ctq is constant and roughly given by M/H where H is the Hubble rate during 
the de Sitter phase. Now the claim is that j227l I228j this prescription for setting the ini- 
tial condition for quantum fluctuations produces observable corrections to the tensor power 
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spectrum of the form 

(6.38) 

where is a numerical constant; Vt is the tensor power spectrum obtained through the 
standard normahzation prescription discussed earUer in this section. Oscillating terms mul- 
tiplying the ratio (H/M)"' have been neglected (see below, however). The ellipses in Eq. 
1)6. 38() stand for terms which are of higher order in H/M. Concerning the parametrization 
given in Eq. (|6.38)) two comments are in order: 

• it is clear that if the prescription 1)6. 37(1 is applied the correction to the power spectrum 
has to be in the form (|6.38|) : the reason is that the power spectrum will now depend 
upon the new scale xq = /ctq ~ M/H ^ 1 that becomes the argument of Hankel 
functions whose limit for xq ^ 1 must indeed produce a result of the type (|6.38|) : 

• the power 7 is crucial for the possible observational relevance of transplanckian effects 
so it becomes crucial to understand what controls the power 7. 

Some authors claim that 7 can be 1 j2271 1225] . In this case the ratio H/M will be larger 
than H/Mp that we know can be as large as 10~^-10~^. If 7 = 2 or even 7 = 3 the 
correction is not relevant. Different arguments have been put forward for values 7 > 1 
|23n| I23H I232| I233| I234j . Notice that often the discussions of these effects take place in 
pure de Sitter space (where, as already pointed out, scalar modes would be absent). In this 
context different vacua can be defined. They are connected by unitary transformations and 
generically indicated as a vacua |235| I236| I237| I238| I239j . Flat space has a global time- like 
killing vector. This allows to define a Hamiltonian whose minimization defines a "vacuum", 
i.e. the lowest energy eigenstate. In de Sitter space there is no global time-like killing 
vector. So a globally conserved energy cannot be defined. It is however still possible to 
find states which are invariant under the connected part of the isometry group of de Sitter 
space, i.e. 5*0(4, 1). In the case of a free scalar field there is an infinite family of invariant 
states called a vacua since they can be distinguished by a single complex number a. 

6.2.1 Minimization of canonically related Hamiltonians 

In the following the roots of the corrections to the power spectrum will be understood 
in terms of a much more mundane feature of the theory of cosmological fluctuations, i.e. 
the possibility, already discussed in section 5, of defining different Hamiltonians related 
by canonical transformations. More specifically it will be demonstrated that the value of 
7 given in Eq. (|6.38|) depends on which Hamiltionian one wishes to be minimized at the 
moment when the given physical frequency cosses the new fundamental scale M as explained 
in Eq. (|6.37j) . By averaging the energy- momentum pseudo-tensor of the fluctuations of the 
geometry over the state minimizing a given Hamiltonian at To(k), the energetic content of 
the fluctuations can be estimated. The result of such a calculation must be smaller that the 

^*Here we discuss the tensor case since, for exact (expanding) de Sitter space only the tensor modes 
are excited. However, the same discussion hold also for the scalar power spectra in the quasi- de Sitter 
space-time. 
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energy density of the background geometry. This chain of calculations has been recently 
presented in j24Uj and j241j and will be briefly reviewed here. 



The first step is to show that, depending on the which Hamiltonian is minimized at 
ro(/c) (defined as in Eq. 1)6. 37() ) a different power of 7 is obtained in Eq. I|6.38() . In order 
to be specific consider the case of conventional power-law inflationary models. We will flrst 
derive the result for the cases of the Hamiltonians ()5.88|) (for the scalar fluctuations) and 
(|5.11j) (for the tensor fluctuations). Then the results of the minimization of the other two 
Hamiltonians (i.e. Eqs. H5.94() - (|5.96|) and Eqs. (|5.13|) - H5.15|) 1 will be given. 



Equation ()5.88|) implies that the canonical field is TZ, i.e. the curvature perturbation. 
The canonical momentum is the density contrast as discussed in Eq. (|5.92|) . The Hamilto- 
nian operator will be 



H{t) 



(Tx 



(6.39) 



Repeating the same procedure outlined earlier in this section, Eq. (|6.39l) can be written as 



The Hamiltonian 1)6. 4U() can be written at to(/c) as 
H{ro) = \ 

where 



(Fkk 



Qdro) 



1 



iz{To)kTZj:{To) 



(6.40) 



(6.41) 



(6.42) 



obeying [Qg, Q^] = 6^^\k—p). Consequently, the state minimizing (|6.4U|1 at tq is the one 
annihilated by Qg, i.e. 

g,Kro)|0«) = 0, Q„fc(ro)|0«)=0. 



(6.43) 



The specific relation between field operators dictated by ()6.43|) provides initial conditions 
for the Heisenberg equations 



in' = [n,H], 

The full solution of this equation can be written as 



n^ir) = ag(ro)A(r) +a^.(ro)/,*(r), 
%(r) = ar{To)gkiT) + ATo)9*k{r), 



(6.44) 

(6.45) 
(6.46) 



where, recalling the explicit solution of the equations in the case of the exponential potential 
(j5.64|) and defining x = kr 
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(6.47) 
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satisfy the Wronskian normalization condition 



fk{T)gUT)-ft{r)gk{T)=i. 



(6.48) 



Notice that the mode functions of Eq. ()6.47() are different from the ones previously wtitten 
earlier in this section. The creation and annihilation operators appearing in 1)6. 12() are 
defined as 



ar(To) 



1 



,,--M{ro) + ikzyi{T,)]Q^^{T,) - [gUr,) - ikz'^fl{T,)]Q^ ^(tq)}, 
z^y Ik ^'^ 

-^^{[9k{ro)-ikzlfk{To)]Ql^{To) - [gkiTo)+ikzlfkiTo)]Q^{To)}{6A9) 



So far two sets of creation and annihilation operators have been introduced: the operators 
Q^{tq) and the operators d^(ro). The state annihilated by Q^{tq) minimizes the Hamilto- 
nian at tq while the state annihilated by a(To) does not minimize the Hamiltonian at tq. 
It is relevant to introduce these operators not so much for the calculation of the two-point 
function but for the subsequent applications to the back-reaction effects. In fact, in the 
standard approach to the initial value problem for the quantum mechanical fluctuations, 
the initial state is chosen to be the one annihilated by a^(''"o) for tq — oo. 

The Fourier transform of the two-point function, 
(0«,ro|7t(f,r)^(y,r)|ro,0«) = 



dk— sin kr 
k kr ' 



\x-y\, 



(6.50) 



can now be computed, and the result is 



|gfc(-^o)l' 



+ k'z{Tof\fk{To)\' 
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(6.51) 



The explicit form of V-ji and Vt can be obtained by inserting Eqs. 1)6. 47() into Eq. 
(|6.51j) . The results should be expanded for |x| = fer ^ 1 and for |xo| = krQ » 1. While 
|fcr| measures how much a given mode is outside the horizon. 
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\xo\ 



defines the moment at which the given mode crosses the scale M. Notice that Eq. 1)6.52(1 
has exactly the same content of Eq. ((6. 37(1 

To have the explicit form of the power spectrum, Eq. ()5.1U7() should be expanded for 
kr <^ 1 and |A;ro| ^ 1 and the result is 
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(6.53) 
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1/2 

where V-jl is the power spectrum of the scalar fluctutations aheady obtained within the 
standard lore [i.e. Eqs. (HM - dCTl) ]. 

The results obtained in the case of Eq. 1)5. 88(1 can be also generalized to the case of the 
tensor modes decsribed by the Hamiltonian 1)5. In this case the calculation is identical 
with the only difference that z{t) is replaced by a(T) and IZ is replaced by h. The result 
will be : 



V. 



1/2 



1 + 



p sin [2xq + ptt/ {p - 1) 



1 



1/2 

where, now, V^ji is given by Eqs. 



2{p 



XQ 



(6.54) 



In Eqs. 1)6. 53(1 and (|6.54|) . on top of the standard (leading) terms there is a correction 
that goes, roughly, as I/xq ~ H^x/M where, as discussed in Eq. 1)6. 52|) . H^^ denotes 
the Hubble parameter evaluated at the moment the given scale crosses M. If M ~ Mp, 
Hex/M ~ 10~^. This is the correction that would apply in the scalar power spectrum if 
quantum mechanical initial conditions were assigned in such a way that the initial state 
minimizes (|5.88j) . 



Having discussed in detail the results for the case of (|5.88j) the attention will now be 
turned to the case of Eqs. ()5.94|) and 1)5. 96() whose tensor analog are Eqs. 1)5. 13() and (|5.15j) . 
In this case the whole procedure of canonical quantization can be repeated with the crucial 
difference that the state minimizing the quantum version of (|5.94jl will not be the same 
minimizing 1)5. 88(1 . If the Hamiltonian 1)5.94(1 is minimized the result is 



1/2 

7e,T 



1 



p COS [2xo + ptt/ {p - I)] 



A{p - 1) 



(6.55) 



where the subscript denotes either the scalar or the tensor power spectrum since the cor- 
rection, in a power-law inflationary background, is the same both for scalars and tensors. 

A comparison of Eqs. (|6.54() and l|6.55() shows two important facts. The first is that the 
leading term of the spectra is the same. This phenomenon simply reflects the occurrence 
that different Hamiltonians, connected by canonical transformations, must lead to the same 
evolution and to the same leading term in the power spectra. The second fact to be noticed 
is that the correction to the power spectrum goes as 1/xq in the case of ()6.55|) . This 
correction is then much smaller than the one appearing in (|6.54|) . If M ~ Mp then the 
correction will be, in the exact de Sitter case, O(10~^^), i.e. six orders of magnitude smaller 
than the correction appearing in Eq. ()6.54|) . 



Finally the case of the Hamiltonians (|5.96j) and ()5.15(1 will be examined. Equation 
(|5.96() can be minimized following the same procedure as already discussed in the case of 
Eqs. (|5.88|) and ()5.94|) . However, again, the state minimizing the Hamiltonian of Eq. (|5.96j) 
will be different from the states minimizing the Hamiltonians of Eqs. (|5.88() and (|5.94|) . 
The result for the power spectra will then be 



' TZ,T 



V. 



1/2 



1 + 



p(^p — 1) sin [2xo + p7r/ (p — 1) 



4^0 



3.56) 



(p-l)2 

In Eq. (|6.56|) the correction arising from the initial state goes as l/a^g ^-iid, again, if M ~ Mp 
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it is 0(10^1^), i.e. 12 orders of magnitude smaller than in the case discussed in Eqs. ()6.53() 
and 

6.2.2 Back-reaction effects 

In order to select the correct Hamiltonian in a way compatible with the idea of assigning 
initial conditions when a given physical frequency crosses the scale M it is desirable to 
address the issue of back-reaction effects. The energetic content of the quantum-mechanical 
state minimizing the given Hamiltonian should be estimated and compared with the energy 
density of the background geometry. The back-reaction effects of the different quantum- 
mechanical states minimizing the Hamiltonians will now be computed. Without loss of 
generality, the attention will be focused on the tensor modes of the geometry. The advantage 
of discussing the gravitons is that they do not couple to the sources and, therefore, the form 
of the energy-momentum pseudo-tensor is simpler than in the case of the scalar modes 

The appropriate energy-momentum tensor of the fluctuations of the geometry will be av- 
eraged over the state minimizing a given Hamiltonian at To{k) and the result compared with 
the energy density of the background geometry, since he energy density of the perturbations 
cannot exceed that of the background geometry. 

The energy density of the tensor inhomogeneities can be computed from the energy- 
momentum pseudo-tensor, written, for simplicity, for one of the two polarizations: 

(-^0°) = ^{{h'h + hh')) + + [dm). (6.57) 

where (...) denotes the expectation value with respect to a quantum mechanical state min- 
imizing a given Hamiltonian and h denotes the field operator corresponding to a single 
tensor polarization of the geometry. 

If, in Eq. (|6.57|) . the average is taken with respect to the state minimizing the Hamilto- 
nian H5.11() . then the energy density is, in the case of exact de Sitter space [241 j : 

.iJU.ro)-^i|^[.S.o(^)] .^(f )'[l.0(g)]. (0.58) 

Since, as already discussed, |rco| = M/H » 1, in the case of de Sitter space, the back- 
reaction effects related to the state minimizing the first Hamiltonian are then large. Recall, 
in fact, that the energy density of the background geometry is 0{H'^Mp). Hence, if M ~ Mp 
the energy density of the fluctuations will be of the same order as that of the background 
geometry, which is not acceptable since, if this is the case, inflation could not even start. 
Let us now turn our attention to the case of the state minimizing the second Hamiltonian, 
i.e. Eq. 1)5. The expectation value of the energy- momentum pseudo-tensor will now 
be taken over the state minimizing the second class of Hamiltonians. Applying the same 
procedure as before the averaged energy density is smaller than the result obtained in Eq. 
(iraj) by a factor {H/Mf. 
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Finally, in the third and last case, we have to average the energy-momentum pseudo- 
tensor over the state minimizing the Hamiltonian (|5.15|) . In this case the result is |/i41| 



2567r2 



27 




(6.59) 



i.e. even smaller than the result discussed in the case of the second Hamiltonian. In this 
case, the averaged energy density is much smaller than that of the background. 

The results reported in j241j and summarized above suggest the following reflections: 

• it is physically interesting to set quantum mechanical initial conditions as soon as a 
given physical frequency crosses a fundamental scale M; 

• to compute "observables" it is crucial to determine which Hamiltonian is minimized 
as soon as the given physical frequency crosses a fundamental scale M; 

• large (observable) effects seems to be ruled out if one accepts that the energy density of 
the fluctuations should be smaller than the energy density of the background geometry. 

Similar considerations, but through different arguments, appeared in |244[ I245[ I246j . It 
is here appropriate to mention that the authors of Ref. [2471 1248| claim to reach more 
optimistic conclusions concerning the back-reaction effects associated with particular classes 
of transplanckian effects. In particular, in Ref. |247 l the analysis has been performed 
just in terms of the scalar field fluctuation (that is not gauge-invariant for infinitesimal 
diffeomorphisms) . It would be interesting to repeat the same analysis in terms of the tensor 
modes of the geometry. 

6.3 Spectator fields 

It is rather plausible that the early Universe is not described by a single scalar degree of 
freedom. On the contrary it seem conceivable that various fields are simultaneously present. 
In the remaining part of the present section the attention will be focused on the possible 
effects of the spectator fields, i.e. fields that are not source of the background geometry but 
whose fluctuations can be amplified with scale-invariant spectra. 

Thanks to some of the considerations developed in section 4 and 5, in the presence 
of some spectator field, a non-adiabatic fluctuation of the pressure density is likely to be 
generated in the system. There are different (but related) concrete examples of such a 
dynamics. The first obvious case is the one where the dynamics of the background is driven 
by a scalar degree of freedom (for instance the dilaton, in pre-big bang models |114| 111?)] ). 
While the Universe inflates, some other spectator fleld (that is not source of the background 
geometry) gets amplified with a quasi-fiat spectrum. The role of spectator field is played, 
in pre-big bang models, by the Kalb-Ramond axion jll5j . When inflation terminates, the 
large-scale fluctuations in the curvature are vanishingly small owing to the steepness of 
the dilatonic spectrum. After the dilaton/inflaton decay, the Universe will presumably be 
dominated by radiation. However, owing to the presence of the spectator field, the dynamics 
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of the inhomogeneities will be much richer than the one of a radiation-dominated Universe. 
Since the radiation peculiar velocity and density contrast are present together with the 
fluctuations of the spectator field, a non-adiabatic fluctuation of the pressure density, 5pnad 
can be expected. Since 6pna,d sources the evolution of the curvature perturbations, the 
quasi-flat spectrum of the spectator/axion field may be converted, under some dynamical 
requirements, into curvature fluctuations after axion's decay. 

Another possible example is the one where inflation occurs at relatively low curvature 
scales, i.e. H <^ 10~^ Mp. In this case one can imagine the situation where the fluctua- 
tions of a scalar field (that is light during inflation and later on decays), can be efficiently 
converted into adiabatic curvature perturbations. 

In |249j the possible conversion of isocurvature fluctuations into adiabatic modes was 
investigated in a simple set up where on top of the inflaton field there is only an extra 
spectator field. The chief objective of Ref. j249j was not the possibility of converting 
isocurvature into adiabatic modes; on the contrary, there was the hope that fluctuations of 
a pseudo-scalar spectator could indeed give rise to isocurvature modes j25Uj after the decay 
(taking place after baryogenesis in the model of |249j ). The possibility of conversion was 
also briefly mentioned in Ref. j251j : the main goal of j251j was however the analysis of 
possible non-gaussian effects arising in isocurvature models. 

In j252.! (see also |253j ) the role of spectator fleld was played by some super-string 
moduli that could be amplifled during a conventional inflationary phase. The authors 
correctly pointed out the possibility of obtaining, after modulus decay, a correlated mixture 
of adiabatic and isocurvature fluctuations. In I254j the conversion of isocurvature modes 
has been analyzed in the context of pre-big bang models. In fact, while the spectrum the 
dilaton and of the graviton is rather steep, the spectrum of axionic fluctuations may be 
rather flat or even red. The axionic fluctuation amplifled during the pre-big bang phase 
can then be converted into adiabatic modes. In Ref. |255j the authors called the spectator 
field curvaton and provided further support to this idea. 

Consider, for simplicity a scalar field a whose potential is quadratic, i.e. 2V{a) = m?a'^. 
The field a, at least initially, is not source of the background geometry. After inflation its 
evolution can be written, 



a 



2 r ^'2 



a'' 



(6.60) 
(6.61) 



dV 

a" + 2na' + a^— = 0, (6.62) 
oa 

p'^ + AHpr = 0. (6.63) 

where 

Pa = ^,+V, P. = ^,-V. (6.64) 

Since at the onset of the radiation-dominated epoch the curvaton is not dominant, its energy 

density is always smaller than the energy density of the background geometry. During 
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radiation Tia is constant; from Eq. (|6.62|) it can be easily derived that 



= (6.05) 

The factor of 5 appearing in Eq. H6.65() comes from the requirement that the potential 
term dominates in Eq. (|6.62|) . Thus the approximate evolution will follow the equation 
TLa' = —Ca?V^a where V^a denotes the derivation with respect to a and C is a constant 
to be determined inserting the approximate equation into the full equation and neglecting 
terms with higher derivatives of the potential, i.e. V^o-o-- This requirement leads to C = 1/5. 
Notice that the factor 1/5 depends on the specific evolution of the dominant component of 
the background geometry. If the dominant component of the background is characterized 
by a perfect fluid with barotropic index w, then Eq. (|6.65|) becomes 

= - -xyr.^o! (6-66) 
7t[6 + ^{w + IjJ 

which reproduces Eq. ()6.65|) if u; = 1 /3. In the following, on top of the radiation-dominated 
case, the stiff-dominated fluid will be specifically considered. In this case pre-factor will be 
1/6. 

By inspecting the evolution equations of o", various relevant scales can be derived: 



• the moment of oscillations: assuming that sufficiently close to the minimum the axion 
potential is quadratic, coherent oscillations of a will start at a typical scale Hm ~ fn; 

• the moment of dominance: when the Hubble rate is comparable with the axion po- 
tential, i.e. MpHu ~ ma, the axion starts dominating the background; 

• the moment of decay: being gravitationally coupled, a will decay at a typical scale 
Hd ~ m^/Mp] this implies that a must decay prior to big-bang nucleosynthesis in 
order not to jeopardize the production of light elements. 



The last requirement implies that m > 10 TeV taking the BBN temperature of the order of 
0.1 MeV. This bound can be made more restrictive by requiring that the decay takes place 
prior the electroweak epoch (for temperatures of the order of 0.1 TeV). 

When i?ni ~ the oscillations of a commence. For instance, in the case of a quadratic 
potential the effective equation of state of the oscillating field will be the one of dust matter 
and a decreases as a~^/^. This occurrence, however, depends upon the specific properties 
of the potential. For instance in the case of a quartic potential, (p^) ~ ^{pa) 256 , 25S1. 

In the following the main features of the conversion between isocurvature and adiabatic 
modes will be investigated in the two conceptually separate but complementary cases of 
conventional inflationary models and pre-big bang models. 
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6.3.1 Pre-big bang models 



During the pre-big bang evolution the Kalb-Ramond axion is amphfied with quasi-flat 
spectrum jll5| . i.e. 



where the spectral index ria depends on the dynamics of the internal dimensions j2561l257| . 
On the contrary, the spectrum of curvature fluctuations is characterized by a rather steep 
spectral index. More precisely, the considerations introduced earlier in this section allow to 
compute the scalar and tensor spectral indices [5B] : 

ns = 4, nx = 3. (6.68) 

In the context of the ekpyrotic models this estimate leads to 

Us = 3, nx = 2, (6.69) 

assuming continuity in TZ |176j . 

At the beginning of the post-big bang evolution the background is characterized by 
a "maximal" curvature scale Hi, whose finite value regularizes the big bang singularity 
of the standard cosmological scenario, and provides a natural cutoff for the spectrum of 
quantum fluctuations amplified by the phase of pre-big bang inflation. In string cosmology 
models such an initial curvature scale is at most of the order of the string mass scale, i.e. 
Fi<Ms ~ 10^^ GeV. 

The evolution of a is different depending upon the initial value of a in Planck units, i.e. 
Gi/M-p. If o"i > Mp the axion dominates earlier since the condition of dominance is saturated 
faster. If cji < M-p the axion dominates later, i.e. its starts oscillating coherently at Hm and 
it will dominate at a redshift {am/cLa) ~ (cii/Mp)^. This observation also implies that if 
cji > 1, the axion dominance precedes the oscillation epoch. For a more specific discussion 
we refer to i25l USHl HSIl I2SH] • 

Consider then, as an illustration, the case cJi < Mp. Initially GH^^^Mp » m?ai since a 
is not dominant. Later on a oscillates at Hm ~ m and, around this scale the approximation 
of Eq. ()6.66|) breaks down and (j(r) ~ cJi(am/o)^/^. At a typical scale H^j ~ {ai/Mp)^m/36, 
the field starts dominating the background inducing, possibly, a short period of infiationary 
expansion. The axion can decay either before being dominant or after. The decay of 
the curvaton has been discussed both analytically and numerically by different authors 
[2571 126(71 1261j . In the following, for sake of simplicity the approximation of sudden decay 
will be adopted. If the decay follows the dominance, i.e. ffd < H^^, then 

(^) » ■ 

The opposite inequality holds if H^^ > H^j. 

The evolution equations for the fluctuations are obtained by considering the simultane- 
ous presence of the the fluctuations of a and of the fluctuations of the radiation background. 
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The Hamiltonian constraint becomes then, in the longitudinal gauge, 
while the i = j component of the perturbed Einstein equations gives 



AttG 



-a p, 



■6r + iTTGa^dpa, 



where 



Spa 



1 



-tpa + cr X + -Q^a X 



(6.71) 

(6.72) 

(6.73) 
(6.74) 



In all the perturbation equations it has been assumed that anisotropic stresses are absent. 
Notice that we denoted with x the fluctuation of a. Finally the evolution of x ^^^^ of the 
fluid variables is given by 



X" + 2nx' - V^x + ^TT^^X - 4ct'$' + 2^a2$ = 0, 



da 



Si 



(6.75) 
(6.76) 
(6.77) 



We are interested in the solution of the system of Eqs. 1)6. 71() - ()6. 72(1 and H6.75() - (|6.77() with 
non-adiabatic initial conditions, i.e. ip{k,Ti) = 5r{k,Ti) = but x(^)''"i) = Xii^)- 

Since for large-scale inhomogeneities 9i- <^ 1, Eq. (|6.76|) implies 



6y ~ Alp, 



(6.78) 



because initially ip{Ti) 0. Using Eq. 1)6. 73() and (|6.78|) . the Hamiltonian constraint ()6.71|) 
can be written, for — > as 



din a 37^2 



a-'ia'x' + V^a'x]- 



(6.79) 



In order to derive Eq. ()6.79|1 we used the fact that (po- + Pa) ^ Pr- Similarly to Eqs. 1)6. 65|) 
and (|6.66|) . Eq. (|6.75|) implies, with isocurvature initial conditions, that 



X 



1 



Hence, from Eq. (|6.79jl . direct integration leads to 

1 a2 



42Mp ^2 



V,aXi + 0{Vt] 



14 Pr 



VaXi + 0{Vl). 



(6.80) 



(6.81) 
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Inserting Eq. (I^^FTj) into Eqs. and l|H?77)) 



The time evolution of TZ in the radiation-dominated, slow-roll regime can finally be 
determined from the definition of TZ in terms of ip 

n{T)c,^^^Xi + 0{Vl), (6.83) 

so that ip and TZ obey the approximate relation iP{t) ~ — (2/7)7^(t) + 0{V^). 

The same result can be obtained directly by integrating the evolution equation for TZ 

i.e. 

dTZ (5pnad „ .^ 

— — ~ , (6.84) 

a ma P + P 

having approximately determined the form of (5pnad- The exact form of (5pnad of the 
sound speed is, in our system, 

<5pnad = /'rQ - + t^ic^ - l)(p, + p^) + ^(1 - c^) - ^^^(l + Cs'), (6.85) 

Apr + 9{p^ + Pa) + 6a'/nV^- 



' 3 



Ap^ + 3(p^ + p^) 



(6.^ 



When the curvaton is subdominant Cg ~ 1/3 and the leading term in Eq. 1)4.104^ is the last 
one. Hence, to leading order, 5pnad — — 41/o-Xi/3; inserting the approximate form of 5pnad 
into Eq. 1)6. 84() we obtain again, by direct integration Eq. 1)6. 83|) since 

-71 - —X- (6.87) 

ama p^ 

After the initial phase, when the oscillation start, we have to assume a specific form of the 
potential which we take to be quadratic. For H < m, x ^ Xi(om/a)'^/^ and a ~ ai{am/a)^^'^ 
and Eq. ()6.87|) then leads to 

where we used that, for H < m, p^ ^ m^Mp(am/a)^. The second equality follows from the 
definition of r(a), i.e. the time-dependent ratio of the axion energy density over the energy 
density of the dominant component of the background. Recalling now that Xi{k) has a quasi- 
flat power spectrum, Eq. (|6.88)) expresses the conversion of the initial isocurvature mode 
into the wanted adiabatic mode. The expression of TZ{a) will be frozen at the moment of 
sudden decay taking place at the scale H^. In the context of pre-big bang models the power 
spectrum of the axionic fluctuations depends on various parameters connected with the 
specific model of pre-big bang evolution. By comparing the generated curvature fluctuations 
with the value experimentally measured various constraints on the parameter of the model 
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can be obtained j2561 1257] . These constraints seem to suggest that the scale Hi should be 
of the order of 10~^Ms, where Mg ~ O.OlMp is the string mass scale. Notice, furthermore, 
that we just treated the case where the axion decays before becoming dominant. However, 
in the context of pre-big bang models the opposite case is probably the most relevant. If the 
axion decays after becoming dominant numerical calculations have to back the analytical 
approach (especially in the case Ui > Mp). 



6.3.2 Conventional inflationary models 



It is conceivable that the curvaton could also play a role in the context of conventional 
inflationary models, provided the inflationary scale is sufficiently low (i.e. -ffinf < lO^^Mp). 
In this case, it is argued j262j . the adiabatic perturbations generated directly by the fluctu- 
ations of the inflaton, will be too small. In |262j a simple model has been discussed where 
the potentials of the curvaton and of the inflaton are both quadratic. During inflation the 
curvaton should be effectively massless, i.e. m <C -f^inf ■ There are some classes of inflation- 
ary models leading naturally to a low curvature scale |263j (see also |264j ). As far as the 
evolution of the curvaton in a radiation dominated epoch is concerned, the results obtained 
in the case of the pre-big bang are also valid provided the scale Hi is now interpreted as 
the curvature scale at the end of inflation. 

To identify a particle physics candidate acting as a curvaton various proposals have been 
made. Different authors [2651 12661 12671 1268j suggest that the right-handed sneutrinos (the 
supersymmetric partners of the righ-handed neutrinos) may act as curvatons. Models have 
been also proposed where the curvaton corresponds to a MSSM flat direction |2691 127U] and 
to aMSSM Higgs [27T] . 

In the following we would like to concentrate on one of the besic questions arising in 
low-scale inflationary models: is it possible to lower the inflationary scale ad libitum ? 

In the standard curvaton scenario the energy density of the curvaton increases with time 
with respect to the energy density of the radiation background. From this aspect of the 
theoretical construction, a number of constraints can be derived; these include an important 
aspect of the inflationary dynamics occurring prior to the curvaton oscillations, namely the 
minimal curvature scale at the end of inflation compatible with the curvaton idea. Suppose, 
for simplicity, that the curvaton field a has a massive potential and that its evolution, after 
the end of inflation, occurs during a radiation dominated stage of expansion. As previously 
discussed the ratio r(a) increases with time during the radiation-dominated oscillations 



where ai is the nearly constant value of the curvaton throughout the later stages of inflation. 
When a decays the ratio r gets frozen to its value at decay, i.e. r{t) ~ ?^(id) = fd for t > id- 
Equation ()6.89|1 then implies 

2 

m = — (6.90) 

rdMp 

The energy density of the background fluid just before decay has to be larger than the 
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energy density of the decay products, i.e. /Or(id) > T^- Since 

'a 



p^itd)^m'M^{^)\ (6.91) 



the mentioned condition implies 

Td V Mp 

which can also be written, using Eq. 1)6. 9U() . as 

3 



m m , ^ 

> 1, (6.92) 



(6.93) 



Equation 1)6. 93() has to be compared with the restrictions coming from the amplitude of the 
adiabatic perturbations, which should be consistent with observations. If a decays before 
becoming dominant the curvature perturbations at the time of decay are 



7^ td - -JT^ 
Pr oa 



rj^. (6.94) 



Recalling that x\ ~ -f^inf/(2'/r) the power spectrum of curvature perturbations 

^1/2 ^ r^H^ ^ ^ ^ ^^_5 

implies, using Eq. 1)6. 93() together with Eq. 1)6. 90() . 



Recalling now that ra < 1 and m < H[^f, the above inequahty implies that, at most, 
^^inf > 10"^^ Mp if Td ~ 1 MeV is selected. This estimates is, in a sense, general since 
the specific relation between Td and i?d is not fixed. The bound (|6.96jl can be even more 
constraining, for certain regions of parameter space, if the condition Td > \/i/d Mp is 
imposed with f^d — m^ /Mp. In this case, Eq. I|6.96jl implies -ffj^j > lO^^m Mp, which 
is more constraining than the previous bound for sufficiently large values of the mass, i.e. 
m > 10~^-ffinf. Thus, in the present context, the inflationary curvature scale is bound to 
be in the interval |M1 HZU 

IQ-^^Mp < Hini < lO^^Mp. (6.97) 



This bound can be relaxed, by a bit, if the evolution of the curvaton takes place, af- 
ter inflation, in a background dominated by the kinetic energy of the inflaton field itself. 
Consider now the case where the inflationary epoch is not immediately followed by ra- 
diation. Different models of this kind may be constructed. For instance, if the inflaton 
field is identified with the quintessence field, a long kinetic phase occurred prior to the 
usual radiation-dominated stage of expansion. The evolution of a massive curvaton field in 
quintessential inflationary models has been recently studied jl38j in the simplest scenario 
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where the curvaton field is decoupled from the quintessence field and it is minimally coupled 
to the metric. In order to be specific, suppose that, the inflaton potential, V{(p) is chosen 
to be a typical power law during inflation and an inverse power during the quintessential 
regime: 

V{ip)=Xiip^ + M^), ip<0, 

yi^) = ^^.^ ^>0, (6.98) 

where A is the inflaton self-coupling and M is the typical scale of quintessential evolution, 
i.e. by appropriately fixing M, the field will be dominant about today. The potential 
of the curvaton may be taken to be, for simplicity, quadratic. In this model the curvaton 
will evolve, right after the end of inflation, in an environnment dominated by the kinetic 
energy of The curvaton starts oscillating at ~ m and becomes dominant at a 
typical curvature scale ~ m(o"i/Mp)^. Due to the different evolution of the background 
geometry, the ratio r{t) will take the form 

r(t)=.m2(^^y(^^^', H<H^. (6.99) 

to be compared with Eq. 1)6. 89() valid in the standard case. For t > td, r{t) gets frozen 
to the value whose relation to <Ti is different from the one obtained previously (see Eq. 
(|6.9U() ) and valid in the case when a relaxes in a radiation dominated environnment. In 
fact, from Eq. (IFrMj) . 



m ~ — ^. (6.100) 



From the requirement 



Pk(td)-"2^M^(^) >^d^ (6.101) 



it can be inferred, using ()6.100|) . that 

Following the analysis reported in jl38j , the amount of produced fluctuations can be com- 
puted. 

The calculation exploits Eqs. (|5.83|) which are valid for a two-field model. Solving Eqs. 
1)5. 83() with the appropriate initial condition, and using that a'^ip'^/H is constant during the 
kinetic phase, Eq. 1)5. 80(1 can be written as |138j 

7^(^) ~ —crq^ ~ ^ ~ rd — . (6.103) 

Pk Pk da \ ai J 

Recalling that ~ ^inf/(2vr), the observed value of the power spectrum, i.e. V^'^ ~ 
5 X 10"^, implies 



(6.104) 
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The same approximations discussed in the standard case can now be appHed to the case of 
quintessential inflation. They imply the absolute lower bound of i^inf > lOOGeV |272l \'2T3\ 

Up to now we always discussed the case of quadratic potentials. However, the potential 
of a should be allowed to have different profiles. There is only one danger. If the potential 
leads to some attr actor solution for the fluctuations, it may happen that the fluctuations 
of the curvaton are erased |2721 1259] . If, for some reason, the curvaton does not decay 
completely, it is possible to have a residual isocurvature fluctuations which will be completely 
correlated with the adiabatic component (see, for instance j274j ). The large-scale anisotropy 
experiments put then interesting constraints on this possibility (see also section 4). 

7 Spectra of relic gravitons 

While the present value of the Hubble rate, expressed in Hertz, is of the order of 10~^^, relic 
gravitons of frequencies between 100 Hz and the GHz can be also produced throughout the 
cosmological evolution. In the present section we will compute the stochastic background of 
relic gravitons in various cosmological models starting with the case of ordinary inflationary 
models and we will then comment on the possible relevance of the obtained signals for 
experiments aimed at the direct detection of backgrounds of relic gravitons. 

7.1 Relic gravitons from power-law inflation 

In order to estimate the relic graviton background in conventional inflationary models in a 
reasonably generic fashion, we shall be assuming that the evolution of the universe consists 
of three distinct epochs: a power-law inflationary phase, a radiation dominated phase and 
a matter dominated phase. In formulae : 



As already discussed, a generic power-law inflationary phase is characterized by a power p. 
In the case /? = 1 we have the case of the expanding branch of de Sitter space. During the 
radiation-dominated epoch the scale factor expands linearly in conformal time while during 
matter it expands quadratically. Notice that the form of the scale factors given in Eqs. 
(|7.H) - (|7.3]) is continuous and differentiable at the transition points, i.e. the scale factors 
and their first derivatives are continuous in — ti and T2. The continuity of the scale factor 
and its derivative prevents the presence of divergences in the pump field, given by a" /a. 
Notice, incidentally that Eqs. (|7.1|) - H7.3() lead to a pump field even simpler than the one 
reported in Fig. [21 in the intermediate phase (between inflation and radiation) disappears. 
Two scales are present in the problem, i.e. ri and T2. Hence, the spectrum of relic gravitons 






(7.3) 
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will be characterized by two branches limited by the comoving wave-numbers fci ~ 1 /ti and 
k2 ~ 1/t2. For k <ti k2, gravitons will be produced both thanks to the transition inflation 
radiation and to the transition radiation-matter. If, on the contrary, k2 < k < ki, the 
modes of the field will only feel the pumping action of the inflation-radiation transition. 
Finally, if /c > ki, the production of gravitons is absent. 

Consider for instance the case when initial conditions on the tensor modes are set by 
minimizing the Hamiltonian of Eq. (|5.15j) . In this case the evolution of the field operators 
in the Heisenberg description can be written, for r < — ti, as 

Mr) = a,Kro)/i(T) + aK^{ro) f* {r) , (7.4) 

T^kir) = a^{To)gi{T) + a^_^{To)g* (t), (7.5) 

where the mode functions during inflation are given by 

y Ik V ^ 



gr{k,T) = -M\ 



'k 



(7.6) 



where v = f3 + 1/2 and, as usual x = kr. 



For r > — Ti the creation and annihilation operators can be expressed as a linear com- 
bination of the ar and at i.e. 

k 

bj; = B+{k)a^{To) + B^{kyal^{ro), 

bl = B+ikTaliro) + i?-a_,,(ro), (7.7) 

and analogously for k — > —k. During radiation the pump fleld changes since, according to 
Eq. (|7.2() . a(.'/or = 0; the mode functions are then plane waves, i.e. 

1 



Mk,T) 

UKr) = -isj^e-'y , (7.8) 

where y = k[f3T + (/? + 1)ti]. The field operators can then be written in terms of the new 
vacuum annihilated by 6^ and i.e. 

^k = h~9r + blf9;. (7.9) 

Inserting Eq. ()7.7|) into Eq. ()7.9|) we can get the expression of the field operators in terms 
of the "new" mode functions and of the "old" creation and annihilation operators. The 
result is simply 

f^^ir) = a^[B+ik)f, + S-(fc)/;] + al^[B+{krf: + S-(A:)*/r], (7.10) 
7rfc(r) = a^[B+{k)g, + B^g*,] + a^_.[i3+(A:)*5; + B*_g,]. (7.11) 
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Since the evokition of the canonical fields must be continuous, Eqs. H7.4|) - H7.5|) and Eqs. 

Tnm - fmn impiy 

5i(-Ti) = B+{k)gA-n) + B_{k)g:{-T,), (7.12) 

which allows to determine the coefficients of the Bogoliubov transformation B±{k). Insert- 
ing Eqs. dZSl) and (f?^ into Eq. ((7?T^ we get 

B±{k) = y e±-i V^|i/«(xi)[l ± ^(^i.- i)] ^iif«,(xi)|. (7.13) 

Notice that because of the unitarity of the evolution |i?_|_(/c)p — \B_{k)\'^ = 1. Since B^{k) 
and B-{k) are complex, the unitary evolution implies that the Bogoliubov transformation 
is completely determined by three real numbers. 

The mean number of produced pairs of gravitons with momentum k2 < k < ki will then 
be given by 

n{k) = {0\N\0) = \B_\^, (7.14) 
where N = \btbf + b^ J) A. Since k < ki 't-T^, we also have that xi < 1. Hence we have 

k 1^ —k 

that the expression of the mean number of produced particles can be written for <C 1, 
with the result that 

921^-5 

n{k) = T{i^f{l-2uf\kTi\-^''-\ k2<k<ki. (7.15) 

TT 

In the case of pure de Sitter space, P = 1, u = 3/2, n{k) = l/(4x^). 
The energy density of the relic gravitons can be written as 

d'^k 

where the factor 2 accounts for the two polarizations. Hence, the logarithmic energy spec- 
trum can be written as 

The logarithmic energy spectrum can be usefully presented as a fraction of the critical 
energy. Dividing Eq. (|7.17|) by the energy density of the radiation background at ri, i.e. 
Pj ~ 3HiMp{ai/a)'^ we obtain 

where we switched from comoving to physical frequencies. 

Equation ()7.18|) shows that f^GW produced in power-law inflation (and in all conventional 
inflationary models) can be (at most) flat. The completely flat spectrum arises in the pure de 
Sitter case (i.e. /3 = 1). In the case of standard power law and slow-roll inflation (3 > 1 and 
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Oqw decreases. Discounting for a possible variation of the relativistic degrees of freedom, 
the maximal frequency of the spectrum, loi ~ 100^ Hi/Mp GHz., while uj2 ~ 

10-^^ Hz. 

Let us finally compute the spectrum in the infrared branch, i.e. for comoving frequencies 
k < k2. During the matter-dominated stage, the evolution of the mode functions is given 

by 

i 



'2k \ z 



9n 




k[ i 1 
2 ~ i ~ ^2 



(7.19) 



where z = [(3{t + T2) + 2(/? + l)ri]/2. Hence, repeating the same steps outlined above, the 
mode functions during radiation can be matched, in T2 to the mode functions during matter 
by means of two (complex) mixing coefficients, i.e. 

fr{T2)=C+Ur2)+C-C{T2), 

ft(T2) = C+5m(r2) + C-9*M. (7.20) 
The mixing coefficients are then determined to be 

_ _g_^ g±2i [a;i + (x 1 +12 )/3] 



c± 



2[xi + (xi + X2)/3]2 



B±{-1 + 2[xi + (xi + X2)/3] [xi±i + (xi + X2)/3]} 

2[xi + (xi + X2)/3]2 ^ ■ ^ 

The mean number of gravitons will then be, for X2 ^ 1 

(2 i /? - 1) {B+ + B^) -AiB^P 

This result implies that for k < k2 the logarithmic energy spectrum decreases faster than 
in the high-frequency branch. 



OGwoc(-^)'"'^(i^)'. (7.23) 

The spectra of relic gravitons arising in conventional inflationary models, as well as in 
other models are reported in Fig. [21 as a function of the present frequency. Notice that in 
the study of stochastic GW backgrounds the logarithmic energy spectrum is given typically 
as a function of the present physical frequency or as a function of the present physical 
wave-number lo. The relation between the two quantities is w = 2tiv. Notice, furthermore, 
that f^GW enters the expression of the signal-to-noise ratio for the correlation between two 
detectors (say two wide-band interferometers). This formula (see below) depends on the 
value of h, i.e. the present indetermination on the Huble rate. So sometimes it is useful to 
discuss directly /i^JIgw- In terms of the physical frequency we have that the logarithmic 
energy spectrum of relic gravitons in the pure de Sitter case can be written as 



^gw{t^:Tq) c^9.^{TQ)yjJ-j y-^j , l^O<V<V^ec (7.24) 
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Figure 5: The logarithmic energy spectrum of rehc gravitons is iUustrated in different models 
of the early Universe as a function of the present frequency, v. 

where fi^ ~ 7.95 x 10~^ is the present fraction of critical energy density attributed to 
radiation and Hi is the curvature scale of the inflationary phase (strictly constant for this 
estimate). 

7.2 Constraints on the spectra of relic gravitons 

If an inflationary phase is followed by a radiation dominated phase preceding the matter 
dominated epoch, the amplitude of the produced gravitons background can be computed and 
the result is illustrated in Fig. I^lwith the flat plateau. The idea that stochastic backgrounds 
of relic gravitons can be produced in the early Universe goes back to the works of Grishchuk 
j275| I276j (see also |277j for a classical review and |278j for a more recent one) . 

The spectrum consists of two branches a soft branch ranging between ~ 10~^^ h 
Hz and fjec — 10~^^ Hz. For v > i^dec we have instead the hard branch consisting of 
high frequency gravitons mainly produced thanks to the transition from the inflationary 
regime to radiation. In the soft branch Ogw(2^)'^o) ~ i^ "^- In the hard branch Q.q,^^{v^tq) 
is constant in frequency (or slightly decreasing in the quasi-de Sitter case) . The soft branch 
was computed for the first time in |2791 128(71 1281j . The hard branch has been computed 
originally in [2H2j (see also [2831 [28^ [2851 [28fi] V 

The large-scale observation of the first multipole moments of the temperature anisotropy 
imply a bound for the relic graviton background. The rationale for this statement is very 
simple since relic gravitons contribute to the integrated Sachs-Wolfe effect as discussed 
in section 2. The gravitational wave contribution to the Sachs- Wolfe integral cannot be 
larger than the (measured) amount of anisotropy directly detected. The soft branch of the 
spectrum is then constrained and the bound reads 287 

f^Gw(z^,To)/i^< 6.9 X 10~", (7.25) 

for 1/ ~ z^o 10~^^Hz. Moreover, the very small size of the fractional timing error in 
the arrivals of the millisecond plusar's pulses imply that also the hard branch is bounded 
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according to 

^^Gw(i^,ro) < 10-8, (7.26) 

for V ~ 10^8 Hz corresponding, roughly, to the inverse of the observation time during which 
the various mihisecond pulsars have been monitored |289j (see also |288j ). 

The two constraints of Eqs. (|7.25|) and (|7.26|) are reported in Fig. [3 at the two relevant 
frequencies. The Sachs- Wolfe and millisecond pulsar constraints are differential since they 
limit, locally, the logarithmic derivative of the gravitons energy density. There exists also an 
integral bound coming from standard BBN analysis j29Uj and constraining the integrated 
graviton energy spectrum: 

/ VLG^{v,T^)d\i^v< 0.2 h^Vi^ ~ 8 X 10"^ (7.27) 

for h ~ 0.72, Vl^ ~ 7.95 x 10^^. In Eq. (|7.27|) t'max corresponds to the (model dependent) 
ultra-violet cut-off of the spectrum and is the frequency corresponding to the horizon 
scale at nucleosynthesis. Notice that the BBN constraint of Eq. (|7.27|) has been derived in 
the context of the simplest BBN model, namely, assuming that no inhomogeneities and/or 
matter anti-matter domains are present at the onset of nucleosynthesis. In the presence 
of matter-antimatter domains for scales comparable with the neutron diffusion scale [2911 
Hm UnS] this bound is relaxed. 

From Fig. |S1 we see that also the global bound of Eq. (|7.27j) is satisfied and the typical 
amplitude of the logarithmic energy spectrum in critical units for frequencies vi ~ 100 
Hz (and larger) cannot exceed 10"^^. This amplitude has to be compared with the LIGO 
sensitivity to a flat f^Gw(j^/5 '^o) which could be at most of the order of h'^QQy/{uj,TQ) = 
5 X 10"^^ after four months of observation with 90% confidence j294j . 

7.3 Relic gravitons from a stiff phase 

Suppose now, as a toy example, that the ordinary inflationary phase is not immediately 
followed by a radiation dominated phase but by a quite long phase expanding slower than 
radiation |84j . This speculation is theoretically plausible since we ignore what was the 
thermodynamical history of the Universe prior to BBN. If the Universe expanded slower 
than radiation the equation of state of the effective sources driving the geometry had to 
be, for some time, stiffer than radiation. This means that the effective speed of sound Cg 
had to lie in the range l/\/3 < < 1. Then the resulting logarithmic energy spectrum, 
for the modes leaving the horizon during the inflationary phase and re-entering during the 
stiff phase, is tilted towards large frequencies with typical (blue) slope given by ;84j 

^In now _ 6c^ - 2 

which is always limited between (matter) and 1 stiff fluid. A situation very similar to the 
one we just described occurs in quintessential inflationary models |112| . In this case the tilt 
is maximal and a more precise calculation shows the appearance of logarithmic corrections 
in the logarithmic energy spectrum which becomes |295| I29(^ , Ogw oc z^ln^ u (see Fig. E}. 
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The maximal frequency fmax(7"o) is of the order of 100 GHz (to be compared with the 100 
MHz of ordinary inflationary models) and it corresponds to the typical frequency of a spike 
in the GW background. In quintessential inflationary models the relic graviton background 
will then have the usual infra-red and flat branches supplemented, at high frequencies (larger 
than the mHz and smaller than the GHz) by a true spike (296j whose peak can be, in terms 
of h"^ f^GW) of the order of 10~^, compatible with the BBN bound and roughly eight orders 
of magnitude larger than the signal provided by ordinary inflationary models. 

An interesting aspect of this class of models is that the maximal signal occurs in a fre- 
quency region between the MHz and the GHz. Microwave cavities can be used as GW detec- 
tors precisely in the mentioned frequency range (see below). This signal is certainly a candi- 
date for this type of devices. The sensitivity of wide-band interferometers to quintessential 
gravitons has been computed in |297j . 

7.4 Relic gravitons from pre-big bang 

In pre-big bang models /i^ ^gw can be as large as 10~''-10^^ for frequencies ranging be- 
tween 1 Hz and 100 GHz |921l93ir286j . The logarithmic energy spectrum can be either blue 
or violet depending upon the given mode of the spectrum. If the mode under consideration 
left the horizon during the dilaton-dominated epoch the typical slope will be violet (i.e. 
^GW ~ up to logarithmic corrections). If the given mode left the horizon during the 
stringy phase the slope can be also blue (i.e. less steep). An example is reported in Fig. 
ini This behaviour is representative of the minimal string cosmological scenarios. However, 
in the non-minimal case the spectra can also be non monotonic. Recently the sensitivity of 
a pair of VIRGO detectors to string cosmological gravitons was specifically analyzed |298j 
with the conclusion that a VIRGO pair, in its upgraded stage, will certainly be able to probe 
wide regions of the parameter space of these models. If we maximize the overlap between 
the two detectors |29S! or if we would reduce (selectively) the pendulum and pendulum's 
internal modes contribution to the thermal noise of the instruments ^299j , the visible region 
(after one year of observation and with SNR equal to one) of the parameter space will get 
even larger. Unfortunately, as in the case of the advanced LIGO detectors, also in the case 
of the advanced VIRGO detector the sensitivity to a flat spectrum will be irrelevant for 
ordinary inflationary models. 

7.5 Detectors of relic gravitons 

GW detectors can be divided in three broad classes: resonant mass detectors, interferome- 
ters and microwave cavities. There are five (cryogenic) resonant mass detectors which are 
now operating: NIOBE 300 (Perth, Austraha), ALLEGRO pUT] (Baton Rouge, Lousiana, 
USA), AURIGA ISnU (Legnaro, Italy), EXPLORER [SOSI (Geneva, Switzerland) and NAU- 
TILUS |304j (Frascati, Italy). They all have cylindrical shape (the are "bars"). They are all 
made in Aluminium (except NIOBE which is made of Niobium) . Their approximate mass 
is of the order of 2200 kg (except NIOBE whose mass is of the order of 1500 kg). Their 
mode frequencies range from 694 Hz (in the case of NIOBE) to the 912 Hz of AURIGA. 
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There are, at the moment, four Michelson-Morley interferometers being built. They 
are GEO [SOni (Hannover, Germany), TAMA (Tokyo, Japan), VIRGO [M! (Cascina, 
Italy), and the two LIGO |3U8j (in Hanford [Washington], and Livingston [Lousiana], USA). 
The arms of the instruments range from the 400 m of TAMA up to the three km of VIRGO 
and to the 4 km of LIGO. The effective optical path of the photons in the interferometers 
is greatly enhanced by the use of Fabry-Perot cavities. 

Microwave cavities have been originally proposed as GW detectors in the GHz-MHz 
region of the spectrum 309]. A first prototype has been built in MIT in 1978 showing that 
this idea could be actually implemented in order to detct small harmonic displacements 
j310j . It is not unreasonable to think that sensitive measurements could be performed 
in the near future. In particular improvements in the quality factors of the cavities (if 
compared with the prototypes of |309j ^ could be foreseen. Two experiments (in Italy [31 Ij 
and in England |312j ) are now trying to achieve this goal with slightly different technologies. 

8 Radiative transfer equations 

The properties of the radiation field at decoupling are affected by the presence of metric 
fluctuations and a preliminary analysis of this problem has been presented in section 3. The 
effect of metric inhomogeneities on the properties of the radiation field will now be ana- 
lyzed using the radiative transfer (or radiative transport) equations that can be integrated 
numerically even if useful analytical results can also be obtained. A classical preliminary 
reference is the textbook of Chandrasekar |313j (see in particular chapter 1 in light of the 
calculation of the collision term of Thompson scattering that is quite relevant for the present 
ends). Another recent reference is |314j . In broad terms the radiative transfer equations 
describe the evolution of the Stokes parameters of the radiation field through some layer 
of matter which could be, for instance, the stellar atmosphere or, in the present case, the 
primeval plasma prior to the recombination of free electrons and protons. 

In the present framework the radiative transfer equations have a further complication 
with respect to the flat space case: the collisionless part of the Boltzmann equation is 
modified by the inhomogeneities of the geometry. These inhomogeneities induce a direct 
coupling of the Boltzmann equation to the perturbed Einstein equations. An interesting 
system of equations naturally emerges: the Einstein-Boltzmann system of equations which 
is, in some approximation, exactly what has been described in section 4 for the discussion 
of the initial conditions of CMB anisotropics. In that case the perturbed Einstein equations 
were coupled to a set of fluid equations for the density contrasts and for the peculiar 
velocities. These are, indeed, the first two terms (i.e. the monopole and the dipole) in the 
Boltzmann hierarchy. Truncated Boltzmann hierarchies are a useful tool for the analysis 
of initial conditions, but their limitations have been already emphasized in connection with 
the description of collisionless particles. 

While the general conventions established in the previous sections will be consistently 
enforced, further conventions related to the specific way the brightness perturbations are 
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defined Denoting by A a brightness perturbation (related generically to one of the 
four Stokes parameters of the radiation field), the expansion of A in terms of Legendre 
polynomials will be written, in this paper, as 

Aik, n, r) = ^(-i)' i2£ + 1) A^fc, r) Pe{k ■ h). (8.1) 

e 

where k is the momentum of the Fourier expansion, k its direction; n is the direction of the 
photon momentum; Pi{k ■ h) are the Legendre polynomials. The same expansion will be 
consistently employed for the momentum averaged phase-space density perturbation (see 
below Eq. 1)8. 19() ). This quantity will be also called, for short, reduced phase-space density 
and it is related to the brightness perturbation by a numerical factor. The conventions 
of |315l 13161 1318j (see also are such that the factor {2£ + 1) is not included in the 

expansion. Furthermore, in j315ll31b1l317j the metric fluctuations are parametrized in terms 
of the Bardeen potential while in |318j the treatment follows the conformally Newtonian 
gauge. Finally, in [1^1 the conventions are the same as the ones of Eq. 1)8.1(1 but the metric 
convention is mostly plus (i.e. — , +, +, +) and the definition of the longitudinal degrees of 
freedom is inverted (i. e. Ref. [1^] calls what we call (j) Sind viceversa). In j32Ul l32Jll322j 
(see also |323|I324] ) the expansion of the brightness perturbation is different with respect to 
Eq. H8.1|) since the authors do not include the factor {—iY in the expansion. In the latter 
case the collision terms are modified by a sign difference in the dipole terms (involving a 
mismatch of (— i)^ with respect to the conventions fixed by Eq. ()8.1|) '). 

8.1 Collisionless Boltzmann equation 

If the space-time would be homogeneous the position variables and the conjugate mo- 
menta Pj could constitute a practical set of pivot variables for the analysis of Boltzmann 
equation in curved backgrounds. However, since, in the present case, the space-time is 
not fully homogeneous, metric perturbations do affect the definition of conjugate momenta. 
Hence, for practical reasons, the approach usually followed is to write the Boltzmann equa- 
tions in terms of the proper moementa, i.e. the momentum measured by an observer at a 
fixed value of the spatial coordinate. 

Consider, for simplicity, the case of massless particles (like photons or massless neutri- 
nos). Their mass-shell condition can be written, in a curved background, as 

gaf^P^'P" = 0, (8.2) 

where ga/3 is now the full metric tensor (i.e. background plus inhomogeneities) . Equation 
(j8.2|) implies, quite trivially, 

gooP'P'' = -gijP'P' ^ S.jpY, (8.3) 

where the second equality is the definition of the physical three momentum pi. Recalling 
that, to first order and in the longitudinal gauge, goo = 0^(1 + 2(1)) and gij = —a'^{l — 2tlj)6ij, 

^^See Eqs. (18.5511 - 18.581 below for the definition of brightness perturbation. 
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then the relation between the conjugate momenta and the physical three-momenta can be 
easily obtained by expanding the obtained expressions for small (j) and ip. The result is 
simply 

PO = ^(1 - </)) = 4(1 - ^'o = apil +</))= qil + cP), 

P' = ^(1 + V) = ^(l + V'), Pi = -api{l-'tP) = -qi{l-i;). (8.4) 

The quantity qi defined in Eq. H8.4|) are nothing but the comoving three-momenta, i.e. 
Pia = qi, while q = pa is the modulus of the comoving three-momentum. Generalization 
of Eq. (|8.4|) is trivial since, in the massive case, the mass-shell condition implies that 
QapP'^P^ = m? and, for instance P'^ = \/ q^ + m?a^{l — (j)). 

In terms of the modulus and direction of the comoving three-momentum |325j . i.e. 

qi = qni, rnv} = riirijS^^ = 1, (8-5) 
the Boltzmann equation can be written as 

Dt dr dr dx^ dq dr drii dr ' 

where a generic collision term, Ccou has been included for future convenience. 

Eq. (|8.6|) can now be perturbed around a configuration of local thermodynamical equi- 
librium by writing 

f{x\q, nj , r) = /o (q) [1 + f^'^ {x\ q, n, , r)] , (8.7) 

where fo{q) is the Bose-Einstein (or Fermi-Dirac in the case of fermionic degrees of freedom) 
distribution. Notice that /o('?) does not depend on n* but only on q. 

Inserting Eq. (|8.7|) into Eq. (|8.6|) the first-order form of the perturbed Boltzmann 
equation can be readily obtained 

/ota)^ + /„b)^r.- + ^gi=C„, (8.8) 

by appreciating that a pair of terms 

9/(1) dq 5/(1) dn^ 



(8.9) 



dq 9r' drii dr ' 

are of higher order (i.e. C'(V'^)) and have been neglected to first-order. 
Dividing by /o Eq. 1)8.8(1 can also be written as 

^ 9/(1) ^ dlnfo dq _ 1 ^ 
dr dq dr fo ^° ' 

Notice that in Eq. (|8.8|) - (|8.1()|) the generalization of known special relativistic expressions 
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has been used. 

To complete the derivation, dq/dr must be written in explicit terms. The geodesic 
equation gives essentially the first time derivative of the conjugate momentum, i.e. 

^ _ _ _pM pa pl3 /n . 2\ 

where s is the affine parameter. As before in this section, F^^ denotes the full Christoffel 
connection (background plus fulctuations). Using the values of the perturbed connections 
in the longitudinal gauge (take Eq. (|A.10() of the appendix and set i? = £" = 0), Eq. 
becomes 

J pi pj pk 

— = -d'cpP^ + 2V'P' - 2nP' - -y^[d'i;5jk - dk^Si - d^^dl]. (8.13) 

Recalling now that q = qin^, the explicit form of dq/dr will be 
dq 



dr 



apt 

■a^{l-iP)+2na^{l-i))P' -a^i)'P^ 



Ui - P'o^djipn^Ui. (8.14) 



dr 

Inserting now Eq. (|8.13|) into Eq. 1)8.14(1 the explicit form of dq/dr becomes'^^ 

^ = q^' - qmd'cp. (8.15) 
dr 

Finally, using Eq. (|8.15|) into Eq. (|8.10|) to eliminate dq/dr the final form of the Boltzmann 
equation for massless particles becomes: 

^ + "'f? + ^1*' - = |C„„. (8^16) 

or ax* omq Jq 

which can be also written, going to Fourier space, as 

^ + + - ikixct>\ = ICcoii, (8.17) 

or d\nq fo 

where we have denoted, according to the standard notation, k as the Fourier mode and 
H = k-h as the projection of the Fourier mode along the direction of the photon momentum 
Clearly, given the axial symmetry of the problem it will be natural to identify the 
direction of k with the z direction in which case = cos 9. 

The result obtained so far can be easily generalized to the case of massive particles 

df^^^ din fo 1 

— h ia{q, m)kfj,f^'^^ + — [i/;' - ia{q, m)kfj,(j)] = — CcoU, (8.18) 

or dlnq fo 



where a{q,m) = q/ \/q^~+~rn?a? and where, now, the appropriate mass dependence has to 
appear in the equilibrium distribution /q. 

^^To derive Eq. (18.151 from Eq. 118. 1411 . the factors P' and P° appearing at the right hand side Eq. 118.131 
have to be replaced with their first-order expression in terms of the comoving-three momentum g' (and q) 
as previously discussed in Eqs. (18.411 . 

^'^Notice that here there may be, in principle, a clash of notations since, in section 5 (see for instance Eq. 
(|5.66|l 'l we denoted with q and fj, the normal modes for the scalar and tensor actions; in the present section 
q and /x denote, on the contrary the comoving three-momentum and the cosinus between the Fourier mode 
and the photon direction. The two sets of variables never appear together and there should not be obvious 
confusion. 
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8.2 Boltzmann hierarchy for massless neutrinos 



The Boltzmann equations derived in Eqs. (|8.17j) and (|<S.18jl are general. In the following, 
two relevant cases will be discussed, namely the case of massless neutrinos and the case of 
photons. In order to proceed further with the case of massless neutrinos let us define the 
reduced phase-space distribution as 



I q^dqfo 

Eq. (|8.17j) becomes, in the absence of collision term, 



dr 



+ ikfiJ^iy = A(ip' — ikficp). 



i.l9) 



.20) 



The factor 4 appearing in Eq. (|8.20|) follows from the explicit expression of the equilibrium 
Fermi-Dirac distribution and observing that integration by parts implies 



3 , 9fo 
J aq 



-4 / q^dqfo. 
a In g Jo 



.21) 



The reduced phase-space distribution of Eq. (|8.19|1 can be expanded in series of Legendre 
polynomials as defined in Eq. (|8.H) 



Mk,h,T) = Y,{-if{2i+l)T,i{k,T)P,{^). 



(8.22) 



Equation (|8.22|) will now be inserted into Eq. (|8.2()jl . The orthonormality relation for 
Legendre polynomials, 

^ ' (8.23) 



together with the well-known recurrence relation 

{£ + l)P,+i(/i) = (2£ + l)fiPe{fi) - £Pe-i{f^), 



.24) 



allows to get a hierarchy of differential equations coupling together the various multipoles. 
After having multiplied each of the terms of Eq. 1)8. 20(1 by /i, integration of the obtained 
quantity will be performed over (varying between —1 and 1); in formulae: 



l'^P,>{fi)J^,dfi = 2{-iY'j^,e, 
ik J fiP£'{ii)J^^dn = 2ik 



nl A g 

4 / ^ilj'Pi/{fj,)dfi = 8ip'5e'o, , -4:i(j) J ^iJ,Pi'{fi)dfi = --ikcpSpi. 



Equation 



follows from the relation 

1 2 
^fiPe{lJ,)Pe{n)dn = 



f + 1 
21' + I 



2f + 1 



(8.25) 
(8.26) 
(8.27) 

(8.28) 
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that can be easily derived using Eqs. (|5?^ and Inserting Eqs. into 

Eq. I|8.2U|) the first example of Boltzmann hierarchy can be derived: 



T',Q = -kT,i + A^', (8.29) 
Ki = ^[-^.0 - + ^k(t>, (8.30) 

Ki = ^^[^-^-.(^-i) - + (8-31) 

Equation (|OT|) holds for ^ > 2. Eqs. and (|On)) are nothing but the evolution 

equations for the density contrast and for the neutrino velocity field. This aspect can be 
easily appreciated by computing, in explicit terms, the components of the energy-momentum 
tensor as a function of the reduced neutrino phase-space density. In general terms, the 
energy-momentum tensor can be written, in the kinetic approach, as 

T; = -j ^ % n^\P^,r). (8.32) 

According to Eq. ()8.4() . establishing the connection between conjugate momenta and co- 
moving three-momenta, the (00) component of Eq. (|8.32|) becomes, for a completely homo- 
geneous distribution, 

P- = ^J d^QQfoiq), (8.33) 

i.e. the homogeneous energy density. Using instead the first-order phase space density, the 
density contrast, the peculiar velocity field and the neutrino anisotropic stress are connected, 
respectively, to the monopole, dipole and quadrupole moments of the reduced phase-space 
distribution: 

^•^ = ^1 dnj^uik, h, r) = T^o, (8.34) 

3i r ^ ^ 3 

0u = -7r dn{k-h)J^y{k,h,T) = -kj^^u (8.35) 
Ibvr J 4 



o"^ = — I d^l 

167r 



{k-hf-l 



J^,{k,h,T) = ^. (8.36) 



Inserting Eqs. ()8.34|) and 1)8.36(1 into Eqs. I|8.29() - ()8.31() . the system following from the 
perturbation of the covariant conservation equations can be partially recovered 

61 = -^9, + 4^1;', (8.37) 

el = —by - k^Gy + k^<^, (8.38) 

cj'y = - ^kj^,3, (8.39) 
15 10 

with the important addition of the quadrupole (appearing in Eq. (|8.37|) ) and of the whole 
Eq. (|8.39() . which couples the quadrupole, the peculiar velocity field, and the octupole .7^3. 
For the adiabatic mode, after neutrino decoupling, J^j^^ = 0. The problem of dealing with 
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neutrinos while setting initial conditions for the evolution of the CMB anisotropies can 
be now fully understood. The fluid approximation implies that the dynamics of neutrinos 
can be initially described, after neutrino decoupling, by the evolution of the monopole and 
dipole of the neutrino phase space distribution. However, in order to have an accurate 
description of the initial conditions one should solve an infinite hierarchy of equations for 
the time derivatives of higher order moments of the photon distribution function. 

Eqs. (|8.29|) - (|8.,'-{ljl hold for massless neutrinos but a similar hierarchy can be derived also 
in the case of the photons or, more classically, in the case of the brightness perturbations 
of the radiation field to be discussed below. The spatial gradients of the longitudinal 
fluctuations of the metric are sources of the equations for the lowest multipoles, i.e. Eqs. 
(|8.29j) and (|8.3fl|) . For £ > 2, each multipole is coupled to the preceding (i.e. {£ — 1)) and 
to the following (i.e. {I + 1)) multipoles. To solve numerically the hierarchy one could 
truncate the system at a certain £max- This is, however, not the best way of dealing with 
the problem since [IH] the effect of the truncation could be an unphysical reflection of power 
down through the lower (i.e. i < £max) multipole moments. This problem can be efficiently 
addressed with the method of line-of-sight integration (to be discussed later in this section) 
that is also a rather effective in the derivation of approximate expressions, for instance, of 
the polarization power spectrum. The method of line-of-sight integration is the one used, 
for instance, in CMBFAST t32fllH^ . 

8.3 Brightness perturbations of the radiation field 

Unlike neutrinos, photons are a collisional species, so the generic collision term appearing 
in Eq. (|8.18|1 has to be introduced. With this warning in mind, all the results derived 
so far can be simply translated to the case of photons (collisionless part of Boltzmann 
equation, relations between the moments of the reduced phase-space and the components 
of the energy- momentum tensor...) provided the Fermi-Dirac equilibrium distribution is 
replaced by the Bose-Einstein distribution. 

Thompson scattering leads to a collision term that depends both on the baryon velocity 
field and on the direction cosinus |313j . The collision term is different for the brightness 
function describing the fluctuations of the total intensity of the radiation field (related to the 
Stokes parameter /) and for the brightness functions describing the degree of polarization 
of the scattered radiation (related to the Stokes parameters U and V). 

The conventions for the Stokes parameters and their well known properties will now 
be summarized: they can be found in standard electrodynamics textbooks (SHI ( see also 
[3261 I33U1 I331j for phenomenological introduction to the problem of CMB polarization 
and |323j for a more theoretical perspective). Consider, for simplicity, a monochromatic 
radiation field decomposed according to its linear polarizations and travelling along the z 
axis: 

E = [Eie, + E2eyy^'''-'''\ (8.40) 

^^Since the electron-ion collisions are sufficiently rapid, it is normally assumed, in analytical estimates of 
CMB effects, that electrons and ions are in kinetic equilibrium at a common temperature T^h 
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The decomposition according to circular polarizations can be written as 

E = [e+E+ + e_^_]e^('=^-'^*\ (8.41) 

where 

e+ = ^{e^ + iey), (8.42) 



1 

72 



^{e^-icy). (8.43) 



Eq. ()8.42p is defined to be, conventionally, a positive helicity, while Eq. 1)8. 43(1 is the negative 
helicity. Recalling that Ei and E2 can be written as 

Ei = Ej^-, E2 = Eye'^y, (8.44) 

the polarization properties of the radiation field can be described in terms of 4 real numbers 
given by the projections of the radiation field over the linear and circular polarization unit 
vectors, i.e. 

(e^-E), (eyE), {e+ ■ E), {e. ■ E). (8.45) 
The four Stokes parameters are then, in the linear polarization basis 

I=\ec,-E\'^ + \eyE\'^ = El + E^, (8.46) 

Q = \e,-E\^-\eyE\^ = El-El (8.47) 

U = 2Re[(e^ • Ey{ey ■ E)] = 2E^Ey cos {5y - 6^), (8.48) 

V = 2Im[(e, • EYiey ■ E)] = 2E,Ey sin {6y - 4). (8.49) 

Stokes parameters are not all invariant under rotations. Consider a two-dimensional 
(clock-wise) rotation of the coordinate system, namely 

e^. = cos ipix + sin ipiy , 

e'y = — simpex + cos ipiy . (8.50) 

Inserting Eq. (|H3n|) into Eqs. it can be easily shown that /' = / and V = V 

where the prime denotes the expression of the Stokes parameter in the rotated coordinate 
system. However, the remaining two parameters mix, i.e. 

Q' = cos 2lpQ + sin 2ipU, 

U' = -sm2ipQ + cos2ipU. (8.51) 
From the last expression it can be easily shown that the polarization degree P is invariant 



P = ^Q^ + m = ^Q''^ + [7'2, (8.52) 
while U/Q = tan2Q transform as U' /Q' = tan2(a — if). 
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Stokes parameters are not independent (i.e. it holds that l'^ = Q'^ + U'^ + V'^ ) , they 
only depend on the difference of the phases (i.e. {5x — 6y)) but not on their sum (see Eqs. 
(|8.46() - (|8.49() ). Hence the polarization tensor of the electromagnetic field can be written in 
matrix notation as 

/I + Q U-^V \_f El ExEye-^^ \ 
^-[u + ^V I-Q ) = [ExEye^^ El j' ^^"^^^ 

where ^ = {5y — 5x)- If the radiation field would be treated in a second quantization 
approach, Eq. (|8.53|) can be promoted to the status of density matrix of the radiation field 

The evolution equations for the brightness functions will now be derived. Consider, 
again Eq. (|8.17j) written, this time, in the case of photons. As in the case of neutrinos we 
can define a reduced phase space distribution .7%^, just changing u with 7 in Eq. (|8.19l) and 
using the Bose-Einstein instead of the Fermi-Dirac equlibrium distribution. The reduced 
photon phase-space density describes the fluctuations of the intensity of the radiation field 
(related to the Stokes parameter /); a second reduced phase-space distribution, be it Q^, 
can be defined for the difference of the two intensities (related to the stokes parameter Q). 
The equations for and can be written as 

dF 

+ ik^F^ — 4:{ip' — ikfjLcl)) = Cj, 

dC 

+ ikfiG^ = Cq, (8.54) 

The collision terms for these two equations are different [3331 1334j and can be obtained 
following the derivation reported in the chapter 1 of Ref. [313| or by following the derivation 
of Bond (with different notations) in the appendix C of Ref. |39| (see from p. 638). Another 
way of deriving the collision terms for the evolution equations of the brightness perturbations 
is by employing the total angular momentum method j332j that will be swiftly discussed in 
connection with CMB polarization. 

Before writing the explicit form of the equations, including the collision terms, it is 
useful to pass directly to the brightness perturbations. For the fluctuations of the total 
intensity of the radiation field the brightness perturbations is simply given by 



/(x\g,nj,r) = /o(^^^^j. (8.55) 
Recalling now that, by definition. 



the perturbed phase-space distribution and the brightness perturbation must satisfy: 

51n/o- 



1 - Ai(a;*,g,nj,r)- 

amq 



.57) 



104 



that also implies 



At 



-/ 



(1) 



91n/o 
d\nq 



-A- 



/ q^dqfo 



dfo 
dlnq 



I q^dqfo 



4Ai, 



(8.58) 



where the second equality follows from integration by parts as in Eq. (|8.21j) . 
The Boltzmann equations for the perturbation of the brightness are then 



A^ + ikfi{Ai + (/)) = ip' + e' 

A'Q + i/c/xAQ = e'|-AQ + ^[l 
A{j + ikfiAu = -e'Au, 



Ai + Aio + fiVb - ^-P2(/^)5'q 



Ay + i/c/xAv 



Av + 2 ^A* ^vi 



^.59) 

^.60) 
^.61) 
^.62) 



where we defined, for notational convenience and for homogeneity with the notations of 
other authors |32Uj 

flu 

(8.63) 



^ ik 



and 



5q = Ai2 + Aqo + Aq2. (8.64) 

In Eqs. H8.60() ~ (|8.61() . P2{^^) = (S/i'^ — l)/2 is the Legendre polynomial of second order, 
which appears in the collision operator of the Boltzmann equation for the photons due to 
the directional nature of Thompson scattering. Eq. (|8.62|) is somehow decoupled from the 
system. So if, initially. Ay = it will also vanish at later times. 

In Eqs. ()8.59|) - H8.6H) the function e' denotes the differential optical depth for Thompson 
scattering 

e = XeUeaT — = — — , (8.65) 

flo z + 1 

having denoted with Xe the ionization fraction and z = ao/a — l the redshift. Defining with 
To the time at which the signal is received, the optical depth will then be 



er, To) 



a(T) 

XeriedT dr. 

ao 



.66) 



There are two important limiting cases. In the optically thin limit e ^ 1 absorption 
along the ray path is negligible so that the emergent radiation is simply the sum of the 
contributions along the ray path. In the opposite case e ^> 1 the plasma is said to be 
optically thick. Notice that the mean free path of CMB photons discussed in Eq. 1)4. 34() is 
nothing but the inverse of the differential optical depth, i.e. ~ 1/^'- 

To close the system the evolution of the baryon velocity field can be rewritten as 



.67) 



v'b + T-Cvb + ik(j)+-[ 3i An + Vb) = 
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having defined with (3 the ratio between the baryon and photon energy densities, i.e. 

3pb 700 

^ = 4^"rTilo:o23j= 

at the decouphng epoch occurring for z ~ 1100, (5 ~ 7/11 for a typical baryonic content of 
h?Q-\y ~ 0.023. Notice that the photon velocity field has been eliminated, in Eq. ()8.67|) with 
the corresponding expression involving the monopole of the brightness function. 

As pointed out in Eq. (|8.52|) . while Q and U change under rotations, the degree of linear 
polarization is invariant. Thus, it is sometimes useful to combine Eqs. ()8.59() and ()8.60() . 
The result of this combination is 

A'p + {ik^i + e')Ap = ^e'(l - ^?)Sp, 

Sp = Ai2 + Apo + Ap2. (8.69) 
With the same notations Eq. ()8.59|1 can be written as 



AJ + {ik^L + e')Ai = ^' - ik^ic^ + e'[Aio + ^xv^, - ]^P2{ii)Sp]. (8.70) 

By adding a 0' and e'cj) both at the left and right hand sides of Eq. (|8.70j) . the equation for 
the temperature fluctuations can also be written as: 

(Ai + (t>)' + {ik^i + e')(Ai + (!>) = (V' + <t}') + e'[(Aio + 0) + M^b - \p2{^^)Sp]. (8.71) 

This form of the equation is relevant in order to flnd formal solutions of the evolution of 
the brightness equation (see below the discussion of the line of sight integrals). 



8.3.1 Visibility function 

An important function appearing naturally in various subsequent expressions is the so- 
called visibility function, /C(r), giving the probability that a CMB photon was last scattered 
between r and t + dr; the definition of /C(t) is 

/C(r) = e'e-^(^-^o), (8.72) 

usually denoted by g{T) in the literature. The function /C(t) is a rather important quantity 
since it is sensitive to the whole ionization history of the Universe. The visibility function 
is strongly peaked around the decoupling time r^ec and can be approximated, for analytical 
purposes, by a Gaussian with variancG of tliG order of few T^iec 

I335j . In Mpc the width of 

the visibility function is about 70. A relevant limit is the so-called sudden decoupling limit 
where the visibility function can be approximated by a Dirac delta function and its integral, 
i.e. the optical depth, can be approximated by a step function; in formulae: 

/C(r) 5{t - Tdec), e-^(^'^o) ^ 9{t - r^,,). (8.73) 

This approximation will be used, below, for different applications and it is justified since the 
free electron density diminishes suddenly at decoupling. In spite of this occurrence there are 
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convincing indications that, at some epoch after decoupling, the Universe was reionized. The 
ionization fraction Xe becomes again of order 1 and the CMB photons can be rescattered. 
The WMAP results seem to imply that the reionization occurred already between redshift 
z ~ 15 and z ~ 20 I 336j . Temperature anisotropics on scales smaller than the angle 
subtended by the horizon at reionization can be affected j337j and polarized fluctuations 
can be generated. The angular scale corresponding to the horizon at reionization can 
be estimated as "dreio ~ ^f^m/^reio in the case of a spatially flat Universe (see |337j and 
references therein). The apparent detection of early reionization detected by WMAP (for 
15 < z < 20) seems to be compatible with a prolonged reionization process. This occurrence 
would be supported by the detection of traces of smoothly distributed neutral hydrogen 
around z 6 via Gunn-Peterson throughs in the spectra of high-redshift quasars j338l I339j . 



8.3.2 Line of sight integrals 

Equations (|8.69|1 and H8.70() - (|8.71() can be formally written as 

M{k,Ty + {ikn + e')M{k,T) =Af{k,T), (8.74) 

where Ai{k, r) are appropriate functions changing from case to case and J\f{k, r) is a source 
term which also depends on the specific equation to be integrated. 

The formal solution of the class of equations parametrized in the form 1)8. 74() can be 
written as 

M{k,To) = e-^^^^^^n e^^'''^Wik,T)dT, (8.75) 

JO 

where the boundary term for r ^ can be dropped since it is unobservable [31511321] . The 
function A(k, r) determines the solution of the homogeneous equations and it is: 

A{k,T) = / {ik/j, + e')dT = ikfiT + / Xen^aT — dr. (8.76) 
Jo Jo «o 

Using the results of Eqs. the solution of Eqs. and (ISTTj) can be formally 

written as 



(Ai + 0)(fc,ro) = / dr e-''=^^"-^("'"o) (</>' + V'O 
Jo 

+ / dT /C(r 
Jo 



(Aio + (t> + flVh) - ^P2{fl)Sp{k,T) 



.77) 



and as 

3 r° 



Ap{k,To) = - IC{T)e-"'^^^l-fi^)Sp{k,T)dT, (8.78) 

4 Jo 

where e(r, tq) is the optical depth already introduced in Eq. I|8.()(ij) and Ar = (tq — r) is 
the (conformal time) increment between the reception of the signal (at tq) and the emission 
(taking place for r ~ Tdoc)- In Eqs. (|8.77j) and (|8.78j) the visibility function /C(r), already 
defined in Eq. (|8.72|) . has been explicitly introduced. 

Equations (|8.77|) and 1)8.78(1 are called for short line of sight integral solutions. There 
are at least two important applications of Eqs. ()8.77|) and (j8.78j) . The first one is numerical 
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and will be only swiftly described. The second one is analytical and will be exploited both 
in the present section and in the following. 

Equation 1)8. 77() and (|8.78|) can be used to get an integral solution of the Boltzmann 
hierarchy. Let us start from Eq. (|8.78() and let us multiply both sides of the equation by 
Pi{^). Then the right and left hand sides can be integrated over /i between —1 and 1. 
Recalling the conventions for the expansion of the brightness perturbations, i.e. Eq. 1)8. 
the left hand side will simply give, up to a numerical factor, the £-th multipole of Ap, i.e. 
Ap£. The right hand side will instead give an explicit expression. Using repeatedly the 
orthonormality and the recurrence relations, i.e. Eqs. (|8.23j) and (|8.24j) . the result of the 
(tedious) calculation can be expressed as 

Ap£(A:,ro) = ^ T lC{T)S^{k,T)[Udt{k^T) +Ue-2ji-2{kAT) +U,+2k+2{kAr)]dT, (8.79) 

4 JO 

where 

_g e+t-l _ _ (^ + 2)(^ + 3) 

^^-^2^-l)(2^ + 3)' ^^-2- (2^ + l)(2£-l)' ^^+2- (2£ + l)(2£ + 3)- ^^'^^^ 

In Eq. (|8.79j) . ji{k/S.T) are the spherical Bessel functions |2n3| 1204] which can be related to 
the ordinary Bessel functions as: 



j,(fcAr) = ^^^J,+i/2(/cAT). (8.81) 

The j^(fcAT) emerge as a consequence of the expansion in a series of Legendre polynomials 
of the plane wave 

g-^Mfc Ar ^ Y,^-if(2i + l)u{k/\T)P,{^,) (8.82) 

I 

appearing in Eq. 1)8.78(1 . 

The same calculation performed in the case of Eq. I|8.78() can now be performed for Eq. 
(|H?77)) with the result 



e' £{£ + 1) 



(V^' + </.')+e'(Aio + ^ 



4 (2£ + 3)(2^-l 
3 
4 



-Spik,T) 



jeikAr) + ivhe' 



+ 1 

-je-i{kAT) - ^——je+i{kAT) 



21 +V ' ' 21 + 1- 



1) ■ n ^ ^ (^ + l)(^ + 2) . ^ , 

' j£_2(A:Ar) + -A— — Zl^^j^^2(fcAr) 



(2^ + 1)(2£- l)-"" ' (2£ + l)(2^ + 3)' 



5p(A:,r)|. (8.83) 



In the tight coupling approximation (which will be one of the subjects of the forthcoming 
analysis) it is possible to determine with reasonable accuracy the evolution of the monopole 
and of the dipole. Thanks to the recurrence relations of the Boltzmann hierarchy, higher 
multipoles can be also estimated. Then, the integral solutions derived above can be used 
in order to estimate the full angular power spectrum. 

Equation ()8.83|) can be used in order to compute (numerically) the An without solving 
the Boltzmann hierarchy for too high L In this approach it is enough to compute accurately 
the source terms appearing in the integrand at a moderately high value value of 
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Bessel functions appearing in Eq. ()8.83|) have a maximum for kAr ~ £. This occurrence 
provides a relation between the wave-number and the multipole so that a multipole i is 
sensitive to a typical inhomogeneity k ~ I /At. More precisely, the relation between a given 
multipole and the comoving wave-number, can be expressed, in the spatially flat 

k ~ i (8.84) 

6000 Mpc ^ ^ 

where h is the indetermination in the Hubble parameter. 

There is another (related) way of obtaining simplified expressions for the integral solu- 
tions of the Boltzmann equations in the form (|8.69|) - (|8.70j) and (|8.7H) . The formal solution 
of Eq. H8.69() can be written in a different form if the term fi'^ is integrated by parts (notice, 
in fact, that the // enters also the exponential). The boundary terms arising as a result of 
the integration by parts can be dropped because they are vanishing in the limit r ^ and 
are irrelevant for r = tq (since only an unobservable monopole is induced). The result the 
integration by parts of the fj? term in Eq. 1)8. 78() can be expressed as 

Ap(fc,ro)= / e-''"'^^Mp{k,T) dT, (8.85) 
Jo 

^Pik,r) = -^[^(Sp + k^Sp) + '^K,'S'^ + 5p/C"], (8.86) 

where, as usual At = (tq — r). The same exercise can be performed in the case of Eq. (|8.70() . 
Before giving the general result, let us just integrate by parts the term —ikfi(j) appearing at 
the right hand side of Eq. 1)8. 70(1 . The result of this manipulation is 

TO 



Ai{k, To) = / e^'^^^^-^f^-^^^'^o) (V^' + 0') dT 







/"TO 

+ / /C(r)e^^'^(^-^°) dT 
Jo 



(8.87) 



Let us now exploit the sudden decay approximation illustrated around Eq. ()8.73|) and 
assume that the (Gaussian) visibility function /C(t) is indeed a Dirac delta function centered 
around rjec (consequently the optical depth e(T, tq) will be a step function). Then Eq. (|8.87jl 
becomes 



Ai(fc,ro) = e*M— o)[^' + + g^M-dec-ro) [^^^ + ^ + 



where the term has been neglected since it is subleading at large scales. Equation (|8.88() 
is exactly (the Fourier space version of) Eq. (|3.20j) already derived with a different chain 
of arguments and we can directly recognize the integrated SW term (first term at the right 
hand side), the ordinary SW effect (proportional to (Aio + </))) and the Doppler term 
receiving contribution from the peculiar velocity of the observer and of the emitter. 

^^Recall, in fact, that because of the relation between brightness and perturbed energy-momentum tensor, 
i.e. Eqs. iS^Ill and (15351 . 4Aio = 5-y. 
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If all the /i dependent terms appearing in Eq. ()8.7U|1 are integrated by parts the result 
will be 



Ai(fc,ro) 

Ni{k,T) = {1C{t) 
3 



(8.90) 



8.3.3 Angular power spectrum and observables 



Equations (|8.88j) together with the results summarized in section 4 for the initial conditions 
of the metric fluctuations after equality allow the estimate of the angular power spectra in 
the case of adiabatic and isocurvature initial conditions. The brightness perturbation can 
be Fourier transformed ^'^ 



Ai(f,n,r) 



1 



(2^)3/2 



d^ke'''-^Ai{k,h,T) 



.91) 



Assuming that the receiver is located at some conformal time r (eventually coinciding with 
To) and at X = the previous formula can be also expanded in spherical harmonics, i.e. the 
following chain of equalities holds 



1 



em{n) 



=0 m=-e 



(27r)3/2 



(fkAi{k,h,T), 



.92) 



where aim are the coefficients to be determined. Now, as previously done in the present 
section, Ai{k,h,ri) can be represented as series of Legendre polynomials 



A{k, h, t) = J2{-iY{2£ + l)Auik, T)Pt 

1=0 



n] 



(8.93) 



An important property of the is expressed by the theorem of addition of spherical 
harmonics which stipulates that 



PAk ■ h) 



2e + 



m=—i 



Inserting now Eq. ()8.94|) into Eq. ()8.93|) and plugging the obtained result into the second 
equality of Eq. ()8.92|) , the coefficients are determined as 



(27r)3/2 



iY / d'kY;^{k)Au{k,T). 



(8.95) 



Notice that as in the definition of the Fourier transforms of all the quantities introduced so far the factor 
(27r)^/2 has been always incuded. Some authors, however, choose, consistently, not to include it. 
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The two-point temperature correlation function on the sky between two directions con- 
ventionally denoted by hi and n2, can be written as C{'d) = (Ai(ni, ro)Ai(n2, tq)) where 
C('i?) does not depend on the azimuthal angle because of isotropy of the background space- 
time and where the angle brackets denote a theoretical ensamble average. Since the back- 
ground space-time is isotropic, the ensamble average of the a£m will only depend upon £, 
not upon m, i.e. 

{aema£tm') = Ce5u'5mm', (8.96) 

where Ci is the angular power spectrum. Thus, using the first equality of Eq. (|8.92j) and 
exploiting Eq. 1)8. 94() . the relation (|8.96() implies 

C{^) = (Ai(n2,ro)Ai(ni,ro)) = ^(2^ + l)QP^(ni • ns). (8.97) 

The Ce spectrum will now be derived for few interesting examples. Consider, for in- 
stance, the adiabatic mode. In this case Eq. (|8.88() (or, Eq. H3.2U() 1 has vanishing integrated 
contribution and vanishing Doppler contribution at large scales (as discussed in section 4). 
Using then the result of Eq. H4.134|l . the adiabatic contribution to the temperature fluctu- 
ations can be written as 

Af{k,To)=e-"'^^-o[Aio + ^]r,^^^e-'^''^^^^i^^{k), (8.98) 

noticing that, in the argument of the plane wave rjec can be dropped since rjec *S ''"o- The 



plane wave appearing in Eq. (|8.98|1 can now be expanded in series of Legendre polynomials 
and, as a result, 

Ai1(^,ro) = ^^<"(^), (8.99) 
where je{kTo) are defined as in Eq. 1)8. 81(1 with their appropriate argument. 

Assuming now that 'ip^{k) are the Fourier components of a Gaussian and isotropic 
random field (as, for instance, implied by some classes of inflationary models) then 

«(fcX(fc')> = -^rf{k)6(^Hk - k'), Vf{k) = ^\r^{k)\\ (8.100) 

where V^^k) is the power spectrum of the longitudinal fluctuations of the metric after 
equality. Then, Eq. (|8.99() can be inserted into Eq. 1)8. 95() : the obtained result has to be 
plugged into Eq. (|8.96|) and from Eq. (|8.100|) (together with the orthogonality of spherical 
harmonics) In this case 

Cr^ = f /f f<(^)^-^(^-o)^- (8-101) 

To perform the integral it is customarily assumed that the power spectrum of adiabatic 
fluctuations has a power-law dependence characterized by a single spectral index n 

n'(^) = ^IV'fcl' = ^ad(^j . (8.102) 



Recall that there are models where the Universe gets dominated by the cosmological constant when, 
approximately, a/ao > {Q.m/OA)^^''^ ■ This generates a time-evolution in the longitudinal fluctuations of the 
geometry: it is the integrated SW effect which will not be specifically addressed here but which has been 
studied in connection with various problems 1771178117^ previously addressed in this paper. 



Ill 



Notice that kp is a typical pivot scale which is conventional since the whole dependence on 
the parameters of the model is encoded in ^-nd n. For instance, the WMAP collaboration 
[3361 l34Uj , chooses to normalize A at 

A;p = A;i = 0.05 Mpc~\ (8.103) 

while the scalar-tensor ratio (defined in section 6) is evaluated at a scale 

ko = 0.002 Mpc"^ = 6.481 x 10"^^ cm~^ = 1.943 x lO"^'^ Hz, (8.104) 

recalling that 1 Mpc = 3.085 x 10^^ cm. 

Inserting Eq. (|8.102j) into Eq. (|8.101|) and recalling the explicit form of the spherical 
Bessel functions in terms of ordinary Bessel functions 



Cf'^ = ^(to A;p)i-" ^ad dyy^-^Jly,{y), (8.105) 

where y = /ctq. The integral appearing in Eq. (|8.105|) can be performed for —3 < n < 3 
with the result 



oo 



dyy^-'jly,iy) = )^ ( )^ (8.106) 



To get the standard form of the Ce use now the duplication formula for the T function, 
namely in our case 

2 / 25/2-"r(^ 

Insert now Eq. H8.107() into Eq. H8.1U6|) : inserting then Eq. 118.1061) into Eq. H8.105|) we do 

get 



2 

36 

n 1 



r(3-n)r(^ + f 



Zin,e) = {Tokp)'~^2"— ^ (8.108) 

where the function Z(n,i) has been introduced for future convenience. Notice, as a remark, 
that for the approximations made in the evaluation of the SW effects, Eq. (|8.1U8|) holds at 
large angular scales, i.e. £ < 30. 

The C^'**^^ have been given in the case of the spectrum of ^jJ. There is a specific relation 
between the spectrum of ip and the spectrum of curvature perturbations which has been 
derived in section 4 (see discussion after Eq. (|4.125() ') which implies, quite trivially, = 
{25/9)Vf. Finally, the spect rum of the longitudinal fluctuations of the geometry may also 
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be related to the spectrum of the same quantity but computed before equality: this entails 
the (9/10) factor discussed in Eq. (|4.126|) . 

The same calculation performed in the case of the adiabatic mode can be repeated, 
with minor (but relevant) modifications for the CDM-radiation non-adiabatic mode. In the 
specific case of this non-adiabatic mode, Eq. (|<S.l()9j) is modified as 

A;7'')(^,ro) = 2je{kTo)r^^{k), (8.109) 



as it follows directly from Eq. (|8.88|) in the case of non-adiabatic initial conditions after 
equality (see also Eq. (|4.145() ). Performing the same computation Eq. ()8.108|) becomes 

Cf"''^=^2^nad2(nnad,^), (8.110) 

where the power spectrum of non-adiabatic fluctuations has been defined as 



VT'' = A„.J'^ ] . (8.1111 



> =-nad|^^^^ 

Again, following the considerations reported in Eqs. 1)4. 143(1 and ()4.145|) . the spectrum of 
non-adiabatic fluctuations can be directly expressed in terms of the fluctuations of S, i.e. 
the fluctuations of the specific entropy (see Eqs. (|4.1() and (|4. 145(1 ). with the result that 
= {l/25)Vg'^'^ . Of course the major difference between adiabatic and non-adiabatic 
fluctuations will be much more dramatic at smaller angular scales (i.e. say between i ~ 200 
and i ~ 350) where the patterns of acoustic oscillations have a crucial phase difference (this 
aspect will be discussed in the context of the tight coupling expansion) . 

There could be physical situations where adiabatic and non-adiabatic modes are si- 
multaneously present with some degree of correlation. In this case the derivations given 
above change qualitatively, but not crucially. The contribution to the SW effect will then 
be the sum of the adiabatic and non adiabatic contributions (weighted by the appropriate 
coefficients) i.e. 

1 



g.<^(fe) + 2C^'i(fc) 
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While taking expectation values, there will not only be the adiabatic and non-adiabatic 
power spectra, i.e. 'P^{k) and 'P^'^'^{k), but also the power spectrum of the correlation 
between the two modes arising from 



(V'ad(fe)V'nad(^')) = '^V''^' {k)b^''\k - }^ ) , ^^(A;) = V^ad^nad 



cos Oc , 



(8.113) 

where the angle parametrizes the degree of correlation between the adiabatic and non- 
adiabatic mode. The total angular power spectrum will then be given not only by the 
adiabatic and non-adiabatic contributions, but also by their correlation, i.e. 



^(cor) TT^ 



= y V^ad^nadCOSac^(nc,^). (8.114) 
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Consider, finally, the specific case of adiabatic fluctuations with Harrison-Zeldovich, i.e. 
the case n = 1 in eq. (|8.1U8jl . In this case 



+ 1) ^(ad) _ ^ad 



(8.115) 



If the fluctuations were of purely adiabatic nature, then large-scale anisotropy experiments 
(see Fig. UJ imply^^ ^ ~ 9 x 10"^°. Notice that in Fig. ^the quantity Cii{i + l)/(27r) is 
directly plotted: as it follows from the approximate equality 



Ciil{i + l)/(27r) is roughly the power per logarithmic interval of £. 

Up to now the large angular scale anisotropics have been treated. In the following 
the analysis of the smaller angular scales will be introduced in the framework of the tight 
coupling approximation. 

8.4 Tight coupling expansion 

The tight coupling approximation has been already implicitly used in section 4 where it has 
been noticed that, prior to recombination, for comoving scales shorter than the mean free 
path of CMB photons, the baryons and the photons evolve as a single fluid (see Eq. (|4.34|) '). 

If tight coupling is exact, photons and baryons are synchronized so well that the photon 
phase-space distribution is isotropic in the baryon rest frame. In other words since the 
typical time-scale between two collisions is set by ~ 1/e', the scattering rate is rapid 
enough to equilibrate the photon-baryon fluid. Since the photon distribution is isotropic, 
the resulting radiation is not polarized. The idea is then to tailor a systematic expansion 
in Tc ~ 1/e' or, more precisely, in kr^ <C 1 and TcTi <C 1. 

Recall the expansion of the brightness perturbations: 



and Aq^ being the l-th. multipole of the brightness function Ai and Aq. 

The idea is now to expand Eqs. (|8.59j) and (|8.6n|) in powers of the small parameter t^- 
Before doing the expansion, it is useful to derive the hierarchy for the brightness functions in 
full analogy with what is discussed in the appendix for the case of the neutrino phase-space 

^■^To understand fully the quantitative features of Fig. Qit should be borne in mind that sometimes the Ci 
are given not in absolute units (as implied in Eq. H8.115|l but they are multiplied by the CMB temperature. 
To facilitate the conversion recall that the CMB temperature is To — 2.725 x 10^ /iK. For instance the 
WMAP collaboration normalizes the power spectrum of the curvature fluctuations at the pivot scale fcp as 
Viz = (25/9) X (SOOtt^/To ) x A where A is not the A defined here but it can be easily related to it. 





Ai(fc, n, r) = 5](-i)'(2^ + l)Ai,(fc, r)P,(//) 





i 
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distribution. To this aim, each side of Eqs. ()8.59|) - H8.6U() and (|8.67|) wih be multipUed by 
the various Legendre polynomials and the integration over fj, will be performed. Noticing 
that, from the orthonormality relation for Legendre polynomials (i. e. Eq. (|8.23|) '). 



and recalling that 



.118) 



Po(/u) = l, A(m) = m, P2(^) = ^(3//2-1), P3(^) = 1(5^3_3^)^ (8^119) 

Eqs. (|8.59|) - (|8.6flj) and (|8.67j) allow the determination of the first three sets of equations 
for the hierarchy of the brightness. More specifically, multiplying Eqs. H8.59() - ()8.6U() and 
(|8.67j) by PoifJ-) and integrating over fi, the following relations can be obtained 



A' + kAQi = -[Aq2 + Ai2 - Aqo], 



^.120) 
^.121) 



v'f, + Hv], = -ikcj) - — (3iAii + Vb)- 



.122) 



If Eqs. (|H3!?|) - (pn)n|) and (P0)7|) are multiplied by Pi(^), both at right and left-hand sides, 
the integration over jj, of the various terms implies, using Eq. (|8.118|) : 



-Ail - \kAi2 + ^Aio 



' + e' 



An + — 

61 



-A^i - ^/cAq2 + ^Aqo = e'Aqi, 
v'b + Tivb = -ik(j) - ^(3iAii + vt). 



The same procedure, using P2{lJ'), leads to 

3 2 
-A'i2 - -kAi3 + -kAn = e' 



^Ai2 - ^(Aqo + Aq2) 



-A'q2 - ^/cAq3 + hAQi = e' 



—A - — 
10 10 



+ A12) 



v'b + Hvh = -ik(j) - ^ ^3iAii + V}}j . 



].123) 
i.l24) 
i.l25) 

i.l26) 
].127) 
i.l28) 



For £ > 3 the hierarchy of the brightness can be determined in general terms by using the 
recurrence relation for the Legendre polynomials reported in Eq. (|8.24|) : 



A[^ + e'Au = 



2£ + 



-[£Aq(,_i) -(^ + 1)Aq(,+i)]. 



^.129) 
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8.4.1 Zeroth order in the tight coupling expansion: acoustic oscillations 



We are now ready to compute the evolution of the various terms to a given order in the 
tight-coupling expansion parameter Tc = After expanding the various moments of 

the brightness function and the velocity field in 

Ak = Ai£ + Tc6u, 
Aq£ = Aq^ + Tc(5q£, 

Vh = v\, + Tc6^^, (8.130) 

the obtained expressions can be inserted into Eqs. (|8.12()j) - (|8.125|) and the evolution of the 
various moments of the brightness function can be found order by order. 

To zeroth order in the tight-coupling approximation, the evolution equation for the 
baryon velocity field, i.e. Eq. (|8.122|) . leads to: 

?Jb = -3iAii, (8.131) 

while Eqs. (|8.121|) and (|8.124|) lead, respectively, to 

Aqo = Ai2 + Aq2, Aqi = 0. (8.132) 

Finally Eqs. (I8.12(ij) and Il8.127j) imply 

9 Ai2 = Aqo + Aq2 , 9 Aq2 = Aqo + Ai2 . (8. 133) 

Taking together the four conditions expressed by Eqs. (|8.132j) and (j8.133j) we have, to 
zeroth order in the tight-coupling approximation: 

Aq^ = 0, £>0, Ak = 0, l>2. (8.134) 

Hence, to zeroth order in the tight coupling, the relevant equations are 

Ufe = -3iAii, (8.135) 
A^o + fcAii = V'. (8.136) 

This means, as anticipated, that to zeroth order in the tight-coupling expansion the CMB 
is not polarized since Aq is vanishing. 

A decoupled evolution equation for the monopole can be derived. Summing up Eq. 
1)8. 123(1 (multiplied by 3i) and Eq. (|8.125|1 (multiplied by /?) we get, to zeroth order in the 
tight coupling expansion: 

f3v'^ - 3iAn + ik(t){(j + 1) - 2iA;Ai2 + ikKio + PHvi, = 0. (8.137) 

Recalling now Eq. (|8.135|) to eliminate Ub from Eq. I|8. 137(1 . the following equation can be 
obtained 

if3 + 1)A;i + nf3 An - ^Aio = 0. (8.138) 
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Finally, the dipole term can be eliminated from Eq. ()8.138() using Eq. H8.136|) . By doing 
so, Eq. (|8.138jl leads to the wanted decoupled equation for the monopole: 



Ql 



ill . 1 1 
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where 



Csb = ^^=^, (8.140) 

is the sound speed which includes the effect of baryons. The term k'^c^AiQ is the photon 
pressure. Defining, from Eq. (|8.140|) . the sound horizon as 

rs(r) = r c,{T')dT', (8.141) 
JO 

the photon pressure cannot be neglected for modes krs{T) > 1. 

At the right hand side of Eq. (|8.139j) several forcing terms appear. The term ^" 
dominates, if present, on super-horizon scales and causes a dilation effect on Aiq. The term 
containing k'^cp leads to the adiabatic growth of the photon-baryon fluctuations and becomes 
important for A;r ~ 1. 

In Eq. (|8.139|1 the damping term arises from the redshifting of the baryon momentum in 
an expanding Universe, while photon pressure provides the restoring force which is weakly 
suppressed by the additional inertia of the baryons. 



8.4.2 Solutions of the evolution of monopole and dipole 

Equation (18.1 39|1 can be solved under different approximations (or even exactly |315j ). The 
first brutal approximation would be to set = /? = 0, implying the the role of the baryons 
in the acoustic oscillations is totally neglected. As a consequence, in this case Cg^b = l/\/3 
which is nothing but the sound speed discussed in Eqs. (|4.13()j) - (|4.133|) for the fluid analysis 
of the adiabatic mode. In the case of the adiabatic mode, neglecting neutrino anisotropic 
stress, = (j) = and ^p' = 0. Hence, the solution for the monopole and the dipole to 
zeroth order in the tight coupling expansion follows by solving Eq. ()8.139|) and by inserting 
the obtained result into Eqs. (|8.135j) and (|8.136|) . i.e. 

— "0 

Aio(/c,r) = — [cos (A;Cs,bT) -3], 

Aii(fc,T) = — ^/cCs,bsin(/i:Cs,bT), (8.142) 

which is exactly the solution discussed in section 4 if we recall Eq. (|8.135|) and the definition 

If = but /3 / 0, then the solution of Eqs. (18.1351) - (I8.13(ill and (18.13911 becomes, in 
the case of the adiabatic mode, 

Aio(fc, r) = ^(/3 + l)[cos (fccs,br) - 3], 

Aii(A;,r) = J ^ sin (kc^^^r) . (8.143) 
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Equation (|8. 143(1 shows that the presence of the baryons increases the ampHtude of the 
monopole by a factor (3. This phenomenon can be verified also in the case of generic time- 
dependent (3. In the case of /3 7^ the shift in the monopole term is {(3 + 1) with respect 
to the case /? = 0. This phenomenon produces a modulation of the height of the acoustic 
peak that depends on the baryon content of the model. 

Consider now the possibility of setting directly initial conditions for the Boltzmann 
hierarchy during the radiation dominated epoch. During the radiation dominated epoch 
and for modes which are outside the horizon, the initial conditions for the monopole and 
the dipole are fixed as 

00 525 + 188i?. + iml ^ 2 2 
= -y - 180(25 + 2i?.) "^'^ " ' 

where 0o is the constant value of (j) during radiation. This result is fully compatible with 
the results obtained in Eqs. (|4.55() ~ (j4.(i3|) . The constant value of ■0, i-e. V'o will be related 
to 00 through i?,^, i.e. the fractional contribution of the neutrinos to the total density (see 
Eqs. ((Oil and dOTll '). 

It is useful to observe that in terms of the quantity Aq = (Ajq — "0)) Eq. ((8.1H9() becomes 



A[)' + fc^c^ Ao = -k" 



V , 2 I 
3 + Cs.bV' 



^.145) 



The initial conditions for Aq are easily obtained from its definition in terms of Ajo and tp. 

The same strategy can be applied to more realistic cases, such as the one where the scale 
factor interpolates between a radiation-dominated phase and a matter-dominated phase, as 
discussed in Eq. ((4.83() . In this case the solution of Eq. ((8.139(1 will be more complicated but 
always analytically tractable. Equation ((8.139(1 can indeed be solved in general terms. The 
general solution of the homogeneous equation is simply given, in the WKB approximation, 
as ^ 

Aio = ^i^^-^y/i COS fcrs + Bs\nkr^\. (8. 146) 

For adiabatic fluctuations, A;^0 contributes primarily to the cosinus. The reason is that , in 
this case, ip is constant until the moment of Jeans scale crossing at which moment it begins 
to decay. Non-adiabatic fluctuations, on the contrary, have vanishing gravitational potential 
at early times and their monopole is dominated by sinusoidal harmonics. Consequently, the 
peaks in the temperature power spectrum will be located, for adiabatic fluctuations, at a 
scale kn such that knrs{T*) = nir. Notice that, according to Eq. ((8.136(1 the dipole, will be 
anticorrelated with the monopole. So if the monopole is cosinusoidal , the dipole will be 
instead sinusoidal. Hence the "zeros" of the cosinus (as opposed to the maxima) will be 
filled by the monopole. The solution of Eq. ((8.139(1 can then be obtained by supplementing 
the general solution of the homogeneous equation ((8.146(1 with a particular solution of the 
inhomogeneous equations that can be found easily with the usual Green's function methods 
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The amplitude of the monopole term shifts as (1 + f))^"^/^ . RecalHng Eq. ()8.68|) . it can 
be argued that the height of the Doppler peak is weakly sensitive to /i^Ob in the ACDM 
model where fib ^ and /?(Tdcc) < 1- 

In £ space the position of the first peak for adiabatic and isocurvature modes is given, 
respectively, by 

rs{Tdcc) 

£W=fn + iV%^' (8.148) 

where, in the spatially flat case 

1 /"^doc dz 

^""''^'^H^oJo Va(l + ^)^ + fim(l + ^)^ + fiA' ^ ' 

and adoc(''"dec — ''"o) is the angular diameter distance. In the case = and for /i^fib = 0.023 
(whose value affects the sound horizon) the position of the first peak will be about 220 
for adiabatic modes and about 330 isocurvature modes. The position of the first peak 
detected in large scale CMB experiments is in good agrrement with Eq. (|8.147() (see the 
introduction). The predicted position of the first adiabatic peak is in good agreement with 
the value inferred from the data reported in Fig. ^ . Notice that when integrating up to 
Zdec ~ 1100 the radiation density (i.e. fir ~ 8 x 10~^) does affect the result since it is 
multiplied by a factor (1 + z)^. Notice, conversely, that the region of small redshifts (i.e. 
z < 10) is sensitive to the presence of a cosmological constant. 

It is difficult to obtain general analytic formulas for the position of the peaks. Degenera- 
cies among the parameters may appear j341| . In j342j a semi-analytical expression for the 
integral giving the angular diameter distance has been derived for various cases of practical 
interest. 

Once the evolution of the lowest multipoles is known, the obtained expressions can be 
used in the integral solutions of the Boltzmann equation (like Eq. (|8.83j) ) and the angular 
power spectrum can be computed analytically. Recently Weinberg in a series of papers 
[3431 13441 1345j computed the temperature fluctuations in terms of a pair of generalized 
form factors related, respectively, to the monopole and the dipole. This set of calculations 
were conducted in the synchronous gauge (see also |3461 13471 1348j for earlier work on this 
subject). Reference |349j also presents analytical estimates for the angular power spectrum 
exhibiting explicit dependence on the cosmological parameters in the case of the concordance 
model. 

The results of the tight coupling expansion hold for krc ^ 1. Thus the present approxi- 
mation scheme breaks down, strictly speaking, for wave-numbers k > r~^. Equation ()8. 139(1 
holds to zeroth-order in the tight coupling expansion, i.e. it can be only applied on scales 
much larger than the photon mean free path. By comparing the rate of the Universe expan- 
sion with the rate of dissipation we can estimate that Tck^ ~ defines approximately the 
scale above which the wave-numbers will experience damping. From these considerations 
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the typical damping scale can be approximated by 

~ 0.3(rJm/i2)"^/2(Ob/i^)"^(a/adoc)^/^ Mpc^. (8.150) 

The effect of diffusion is to damp the photon and baryon oscillations exponentially by the 
time of last scattering on comoving scales smaller than 3 Mpc. For an experimental evidence 
of this effect see ^8| and references therein. 

In order to have some qualitative estimate for the damping scale in the framework of 
the tight coupling approximation, it is necessary to expand the temperature, polarization 
and velocity fluctuations to second order in t^. Since for very small scales the role of gravity 
is not important the longitudinal fluctuations of the metric can be neglected. The result of 
this analysis j321j shows that the monopole behaves approximately as 

where 

kl Jo 6(/3 + 

The factor 16/15 arises when the polarization fluctuations are taken consistently into ac- 
count in the derivation |321j . 
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8.4.3 First order in tight coupling expansion: polarization 

To first order in the tight-coupling limit, the relevant equations can be obtained by keeping 
all terms of order Tc and by using the first-order relations to simplify the expressions. From 
Eq. (18.1241) the condition 5qi = can be derived. From Eqs. (|8.121B and (I8.12eil) - (l8.127j) . 
the following remaining conditions are obtained respectively: 

-Sqo + Si2 + 6q2 = 0, (8.153) 

^5i2 - ^[6qo + 6q2] = IkKn, (8.154) 

^Sq2-^[Sqo + 6u]=0. (8.155) 

Equations ()8.153() - ()8.155|) are a set of algebraic conditions implying that the relations to 
be satisfied are: 

5qo = lSi2, (8.156) 

'^Q2 = \6i2, (8.157) 

6i2 = ^kAn. (8.158) 
15 

Recalling the original form of the expansion of the quadrupole as defined in Eq. ()8.130|) . 
Eq. (|8.158|1 can be also written 

g 

Ai2 = TcSu = — /cTcAii, (8.159) 
15 
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since to zeroth order the quadrupole vanishes and the first non-vanishing effect comes from 
the first-order quadrupole whose value is determined from the zeroth-order monopole. 

Now, from Eqs. (|8.156() and H8.157|) . the quadrupole moment of Aq is proportional to 
the quadrupole of Aj, which is, in turn, proportional to the dipole evaluated to first order 
in Tc- But Aq measures exactly the degree of linear polarization of the radiation field. So, 
to first order in the tight-coupling expansion, the CMB is linearly polarized. Notice that 
the same derivation performed in the case of the equation for Aq can be more correctly 
performed in the case of the evolution equation of Ap with the same result |321j . Using 
the definition of Sp (i.e. Eq. (|09|) ). and recalhng Eqs. (|8.156j) - (|8.158|) . we have that the 
source term of Eq. (|8.78|) can be approximated as 

Sp ~ ^A:rcAii. (8.160) 

Since Tc grows very rapidly during recombination, in order to have quantitative estimates of 
the effect we have to know the evolution of Sp with better accuracy. In order to achieve this 
goal, let us go back to the (exact) system describing the coupled evolution of the various 
multipoles and, in particular, to Eqs. H8.121|) and H8.126|) - ()8.127|) . Taking the definition of 
Sp (or Sq) and performing a first time derivative we have 

S^ = A'j^ + A'p^ + A'po. (8.161) 

Then, from Eqs. 1)8.121(1 and (|8.126() - ()8.127|) . the time derivatives of the two quadrupoles 
and of the monopole can be expressed in terms of the monopoles, quadrupoles and octupoles. 
Simplifying the obtained expression we get the following evolution equation for Sp, i.e. [321j 



3 

Sp + Y^e'Sp = k 



^Aii - ^ ( Api + Ap3 + Ai3 



i.l62) 



This equation can be solved by evaluating the right hand side to zeroth-order in the tight 
coupling expansion, i.e. 

Sp{t) = -k TdxAne-^^^^'^). (8.163) 



5 Jo 

This equation, giving the evolution of S'p, can be inserted back into Eq. ()8.78|) in order to 
obtain Ap. The result of this procedure is 321] 

Ap ~ (1 - /u2)e*'=^("'d--^«)p(A:), 

P(/C) ~ (0.51) k CTdecAii(rdec), (8.164) 

where Udec is the width of the visibility function. This result allows to estimate with rea- 
sonable accuracy the angular power spectrum of the cross-correlation between temperature 
and polarization (see below), for instance, in the case of the adiabatic mode. 

While the derivation of the polarization dependence of Thompson scattering has been 
conducted within the framework of the tight coupling approximation, it is useful to recall 
here that these properties follow dirrectly from the polarization dependence of Thompson 
scattering whose differential cross-section can be written as 

^ = r2|eW . ^ . ,{a')|2. ^8.165) 

df] °' ' 87r ' ' ^ ^ 
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where e^"^ is the incident polarization and e^"'-* is the scattered polarization; tq is the 
classical radius of the electron and ctt is the total cross-section both defined in Eq. (|4.32|) . 

Suppose that the incident radiation is not polarized, i.e. U = V = Q = 0; then we can 
write 
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Defining the incoming and outgoing polarization vectors as 

= (1, 0, 0), ey = iO, 1, 0), A; = (0, 0, 1). 
e'^ = {—smip, — COS93, 0), = (cos t9 cos — cosi?sin(^, — sin??), (8. 167) 

the explicit form of the scattered amplitudes will be 

3(Tt 



I' 



^a; '^x\ ^ ' ^x\ -'"I 



' |2 

y ' ^x\ -^y 



3aT 



X, 



3(Tt 
"8^ 



If .f'pT +|f .f'pT 

\tx ty\ J.X ^ \ty ty\ -Ly 



167r 

— -^Tcos^ -d. 
Wit 



.168) 



Recalling the definition of Stokes parameters: 



I'=I'^+I' = —arAl + cos^ ??), 
^ 167r 



3 

16^' 



^.169) 



Even if [/' = the obtained Q and U must be rotated to a common coordinate system: 

Q' = cos 299Q, U' = - sin 2LpQ. (8.170) 
So the final expressions for the Stokes parameters of the scattered radiation are: 

= t|-^tX(1 + cos2i?), 



Q' = (TT^sin^ i?cos2(/7, 

167r 

U' = cTT^^sin'^ t?sin2(i3. 

IGvr 

We can now expand the incident intensity in spherical harmonics: 

£m 

So, for instance, Q' will be 

3 f 

Q' = TT^f^T / ^Yem{'&, 'P) aim sin^ '& COS 2(fdQ. 
Ibvr J ^ 

Ira 



5.172) 



(8.173) 



By inserting the explicit form of the spherical harmonics into Eq. H8.173|) it can be easily 
shown that Q' 7^ provided the term 022 7^ in the expansion of Eq. (|8.172j) . 
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Figure 6: Temperature-polarization correlation measured by the WMAP collaboration 
(adapted from Ref. j336j . 



Recall now that under clockwise rotations the Stokes parameters Q and U transform as 
in Eq. H8.51() . As a consequence 

{Q±iUy = e'^^''^{Q±iU) (8.174) 
where is the rotation angle. From this observation it follows that the combinations 

(Aq ±iAu)(n) = ^a±2/m ±2yem{n) (8.175) 

im 

can be expanded in terms of the spin-2 spherical harmonics, i.e. ±23^™(^) |35Ul 1351] . 
The expansion coefficients are 

a±2/m = Jdn ±23^;„(AQ±iAu)(n). (8.176) 

In Ref. |35U1 1351j . the authors introduce the following linear combinations of a±2/m to 
circumvent the impasse that the Stokes parameter are not invariant under rotations; 

1 i 

dE/m = — 2 + a-2/ni)aB/m = ^ — a-2/m)- (8.177) 

These newly defined variables are expanded in terms of ordinary spherical harmonics, 
^i^) = ^dEimyimin), B{n) = ^aB/mYimin)- (8.178) 

Im Im 

The spin-zero spherical harmonics, Yim{h), is free from the ambiguity with the rotation of 
the coordinate system, and therefore E and B are rotationally invariant quantities. The 
i?-mode has (—1)^ parity and the B mode (— l)(^+i) parity in analogy with electric and 
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magnetic fields. Scalar perturbations generate only the E mode |353| . While the scalar 
fluctuations only generate an E mode, tensor fluctuations also generate a B mode. The 
Boltzmann equation for the tensor modes can be easily derived following Refs. [3541 1355j 
(see also |H^rM] V 

Consider now, specifically, the adiabatic mode. While the temperature fluctuation Ai 
oscillates like cos (fecs,bTdec)) the polarization is proportional to the dipole and oscillates 
like the sinus of the same argument. The correlation function of the temperature and 
polarization, i.e. (AjAp) will then be proportional to sin (/ccs^b'^'dcc) cos (fecs^bTdec)- An 
analytical prediction for this quantity can be inferred from Eq. (|8.164|) (see |321j ). The 
spectrum of the cross correlation must then have a peak for fccg^b'^dec ~ 37r/4, corresponding 
to £ ~ 150. This is the result suggested by Fig. Eland taken from Ref. ^333 reporting the 
measurement of the WMAP collaboration. In Fig. El the temperature-polarization angular 
power spectrum is reported for adiabatic models (solid line) and for isocurvature models 
(dashed line). 

We should mention here that a rather effective method in order to treat on equal footing 
the scalar, vector and tensor radiative transfer equations is the total angular momentum 
method j332[l357|l358j . Within this approach, the collision terms couple only the quadrupole 
moments of the distributions and each moment corresponds directly to observable patterns 
in the microwave sky. In this language the analysis of the polarization of the radiation field 
becomes somehow more transparent. 
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A Metric fluctuations 



In this appendix, some technical aspects concerning the scalar, vector and tensor fluctua- 
tions of the geometry will be summarized with the purpose of making the present review 
self-contained. The fluctuations of the geometry and of the sources are parametrized as in 
section 2. The background geometry will be assumed to be spatially flat. To first order in 
the amplitude of the metric fluctuations, the fluctuations of the Christoffel connections can 
be computed 

^'^afs = ^9^''i-duSgai3 + dpSg^a + daSg^u) + ^^g^^i-dugap + dpg^a + dagf3„), (A.l) 

where the overline denotes that the corresponding quantity is to be evaluated on the back- 
ground geometry and the usual convention of summation over repeated indices is adopted. 
The first-order fluctuation of the Ricci scalar can also be computed as 

SRf,. = 5„<5r-, - d,6T% + 5r^, T% + t;, 5rf^ - ^r^^r;, - ^<5r^,. (A.2) 

The fluctuations of the Ricci tensors with one controvariant index and one covariant index, 
as well as the fluctuations of the Ricci scalar, are then 

5R; = 6g'"'R^^+r''6Ra^,, (A.3) 
5R = 6g^PR^p+T^5R^p. (A.4) 

Finally, the fluctuations of the components of the Einstein tensor can be easily deduced 
once Eqs. ()A.3|) - (|A.4|) are known in explicit terms: 

Formally, the perturbation of the covariant conservation equation is: 

d^ST^'" + T'^^JT'"' + JF^l^T"" + r^pST""^ + ST^^pT"^ = 0. (A.6) 

To obtain explicit formulae for scalar, vector and tensor modes it is necessary to have the 
values of the Christoffel connections of the background: 

The components of the Ricci tensor and the Ricci scalar of the background are then: 

^ 3 
Rqo = — SW', Rq = — kTC , 

Rij = in' + 2n'^)6ij, Ri = -^{n' + 2n^)5{, 

R = -^{n^ + n'). (A.8) 



125 



A.l The scalar modes 



The scalar fluctuations of the geometry with covariant and controvariant indices can be 
written to first-order , as 

4500 = 2a^0, 
6sgij = 2a^{ii6ij 
Ssgoi = -o^diB, 

Inserting Eqs. ()A.9|) into Eq. (|A.H) the various first-order components of the (scalar) 
Christoffel fluctuations can be easily computed: 

5sro. = -[ip' + 2H(0 + ^pMj + didj[{E' + 2nE) - b], 

5,Tio = d\<p + B' + nB), 

5, rf^. = {d^i^Sij - dii^S^ - dji^Sf) + d.djd^E - m^B5,j, 

6, Ti^ = -^'5{ + d,d^E'. (A.IO) 

The first-order fluctuations of the Ricci tensors are 

S.Roo = V^[^+iB-E'y + niB-E')] + 3[t{;" + ni^' + ij')], 

dsRij = -5ij{[^" + 2{n' + 2n^){^ + ^) + n{(i)' + 5i^')-v^'iP]+nv\B- E')}, 
+ didj [{E' - B)' + 2{n' + 2n^)E + 2niE' -b) + {^- 0)] (a.ii) 

while the first-order fluctuation of the Ricci scalar are determined to be, 
5^R = A{3^"+6(7^'+7^2^^^3-^(^^/^3^/^^y2^^^_2^,)^(^_^/y+3-^(^_^/)]}^ (A.12) 

Using Eqs. ()A.11|) into Eqs. ()A.8|) and (|A.4|1 . the fluctuations of the components of the 
Ricci tensor with mixed indices will then be: 

s.R'l = ^{v^[(i) + {B-E')' + n{B-E') + 'i[i:" + n{<i)' + ^') + 2n'4>\}, 
6,4 = ^[ij" + 2{n' + 2n'^)<i) + n{<t>' + bi)')-v'^i) + nv'^{B - E')]5{ 

- \did^[{E' - B)' + 2n{E' -B) + {^- 0)], 

6sR}) = ^d'[-{^' + n<p) + {n' -n^)B]. (a.is) 



Os9 



00 



didjE) , = - ^ {iP5'^ - d'd^E) , 



Pi 



-^d'B. 



(A.9) 
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Finally the fluctuations of the components of the Einstein tensor with mixed indices are 
computed to be 

= ^{v^i'-nv\B-E')-m{ij' + n(i))}, (a.m) 

5,gi = J-{[-2<- 2(7^2 + 2W')0-2W(/''-4W^'] 

- V\(l)-iP) + {B-E'y + 2n{B-E')]}5i 

- ^did^[{E' -B)' + 2n{E' -B) + {i^-4>)], (A.15) 
,5sg° = 5Rl (A. 16) 



A. 2 The vector modes 



The vector fluctuations of the geometry can be written as 
o^goi = -a Qi, d^g = 



S.gij = a\diWj + djWi), 5^g'^ = -^{d'W^ + dm'), (A.17) 



and the related fluctuations of the Chrristoffel connections can be deduced, after some 
algebra, by inserting Eqs. HA.17|) into Eq. (|A.1|) . The result is: 

5, E°. = -^{d^Qj + djQi) - H{diWj + djWi) - \{d,w] + djW;), 

6. ri^, = Q'' + nQ\ 

,5,E^, = -HQ'^d,, + Id'^idiWj + djWi) - Idjid'^w^ + diW'^) - IdiidjW'^ + d'^Wj), 

■' z z z 

-^vFio = \m' - d'Qi) - lid'Wl + diW^'), (A.18) 
where, the conditions of Eq. (|2.15() on Qi and Wi have been extensively used. 
The fluctuations of the Ricci tensors are then 

6,Ro^ = [n' + 2n^]Qi-^V^Qi-^V^Wl 

6.Rij = - ^ I {diQj + djQiY + 2n{diQj + 9, Qi) } - ^ (d^W^' + djW'l) 

- n{d^w'^^d,w'i)-\{diWj^d,Wi){2'H! ^m"), (a.i9) 



A. 3 The Tensor modes 

Consider now the case of the tensor modes of the geometry, i.e., according to Eq. (|2.15p . 
the two polarization of the graviton: 

Stgij = -a^hij, 6tg'^ = -j. (A.20) 
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From Eq. ()A.1|) the tensor contribution to the fluctuation of the connections can be ex- 
pressed as 



-0 _ }i(i.> 

2 

StT^j = ^[djh'l + dih] - d%j]. (A.21) 
Inserting these results into the perturbed expressions of the Ricci tensors it is easy to obtain: 

^^t^ = \[hlj + 2nhij + 2{n' + 2n'^)hij - v'^Kj], (A.22) 

5tRl = -^[h^' + mh( - V^/.!]. (A.23) 
A. 4 Fluctuations of the energy momentum tensor 

The energy-momentum tensors usually discussed in the theory of cosmological perturba- 
tions range from the fluid case to the case of scalar fields and possibly vector fields (i.e. 
electromagnetic energy- momentum tensor). The background energy-momentum tensor of 
a perfect fiuid is 

T^,u = ip + p)Uf,Ui, - pg^y. (A.24) 

From the normalization condition Qfj.uU^u'^ = 1 it can be concluded uq = a and 5u^ = —(j)/a. 
Hence, 

6sT^ = 6p, JsTOO = ^{Sp- 2p(t>), (A.25) 



and 



= -5p6{, (A.26) 

d,T^ = 4 \SpS'^ + 2p('iP6'^ - d'd^E)] , (A.27) 

6,T^ = {p + p)v\ (A.28) 

5,ro^ = h(p + py + d'B]. (A.29) 

where we have chosen to define 6u^ = v^/a. Notice that the velocity field can be split into 
divergenceless and divergencefull parts, i.e. 

v' = d'v + V\ d^V' = (A.30) 

Notice that in the case of neutrinos (and possibly also in the case of other collisionless 
species) the anisotroopic stress has to be considered both in the Einstein equations and in 
the covariant conservation equations. The anisotropic stress is the introduced as 

5sTi = -6p6{+-ni, (A.31) 
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where Uj = T/ - 6{t^/3. In the case of scalar fluctuations it is practical to adopt the 
following notation 

-(p + /9)VV, (A.32) 



which is equaivalent, in Fourier space, to the following identity 



{p + p)a = kjk 



(A.33) 



There could be a potential confusion in the definition of the perturbed velocity field. It is 
plausible to define the peculiar velocity in two slightly different ways, namely 



SsTq = {p + p)uo6u' = {p + p)v'' = ip + p)d''v 
SsT° = {p + p)vP5ui = {p + p)vi = {p + p)diV, 



(A.34) 
(A.35) 



where, from the normalization condition g^^u^^u'^ = 1, uq = a and = 1/a. The velocity 
fields V and v defined in Eqs. (|A.34|) and (|A.35|1 are not equivalent. In fact 5ui = avi and 
6u^ = ja. Recall, now, that 



<5n^ = %*-n„)^55*%o + 5'''J^^fc. 



(A.36) 



Inserting now the explicit definitions of (5n* and bu}^ in terms of and and recalling that 
(5gr*o = —d^Bjc?, we have that 

Vi = -Vi - diB. (A.37) 

This difference in the definition of the velocity field reflects in the gauge transformations 
which are different for Vi and Uj, namely 



Vi ^ Vi = Vi + die' , 
Vi ^Vi = Vi - dieo- 



(A.38) 



In this paper we always define the velocity field as in Eq. HA.34|) however, this remark 
should be borne in mind since different authors use different definitions which may only be 
equivalent, up to a sign, in specific gauges (like the class of gauges where B = 0). 

Let us now consider the fluctuations of the energy-momentum tensor of a scalar field (p 
characterized by a potential W{if): 



1 



(A.39) 



Denoting with x the first-order fluctuation of the scalar field 99 we will have 



and , in explicit terms 



1 



,(A.40) 



<5sroo = xV + 2a^^W + a- 



dip 
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6,Toi = if'dix + a^diB 



2a2 



W 



•P 

2a? 



w 



didjE. 



Recalling that 



/2 

Too = ^ + a^W, 



/2 



(A.41) 



(A.42) 



the perturbed components of the scalar field energy-momentum tensor with mixed (one 
covariant the other controvariant) indices can be written as 



i.e., in explicit terms, 



r) tO - — 



6.T^ = ^ 



.'2 



dW 



,2 dW 



+ 



dip 



- --2 



1 / 



,/2 



(A.43) 

(A.44) 
(A.45) 
(A.46) 



The covariant conservation of the energy- momentum tensor implies, in the scalar field 
case, the validity of the Klein-Gordon equation which can be written as 



From Eq. ()A.47|) the perturbed Klein-Gordon equation can be written as 
— — d^W 



(A.47) 



(A.48) 



which becomes, in explicit form, 



dip 



X" + 2nx' - y\ + -Q^^\ + 20— a' - 'f'{4>' + 3V') + 'p'^He' - B) = 0. (A.49) 



A. 5 Generalized scalar perturbation equations without gauge fixing 

Different gauges are suitable for different problems. Hence, it is practical to collect here 
the generalized evolution equations of the scalar fluctuations obtained without fixing a 
particular coordinate system. The sources will be assumed to be barotropic fluids but this 
is not a severe limitation since, with the appropriate identifications, this set of equations 
can even describe more general situations. 
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From Eqs. HA.14|) . HA.25|) and from Eqs. HA.16|) . (|A.28|) the general expressions for the 
Hamiltonian and momentum constraints reads 

VV - nV^{B - E') - Sniip' + n^) = AirGa^Sp, (A.50) 
d'[{^' + n^p) + (W^ - n')B] = -inGa'^ip + py, (A.51) 

From Eqs. (|A.15|) and ()A.27|) the generahzed expression of the (ij) component of the 
perturbed equations is: 

+ (n^ + 2n')<P + n{<l}' + 2^l:')]5l + \5lV\<P -^) + {B- E')' + 2n{B - E')] 
-^did^Kcl) -ij) + {B- E'Y + 2n{B - E')] = 4TrGa^[6p6l - U^]. (A.52) 
Separating, in Eq. HA.52|) the tracefull from the traceless part, we obtain, respectively, 

V'" + n{(j)' + 2iIj') + {Ti^ + 2n')<i)+ \v\(l)-i^) + {B-E')' + 2n{B - E')] = AnGa^Sp. (A.53) 
and 

V2[(0 - V) + (5 - E'Y + 2n{B - E')] = UnGa^ip + p)a. (A.54) 

In general terms the the (0) and (i) components of Eq. HA.6|1 can be written, for scalar 
fluctuations as 

+(2r°o + r^o)4r°o + r^-^T*^' + 4r°r*^' = o, (A.55) 

+(rJo + Tl)6,T^^ + 4r^oo^°° + '^sFi^T'"^ + 2r^,4T0'= = 0. (A.56) 

Inserting now the specific form of the perturbed connections of Eqs. HA.10|) into Eqs. ()A.55|) 
and (|A.56|) the following result can be, respectively, obtained: 

5p' - 3V'(p + p) + ip + p)0 + m{5p + 5p) + (p + p)V^E' = 0, (A.57) 

for the (0) component, and 

(p + p)9' + 9[{p' + p') + m{p + p)] + {p + p)V^B', 

+ \p' + n{p + p)]V^B + V^Sp + {p + p)V^(l) = 0, (A.58) 

for the (i) component. In the above equations, as explained in the text, the divergence of 
the velocity field, i.e. 6 = diV^ = did^v, has been directly introduced. In the text, several 
gauge-dependent discussions are present. The longitudinal gauge equations are obtained, for 
instance, by setting everywhere in the above equations E = B = 0. The off-diagonal gauge 
equations are obtained from the above equations by setting everywhere E = ^ = 0. The 
sychronous gauge equations can be derived by setting everywhere in the above equations 
(j) = B = 0, and so on for the gauge that is mostly suitable in a given calculation. 
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